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Resumen 



Existen numerosas evidencias de la presencia de polvo en galaxias activas. Al 
estudiar la naturaleza del niicleo activo galactico y las condiciones fisicas del medio 
interestelar de estas galaxias es necesario entender como el polvo interacciona con 
la radiacion y las particulas del gas. Si no tenemos en cuenta en nuestro analisis 
los efectos derivados, posiblemente llegaremos a interpretaciones erroneas. Existen 
muchas preguntas relacionadas con la existencia de polvo en galaxias activas como: 
^En que condiciones permite el continuo extremadamente duro del micleo activo 
la supervivencia de los granos de polvo? ^Es la naturaleza del polvo en galaxias 
activas la misma que en nuestro medio interestelar? /,C6mo se distribuye el polvo 
en relation con el gas altamente ionizado por el AGN central? ^Existe polvo en 
radio galaxias con desplazamientos al rojo muy altos? El trabajo desarrollado 
en esta tesis trata de encontrar respuesta a algunas de estas preguntas, mediante 
un estudio detallado, tanto observacional como teorico, de los mecanismos que 
controlan la interaction entre el polvo, la radiacion y las particulas del gas. Los 
efectos observables del polvo sobre el espectro de lineas de emision son tambien 
analizados en profundidad. El objetivo es entender cuestiones mas generales como 
el origen del gas emisor, el/los mecanismos de ionization, la geometria, la conexion 
entre galaxias activas a alto y bajo redshift o la validez del modelo de unification. 
En pocas palabras, este trabajo intenta alcanzar una vision mas clara del mundo 
de las galaxias activas. 
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Summary 



There are strong evidences which favour the existence of dust in active galax- 
ies. Understanding the way in which dust interacts with the radiation and influ- 
ences the physical conditions of the gas is crucial if we want to learn about the 
nature of the central active nucleus and about the physical conditions of the ISM 
in such galaxies. Not taking into account such effects may lead us towards mis- 
leading interpretations. Many intriguing questions concerns to the nature and the 
existence of dust in active galaxies: for instance, under which conditions does the 
very hard ionizing continuum of an AGN allows the survival of dust grains? Is the 
composition and size distribution of the dust the same as in our local interstellar 
medium? How is dust distributed compared to the gas which is at least in part 
highly ionized by the central AGN? Does dust also exist in radio galaxies at very 
high redshifts? The work developed in this thesis tries to find answers to some of 
these questions, through a detailed theoretical and observational research of the 
mechanisms which control the interaction of dust with the radiation and with the 
ions. The observable effects of the dust on the emission line spectrum are also ana- 
lyzed in detail. The final goal has been to give clues about more general questions: 
origin of the emitting gas, ionization mechanisms, geometry, connection between 
low and high redshift active galaxies or the validity of the unification scenario. 
This thesis tries, in summary, to provide a clearer understanding of active galaxies 
in general. 
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Chapter 1 



Introduction 



1.1 The importance of studying dust in active 
galaxies. 

Dust is an important constituent of active galaxies and must be included in any 
physical description of them. It is necessary to reach a clear understanding of its 
influence, not only on the radiation, but also on the physical conditions of dusty 
nebulae. Ignoring dust in the interpretation of the observations is an oversimpli- 
fication which may lead us to false conclusions. An example will describe clearly 
the kind of misunderstandings which can be reached when skipping dust in our 
interpretations. 

High z radio galaxies (z>0.5) show very blue colors and emit a very strong 
Lya line. The initial interpretation has been that these objects were galaxies in 
the process of formation and that the blue young stars could explain the observed 
properties. In this scenario, the continuum is produced by the stars while Lya 
results from recombination of gas from HII nebulosities. Further investigations 
revealed that the continuum was extended and closely aligned with the radio axis. 
This indicated a common mechanism for both the radio and the UV continuum 
emission. Polarization studies (e.g. Tadhunter, Fosbury & di Serego Alighieri 
1988; Cimatti et al. 1993) revealed that the continuum was highly polarized and 
demonstrated that the extended radiation could be scattered nuclear radiation by 
dust within the interstellar medium. This implies that the star formation rate is 
lower than previously believed, also lower than we thought if these objects were 
protogalaxies. There is no clear evidence that these objects are galaxies being 
born, not even at redshift ~ 4. Some evidences rather demonstrate the existence 
of an old stellar population (this will be treated in detail in Chapter 6). The 
redshift corresponding to galaxy formation seems to be much further back than 
the currently highest redshift quasar known at z — 4.9 (Shneider, Schmidt & Gunn 
1991). 

The regions responsible for the emission of the narrow lines in active galaxies 
might well be mixed with dust resulting in selective absorption of the impinging 
ionizing photons and, therefore, affecting the ionization structure of the gas. The 
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physical conditions of the gas are also directly influenced: the temperature and the 
(gas phase) abundances of the heavy elements are sensitive to the presence dust. 
All these effects will manifest themselves in the emission line spectrum. 

In this work we have developed a complete and detailed study (observational 
and theoretical) of the effects of internal (i.e. mixed with the ionized gas) as well as 
external (i.e. intervening) dust, on the emission line spectra of narrow line active 
galaxies. We recall that the emission lines observed in the spectra of this type of 
objects are not contaminated by a broad component emitted within the so-called 
Broad Line Region, which makes our study much simpler. 

This thesis has consisted of three main projects and attempts to address and 
solve some of the important problems encountered nowadays in the field of active 
galaxies. These often involve the presence of dust: 

1) Effects of dust on the Lya and Balmer decrements in the Narrow Line 
Region of Seyfert 2 galaxies. What is the proper geometry describing the Narrow 
Line Region (NLR)? 

2) Does dust exist as a component which is mixed with the ionized gas observed 
in Extended Emission Line Regions of Radio galaxies? What is the origin and fate 
of this extended gas ? 

3) Effects of dust and of resonant scattering on the UV spectrum of high z 
radio galaxies. Which factor plays a major role on the UV line ratios at high z: 
geometry or internal dust? 

These specific objectives are quite different in terms of the type of object 
considered. In effect, we move from low redshift (z) to very high redshift objects 
(which correspond to much earlier epochs) and consider regions which strongly 
differ by their physical conditions. However, the goal remains the same: to under- 
stand the nature of active galaxies. Many researchers are working on this thrilling 
field, trying to answer many questions which remain open, using different theoret- 
ical and observational techniques. In our work, we have modeled the emission line 
spectra of the regions under study. Starting from a given set of input parameters 
(some of them well established by the observations), we create theoretical spectra 
which are then compared with the observations. The discrepancies or agreements 
resulting from the comparison provide important information about the nature of 
the regions and objects under study. 



1. INTRODUCTION 



5 



1.2 Evidences of dust in active galaxies. 
1.2.1 Absorption 

Extinction occurs whenever electromagnetic radiation is propagated through a 
medium containing small dust particles. In general, the transmitted beam is re- 
duced in intensity by two physical processes, absorption and scattering. The energy 
of an absorbed photon is converted into internal energy of the particle, which is 
thus heated, while a scattered photon is deflected from the line of sight. There is 
direct evidence for dust in our Milky Way in the form of dark lanes and rifts that 
we can see with the unaided eye. Also in some active galaxies, we can detect the 
presence of dust from the existence of obscured regions, like Cen A (NGC5128), the 
nearest giant double radio galaxy. In 1847, Sir John Herschel observed for the first 
time the nowadays familiar dark lane contrasted against the diffuse background of 
light. Currently, it is known that the galaxy consists of an ellipsoidal stellar com- 
ponent bisected by a dense dusty layer across the apparent minor axis (see Fig. 1) 
(e.g. Bland et al. 1987; Morganti et al. 1991). Recent HST WFPC2 images of 
3CR radio galaxies in the redshift range 0.0< z <0.5, show lanes, patches and 
wisps that the authors attribute to dust out to a redshift of ~0.48 in ~ 30-50% of 
the sample (Baum et al. 1995). At z >0.1, dust is more difficult to detect and one 
cannot easily separate dust obscuration from emission dumpiness. 



Figure 1: 

The dusty lane in the active galaxy Centaurus A. A deep (1 hr) image taken at the prime 
focus of the Anglo-Australian Telescope. (Reproduced from Morganti et al. 1991) 

Selective extinction (or reddening) is also evident in the spectra of Narrow 
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Line Regions (NLR) in Seyfert 2 galaxies. In effect, the Lya/H/3 and Ha/H/3 
ratios indicate a mixture of reddening due to dust mixed with the emitting clouds 
and of resonant scattering of Lya by neutral hydrogen (Binette et al. 1993, see 
also Chapter 2) (some authors propose reddening by pure dust (i.e. gas-free dust!), 
like Ferland and Osterbrock 1986, and Wills et al. 1993). 

The narrow emission line profiles provide another indication of dust in the 
NLRs. Many lines have a noticeable blue asymmetry which is interpreted as a 
combination of radial motion and dust obscuration (e.g. de Robertis & Osterbrock 
1984). The asymmetry parameter (defined by Heckman, Miley & Breugel 1981) 
measures the asymmetry of the line profiles. Dahari & de Robertis (1988) searched 
for correlations between this parameter and features attributed to dust such as 
reddenning (meaasured from the Balmer decrement, Ha/H/3) far IR luminosity 
(60-100 fim) and the flux ratio F(60/mi)/F(A5007). The results showed a good 
correlation in Seyfert 2 galaxies between the asymmetry parameter and these three 
variables. The authors propose that such correlations suggest a direct relation 
between the asymmetry parameter of the lines and the presence of dust in the 
NLR. 

1.2.2 Scattering and Polarization 

The existence of scattered light from the interstellar dust grains was discovered 
by means of long exposure photographs of systems such as the Pleiades cluster of 
stars. The associated nebulosity is referred to as reflection nebulosity. The light 
of the stars in the cluster illuminates the dust grains which then scatter the light, 
most efficiently at blue wavelengths, explaining the blue colors associated with the 
nebulosity. The scattering of radiation has a useful observational consequence: the 
scattered light is polarized. This was first demonstrated by Hall (1949) and Hiltner 
(1949) who showed that the light of reddened stars is partially plane polarized, 
typically at the 1-5% level. 

Polarization has often been observed in active galaxies and the interpretation 
generally proposed has been scattering of the radiation by an agent which nature 
does not reveal clearly: either electrons or dust grains. Antonucci and Miller (1985) 
discovered that in the Seyfert 2 galaxy NGC1068, both the nonstellar nuclear 
continuum and the weak broad wings to the Balmer lines and the Fell emission 
are linearly polarized with the polarization plane perpendicular to the symmetry 
axis of the nuclear radio morphology. The most plausible interpretation is the 
existence of a Seyfert 1 type nucleus and of a Broad Line Region hidden from our 
line of sight by an obscuring torus. The continuum and broad lines which are 
detected in polarized light, are emitted originally within these hidden regions and 
later scattered towards the observer. The most plausible scattering agent is, in 
this case, electrons, although dust could also play a role. 

Since Antonucci and Miller's discovery, the number of Seyfert 2 galaxies that 
harbor hidden BLRs visible only in polarized flux spectra has noticeably increased: 
Mrk 3, Mrk 48, NGC 7212, etc (e.g. Miller & Goodrich 1990; Tran 1995). The 
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nature of the scattering material is not clear, dust and/or electrons? But, at least 
for some objects, scattering by dust seems to be the dominant mechanism (e.g. 
NGC 7674, Tran 1995). Recently, Inglis et al.(1995) have mapped the polarization 
of NGC1068 also in regions out of the nucleus. In off-nuclear regions, dust is also 
responsible for the scattering. 

The low redshift radio galaxy PKS2 152-69 (see Fig. 2) shows extra nuclear 
scattering light too. It contains a very blue cloud at 8 kpc from the nucleus which 
is aligned with the radio axis and emits highly polarized flux. Dust (and electrons) 
in this cloud could be the scattering agent of the central UV radiation emitted by 
the nuclear source. 



Figure 2: 

A WFPC2 UV-continuum (F300W) image of PKS2152-69. The angular separation be- 
tween the nucleus (lower right) and the jet-cloud interaction front is about 8 arcsec. 
The interaction site emits X-rays (ROSAT HRI), optical/UV emission lines, a blue/UV 
continuum as well as radio waves. (Fosbury et al. in prep) 

Scattering and polarization (maybe by dust) are also detected at very high 
redshifts. High redshift radio galaxies (z>0.5) present elongated structures in the 
UV rest frame, which are aligned with the radio axis. This extended UV radiation 
is polarized, with the electric vector perpendicular to the axis of the optical struc- 
tures. The observations indicate that the radiation emitted by a hidden quasar is 
being scattered by dust and/or cool electrons in the extended ISM. The discovery 
in several radio galaxies (e.g. di Serego et al. 1996; Dey & Spinrad 1995) of broad 
polarized MgIIA2798 emission line (originally emitted in the Broad Line Region 
which is also hidden from our line of sight), supports this scenario. 

1.2.3 Thermal emission of warm dust 

Diffuse emission from interstellar dust was predicted by van de Hulst (1946) as 
a consequence of absorption: the energy absorbed by dust grains must re-emerge 
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in the infrared. In a typical interstellar environment, a dust particle gains energy 
mainly from the absorption of ultraviolet photons from the ambient interstellar 
radiation field. The dust grains will emit in the IR a power equal to that absorbed. 
These effects are also clearly seen in active and starburst galaxies. 

Within the framework of unification theories of AGNs (see Antonucci 1993 for 
a review), it is believed that a dusty torus surrounds the central accretion disk and 
the broad line region of active galactic nuclei (AGNs). The main continuum source 
and the broad lines would be obscured from direct view by a torus in narrow line 
active galaxies. The dust in the torus should absorb much of the nuclear radiation 
and re-emit it in the infrared. Measurements by Rieke & Low (1975) and Stein & 
Weedman (1976) showed already in the 70 's that many Seyfert galaxies have large 
infrared excesses. In fact, the overall characteristics of the spectra of AGN show a 
certain gross uniformity with much of the energy emitted in the IR. This emission 
is attributed to warm dust, heated either by the central AGN or by recently formed 
stars, or both. At very high z, the dust thermal emission is redshifted to the far 
IR and the submm waveband, depending on z. Several high redshift radio galaxies 
have been detected with excesses in these spectral ranges, which indicates the 
existence of warm dust also at much earlier epochs of the Universe (e.g. Downes 
et al. 1992; Isaak et al. 1994). 

The unified models of AGNs propose that the diversity of AGNs is smaller 
than believed and that many of the observed differences are due to orientation 
effects: the obscuring torus determines a radically different appearance of the 
AGN at different viewing angles, AGNs of different orientation will be assigned to 
different classes. This scenario has been tested by studying dust emission from the 
torus (Hess 1995): for instance, if radio galaxies and quasars are intrinsically the 
same objects, the torus should intercept the same amount of ultra-violet photons 
in both types of objects and we would expect that they emit similar output of 
reprocessed infrared emission. IRAS data indicate that quasars are brighter at 60 
jum than narrow line radio galaxies, which contradicts the unified models. This 
could be explained if the obscuring circumnuclear torus is still optically thick at 
60/zm. Other mechanisms dependent on orientation could explain the differences 
if for instance the IR emission was not only due to thermal dust emission but also 
due to a non thermal component. In this case, relativistic beaming of the core 
component might influence the 60 \i emission. 

1.2.4 Spectral absorption features 

Interstellar absorption features in the spectra of reddened stars and infrared sources 
are attributed to solid particles on the basis of position, width, shape and continuity 
of profile. The absorption lines due to gas phase atomic or molecular species in the 
interstellar medium are, in general, extremely sharp (FWHM of 1 km s _1 or less), 
reflecting the conditions of low temperature and pressure in which they are formed. 
Solid state spectral features are intrinsically broad and continuous and cannot be 
resolved into discrete lines, in contrast to the vibration-rotation bands of many 
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gas phase molecules. Of the observed absorption features attributed to interstellar 
dust, the A1275 feature in the mid-ultraviolet is the strongest and, in frequency 
units, the broadest, forming a prominent peak in the interstellar extinction curve. 
In principle, dust-related spectral features provide a direct means of identifying 
the chemical composition of interstellar grains. 

Active galaxies also show absorption features due to dust. Mrk 231 is a Seyfert 
1 galaxy with very strong IR emission. Its optical spectrum is very heavily red- 
dened. Relatively high resolution infrared measurements by Rieke (1976) show the 
redshifted interstellar absorption feature due to silicates at 10 /im, which means 
absorption certainly arises in Mrk 231. 

However, the absence of spectral features which are customarily attributed to 
dust in our galaxy should not be interpreted necessarily as indicating the absence 
of dust. The dust in external galaxies might not have exactly the same extinction 
properties as in our interstellar medium. Even in our own Galaxy, there are regional 
variations in the optical properties of the interstellar dust. Already in 1937, Baade 
and Minkowsky found that the extinction curves for stars in the Orion Nebula 
differ from the "mean" galactic extinction curve. Selective removal of the small 
particles is the reason. Environmental influences (grain growth by coagulation, 
size dependent destruction, etc) can strongly influence the composition and size 
distribution of the dust grains and, therefore, the extinction properties. Dust in 
active galaxies may exist under very different conditions than in the ISM of our 
galaxy. 

1.2.5 Spectral emission features 

The principal emission signatures attributed to cosmic dust arise in silicates and 
hydrocarbons. The PAH spectrum is a series of narrow emission features with 
principal wavelengths at 3.3, 6.2, 7.7, 8.6 and 11.3 fim. They are now widely 
attributed to polycyclic aromatic hydrocarbons, a class of organic molecules com- 
posed of benzene rings. Their infrared spectra are characterized by a number of 
resonances, excited by absorption of ultraviolet photons (Duley and Williams 1981; 
Leger and Puget 1984). The PAH features are commonly observed in the lines of 
sight to planetary nebulae, HII regions and reflection nebulae. 

PAH features are observed in the nuclei of many external galaxies. They are 
particularly prominent in systems with nuclear HII regions and are therefore in- 
dicative of recent OB star formation in which an intense ultraviolet radiation field 
very soft beyond the Lyman limit can excite the different transitions. Indeed, the 
occurrence of PAH features appears to be a useful discriminator between starburst 
galaxies (dominated by rapid star formation) and active galaxies (dominated by 
compact, energetic nuclei). PAH emission occurs in the spectra of starburst galax- 
ies because many luminous HII region are included in a resolution element. In 
contrast, PAH features are generally absent from the spectra of active galaxies, 
suggesting that the carriers are destroyed by hard UV radiation in the vicinity of 
the nuclear source (Aitken & Roche 1985; Desert & Dennnefeld 1988). 
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1.2.6 Depletion of heavy elements 

The term depletion refers to the underabundance of gas phase elements with res- 
pect to the solar standard as a result of their presence in the dust grains. Ultraviolet 
spectroscopy of interstellar gas indicates a depletion in the abundances of many of 
the heavy elements with respect to solar values, which is most readily explained if 
the missing atoms are tied up in solid particles. 

Ferland (1993) demonstrated the presence of dust grains in the narrow line 
region gas. He showed that in the absence of dust the [Call] near infrared forbidden 
lines (AA7291,7324) would be among the strongest lines in the optical to infrared 
spectrum of the NLR, assuming cosmic Ca abundance. He interpreted the non 
detection of these lines as due to the depletion of calcium in the NLR, showing 
that grains must be present, mixed with the ionized gas. 
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1.3 The projects developed in this thesis 

1.3.1 Effects of dust on the Lya and Balmer decrements in 
the Narrow Line Region of Seyfert 2 galaxies. (Chap- 
ter 2) 

Despite its simple atomic configuration, the hydrogen line spectrum of AGNs is 
not well understood. In 1984, Gaskell and Ferland showed that the theoretical 
spectra of many NLR models of Seyfert 2 define a narrow band in the Lya/H/3 vs. 
Ha/H/3 diagnostic diagram named the intrinsic band. The displacement inside this 
band is explained by variations in metallicity, hardness of the ionizing continuum 
and electronic density. However, the observed data locate many Seyfert 2 galaxies 
outside the theoretical band, sometimes quite far away from it, being Lya much 
fainter than expected with respect Ha and R(3. Previous studies (e.g. Ferland & 
Osterbrock 1986; Wills et al. 1993) propose a combination of atomic processes and 
dust absorption to explain the observations. 

In this work we propose a new scenario that solve the discrepancies between 
the HI observed and predicted ratios of the NLR, a scenario where dust plays an 
important role, but also geometry, a factor which has not been considered before 
and which strongly influences the line emission, specially the resonant line Lya. 

1.3.2 Does dust exist within the ionized gas of the Ex- 
tended Emission Line Regions of Radio galaxies? ( Chap- 
ters 3 and 4 ) 

The main goal of this project has been to investigate the origin of the gas which 
forms the extended emission line regions (EELRs) and which is observed to extend 
up to several hundreds of kpc in powerful radiogalaxies. In some cases, morpho- 
logical similarities with filamentary systems surrounding galaxies near the centers 
of rich clusters suggest that the origin of the EELR consist of gas cooling from 
the hotter phase, similar to the X-ray emitting corona which surround many cD 
galaxies in clusters. In other cases, morphological and kinematical studies of the 
gas often suggest the existence of a collision or merger in the recent past. 

The existence or not of dust in such regions has important consequences for 
the origin of the gas: if the gas had cooled from a hotter X-ray phase, any dust 
introduced in the intracluster medium would have been sputtered and rapidly 
destroyed. In the cooling process there would have been not enough time for 
the dust to form. On the other hand, if the gas consists of galactic debris, the 
gas/dust ratio is expected to have a value appropriate to the chemical composition 
of a normal galaxy. 

We have investigated the existence of dust in the EELRs of a sample of low 
z radio galaxies in order to discriminate between these two main theories relative 
to the origin of the gas. 
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1.3.3 Effects of dust and resonant scattering on the UV 
spectrum of high z radio galaxies. (Chapter 5) 

Radio galaxies at very high z (z >2) have been until very recently the youngest 
galaxies we know of and, therefore, good laboratories to study galaxies at early 
epochs. In order to determine correctly the epoch when they were formed, we need 
to establish the age of the oldest stars in such systems. However, we don't have 
yet a clear picture of the processes which are going on in these special objects. 
The presence of an active (very powerful) nucleus influences strongly the ISM 
of the galaxy and it is difficult to separate the different components (including 
stars) which contribute to the observed radiation. We need, therefore, a clear 
understanding of the physical processes involved in the formation of the various 
lines and of the continuum in order to be able to disentangle the stellar from the 
AGN-related sources. 

In our work we interpret the UV (rest-frame) emission line spectrum of a 
sample of high redshift radio galaxies. Due to the short wavelengths and the 
presence of several resonant lines, this spectral range is very sensitive to geometry 
and dust. Comparing our photoionization models with the observed line ratios, 
we study in detail the effects that both, dust and geometry have on the UV line 
ratios. Our study provides results about the presence of dust and the 3-dimensional 
structure of these galaxies, which can be very useful for understanding the nature 
of such extreme objects. 



Figure 3: 

A composite radio galaxy spectrum in the UV rest-frame constructed from observations 
of 0.1< z <3 galaxies (from McCarthy 1993). 
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Effects of internal dust on the 
Narrow-Line region Lyman and 
Balmer decrements. 

Binette, Wang, Villar- Martin, Martin & Magris 1993, ApJ, 414, 535 

Abstract 

We present detailed calculations on the effects of internal dust on the Balmer 
and Lyman decrements for a spherically symmetric distribution of low covering 
factor clouds represented as slabs and consider the effects of Lja resonance and 
absorption by dust. We consider the important effects of perspective on the emer- 
gent fluxes, which in our simplified scheme present either the photoionized face to 
the observer ("f") or the back ("b") face. We adopt canonical values for the gas 
excitation (Uf) and for the ionizing energy distribution (a = —1.4, F„ ~ u +a ) and 
compute sequences of photoionization models in which the relative internal dust 
content by mass (/i) is progressively increased to values comparable to the local 
ISM (i.e., /i=l). We find that for moderate amounts of internal dust /i=0.2-0.3, 
radiation-bounded clouds result in Lyman and Balmer decrements in the range 
32-37 and 3.0-3.1, respectively. Our main result, valid for the calculations with 
standard NLR input parameters and for an open geometry in which multiple cloud 
covering is negligible, is that even with internal dust, while Ha/H/3 is reddened, 
Lya/H/3 turns out not very different from recombination case B. The Seyfert 2 
observations show much lower Lyman decrements than predicted from the spheri- 
cally symmetric model and a critical study is made of various explanations: dust 
-or line scattering screens, a semiopen geometry can span the region covered by 
the Seyfert 2 observations. We also study the possibility that dust reddening of 
the continuum may account for the apparent deficit of ionizing photons seen in 
many Seyfert 2's. This reddening may be present in addition to, or in place of, 
the anisotropic beaming/occultation of ionized radiation (cf. "occultation/reflection 
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picture") that is generally invoked to explain the deficit. 

2.1 Introduction 

The calculation of the intrinsic Balmer decrement emitted by the narrow-line region 
(NLR) of Seyferts has attracted much interest in the past (e.g., Gaskell 1984; 
Malkan 1983; Halpern & Steiner 1983; and Ferland & Netzer 1983) because it 
is the most direct means of determining how much (intervening) extinction the 
observed emission-line spectrum has undergone. This information can be used to 
deredden the whole line spectrum for the effect of external dust along the line of 
sight. It has been recognized early on by Halpern & Steiner (1983) and Ferland 
& Netzer (1983) that the NLR intrinsic Balmer decrement was likely to be higher 
than that of recombination case B as a result of collisional excitation within the 
partially ionized zone. The resolution of the problem of intrinsic HI line ratios was 
further advanced by combining the information provided by the Lyman decrement 
with that of the Balmer decrement (e.g., Ferland & Osterbrock 1985). Gaskell & 
Ferland (1984) studied in details how the NLR metallicity, density, and hardness of 
the ionizing continuum affected both decrements. The most probable range for the 
NLR decrements which these authors favor is Ha/H/?=2. 8-3.1 and Lya/H / 9=30- 
50 (similar to values determined earlier by Ferland & Osterbrock 1985). The 
observational determination of the intrinsic HI line ratios is facilitated in Seyfert 
2 galaxies since one need not decompose the line profiles into BLR and NLR 
components. However, Seyfert 2's are on average significantly more reddened than 
Seyfert l's (Gaskell 1984). De Zotti & Gaskell (1985) concluded that the dust 
was probably associated with the NLR clouds or filaments owing to the extremely 
weak correlation of the extinction with disk inclination angle and also due to the 
significantly lower extinction found in Seyfert 1.5. 

If there is so much dust associated with the nuclear emission region, we might 
expect some dust to be associated with the photoionized gas itself. The effect 
of internal dust mixed with ionized gas for physical conditions pertaining to the 
broad-line region (BLR) has been investigated by Ferland & Netzer (1979). More 
recently, however, Netzer (1993) concluded that the BLR was devoid of dust (but 
see Crosas & Weisheit 1993). From considering the ratios Ha/H/3 vs. H7/H/3 in 
Seyfert 2's, Binette et al. (1990) suggested the possibility of an intrinsic (NLR) 
Balmer decrement as high as Hct/H/3 m 3.4, which they modeled using photoion- 
ization calculations with dust internally mixed with the ionized gas. 

Internal dust results not only in somewhat higher Balmer decrements but can 
furthermore increase the equilibrium temperature of the ionized gas. This property 
was used by Magris, Binette, & Martin (1993) to solve the problem of the high 
electronic temperature (based on [OIII] lines) observed in the extended ionized 
nebulosities of powerful radio galaxies (cf. Tadhunter, Robinson & Morganti 1989). 
Furthermore dust can reflect towards the observer a nonnegligible fraction of the 
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impinging optical radiation, an interesting property which appears to explain quite 
well the extended blue polarized continuum observed in distant radio galaxies (cf. 
Fosbury 1993; Cimatti et al. 1993; Binette et al. 1993a). Dust may also explain 
the absence of any intermediate region between the BLR and the NLR according 
to Netzer & Laor (1993), not to mention its relation with the infrared emission 
(Clavel, Wamsteker, & Glass 1989; Barvainis 1992; Sanders et al. 1989; Pier & 
Krolik 1992). 

In this work , we expand upon the calculations of Binette et al. (1990) and 
study the effects of internal dust on the intrinsic Hct/H/3 by taking into account the 
observer's perspective of the emitting clouds. In addition, we consider the effects 
of dust on the Lyman decrement. Kwan & Krolik (1981) and Puetter & Hubbard 
(1987) have shown using dust-free BLR calculations how significant is the effect 
of perspective on resonant Lya. With dust mixed in, it becomes essential to take 
into account satisfactorily the cloud's perspective and the global NLR geometry 
when deriving the intrinsic decrements, as described in §2.2. We first investigate in 
§2.3 the magnitude of the Lya destruction due to internal dust in photoionization 
calculations and compare in §2.4 our calculations with the Seyfert 2 data of Kinney 
et al. (1991) (hereafter KAW3). 

2.2 Computational method 

We describe the photoionization code and the line transfer of Lya, Ha, and Wf3 
across the dusty medium and review the input parameters that enter our calcula- 
tions. 

2.2.1 The photoionization code 

To compute HI lines we have employed the multipurpose photoionization-shock 
code MAPPINGS (cf. Binette, Dopita, & Tuothy 1985). To compute the hydrogen 
lines, we treat the hydrogen atom clS db SIX - level system (Is, 2s, 2p, 3, 4 and 5) 
plus continuum. The hydrogen levels are populated through direct recombination 
(followed by cascade) as well as by collisional excitation either by thermal electrons 
(using Johnson 1972 rate coefficients) or by suprathermal electrons (adopting the 
method developed by Shull & Van Steenberg 1985) as earlier described in Binette 
et al. (1993b, hereafter BWZM). 

In the computation of the ionization structure, we adopt a simple-plane- 
parallel geometry whereby our putative thick gas cloud is represented by a radiation- 
bounded slab of gas containing interstellar dust. The ionized slab is subdivided 
into many small layers in which the equations of photoionization and thermal 
equilibrium are solved by standard methods. 
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The escape probability formalism is used to solve for the transfer of the reso- 
nance lines. The effects of dust on the ionization structure as well as on the thermal 
balance of the plasma are considered in detail by implementing the physical pro- 
cesses relevant to dust mixed with emission plasma as described in Appendix C of 
Baldwin et al. (1991). For instance, heating of the plasma by dust photoelectric 
emission is considered. The dust grain charge is calculated self-consistently and 
the formula describing the photoelectron energy distribution and the yield are from 
Draine (1978), but with a cap of 0.2 for the yield at high photon energies as in 
Baldwin et al. (1991). 

One interesting aspect of the new code MAPPINGS is that the effect of dust 
scattering on the line transfer is explicitly solved using the numerical solution 
of Bruzual, Magris & Calvet (1988). This aspect of the transfer is particularly 
relevant to the "open geometry" adopted for most of the calculations presented 
in this thesis. Destruction of the resonant line by dust absorption is taken into 
account using the results of Hummer & Kunasz (1980). We describe in more details 
in Appendix A of BWZM the method implemented in MAPPINGS for solving the 
line transfer of emission plasma with internally mixed dust. 

Trace elements and dust content 

Even though this work only addresses the intensities of hydrogen lines, the abun- 
dance of trace elements remains an important factor since they affect directly the 
thermal balance of the photoionized plasma. Abundance gradients observed in 
spiral galaxies point towards a metallicity even higher than solar for the nuclear 
region. Many models of the NLR on the other hand have favored more often than 
not the exploration of abundances lower than or equal to solar (e.g., Gaskell & Fer- 
land 1984; Ferland & Netzer 1983). As this work emphasizes the role of internal 
dust, for defmiteness we will assume that the abundances of metals appropriate 
to the NLR are close to solar and adopt the solar system abundances of Anders 
& Grevesse (1989) (see list in Table 4 of Appendix A). We define Z as the total 
(gas+dust phases) metallicity relative to solar which implies Z=l for this work. 

The absorption and scattering cross sections of dust (from the infrared to the 
soft X-rays) were computed by Peter Martin (the curves are displayed in Fig. 6 
of Martin & Rouleau 1991) and corresponds to a dust models of the extinction 
curve of the solar neighborhood interstellar medium. In MAPPINGS, we scale the 
scattering and absorption cross sections (and, therefore, the dust content) by a 
dimensionless factor p. The grain model of Martin & Rouleau (1991) results in a 
dust-to-gas mass ratio of 0.0064 p (p dust / P gas — 0-7^-Pdust/ Ph) and to an extinction 
opacity at 5500 A r y (=0.92lA y ) = 4.8 x !0~ 22 pN H cm 2 , where iV ff (cm" 2 ) is the 
total hydrogen column density. When p = 1, these quantities are all consistent 
with the standard extinction curve of the local ISM (interstellar medium). 

To be consistent with the depletion of metals unto dust grains, we deplete 
the abundances of metals in the gas phase, Z gas , according to the dust content 
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fi. The adopted scheme described in Appendix A uses the depletion indices listed 
in Whittet (1992) and if necessary scales them with /x/Z. Renormalization of the 
derived gas phase abundances is performed until the mass locked into dust grains 
becomes consistent with the dust-to-gas mass ratio implied by ji. For instance, 
when ii— 1, Z gas = 0.6 as a result of depletion. In this work, the total metal 
abundances are always solar and, therefore, Z = Z gas + Zdust = 1 

The slab geometry and the input parameters 

All our calculations assume the same standard values for the input parameters, 
that is: an isobaric density behaviour, a power law ionizing continuum of index 
a = —1.4 (F„ oc u +a ) and solar total abundances (but depleted gas abundances: 
Z gas = Z — Z dust = 1 — Z dust ). Our constant pressure models are parameterized 
in terms of the so-called ionization parameter, that is, the ratio of the density of 
ionizing photons impinging on the slab to the density of the outermost gas layer 
of the slab: 




where c is the speed of light and z/ is the Lyman limit frequency, (p[ is the 
ionizing photon flux (in photons cm' 2 s" 1 ) impinging on the slab (quantities with 
superscripts / or b refer to their values at the front or at the back of the slab, 
respectively). For the reference dust-free model (filled star in the figures) which is 
referred below, we adopted n^=5000 cm~ 3 and £//=0.0015 as justified in §§ 2.3.1 
and 2.3.3. 

Using the code MAPPINGS, the calculation of the ionization and emissivity 
structure is carried inward in the slab either up to a column depth (i.e., J N^dx) 
N^ ab < N H * in the case of matter-bounded calculations or, in the case of radiation- 
bounded calculations, up to Nh*- The column density Nh* gives the depth of the 
complete "photoexcited" region, that is, the depth at which the incoming ionizing 
flux is exhausted. The operative definition of the boundary of the photoexcited 
region (at depth Nh*) is defined as the depth where the following two conditions 
are simultaneously satisfied: (1) the inabsorbed ionizing flux <fil < 1 %, of the 
impinging flux <fil , and (2) the ionized fraction n b H+ /n H < 1%. The code computes 
successively the emergent flux seen directly from the phtoionized face (perspective 
"f"; see Fig.4) as well as seen from the back (perspective "b"). The transfer 
solution used across the dusty medium not only considers absorption but dust 
scattering as well (see BWZM). 

Because the harder photons of the power-law ionizing continuum (with a = 
— 1.4) creates quite a large partially ionized region (hereafter PIZ) in which the 
ionized gas very gradually becomes neutral, Nh* is a factor ~ 10 higher than the 
Stromgren depth (in cm~ 2 ) which is technically defined as 
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Figure 4: 

Adopted slab geometry for the constant pressure photoionization calculations. The slab 
comprises a fully ionized zone and a partially ionized zone (PIZ) which adds up to the 
column density of the photoexcited zone Nh* . Beyond the photoexcited regions a zone of 
neutral gas may exist of column density N^ . Perspective "f " considers the emergent line 
fluxes seen from the side which is photoionized (front), while perspective "b" corresponds 
to the opposite back side. The dust extinction opacity of any zone is proportional to its 
relative dust content [i and total hydrogen column density Nh- 

N H+ = cU f /a B w 10 23 L7 / cm" 2 (T e = 10 4 K) [2] 

where as is the recombination coefficient to excited states of hydrogen. In an 
HII region, because of the softness of the ionizing continuum, the PIZ is negligible 
and N H * ~ N H+ ■ For radiation-bounded clouds, equation (2) entails that the fully 
ionized zone [Njj+(tih+ /nuo > 10) ~ 0.8iVj| + ] which contributes most of the line 
luminosity nevertheless contributes only a small fraction of the dust opacity inside 
Nh* and, therefore, the spectrum seen from perspective "b" is more reddened than 
the one emerging from perspective "f". 

By portraying the emitting NLR cloud as a slab, a determinant parameter with 
which the dust opacity scales is the total column density N^1 b (= N^o b + iVj/ ab ) 
of the slab. In this work, in most instances, we assume the slab to be not only 
radiation-bounded but in some cases to exceed Nh* due to a neutral H° zone of 
column density N H0 (beyond the PIZ) as depicted in Figure 4. The dust inside this 
region can play a determinant role in reducing the emergent line fluxes seen from 
perspective "b". Apart from absorption of the lines by dust, the neutral zone also 
is assumed isothermal with a temperature taken to be equal to the temperature at 
the last computed ionized layer. Changing this temperature even arbitrarily has 
practically no effect on the zone's reflectivity to Lya since even without such a 
zone the already huge Lya line opacity within the PIZ would cause most of the 
Lya to emerge from perspective "f" rather than from "b" (cf.2.3.2). 
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In our nomenclature, the total column density of the slab is given by JVjy = 
N H * + N^o (see Fig. 4) if the slab is radiation-bounded (otherwise, N^ ah < N H * 
if the slab is matter-bounded). We assume that dust-to-gas ratio ji is uniform 
within each zone but may take different values in the photoexcited zone from that 
in the neutral zone. Assuming that the same dust size distribution and composition 
applies to the different zones (that is, the shape of the extinction curve remains 
constant), we obtain that the total slab opacity in the V band (5500 A ) is given 
by 

T# ab = T V (N H .)+T V (N° ) = 

= 4.8 x lQ- 22 [u,{H*)N H *) + u,(H°)N° H0 ] [3] 

Due to the high reflectivity of Lya by the H° zone as well as by substantial 
amount of H° inside the transition zone, the emergent Lya originating from per- 
spective "b" is generally negligible. The Balmer lines are generated isotropically 
but are absorbed differently when emerging from perspective "b" or "f" because 
of the asymmetry in dust opacity introduced by the N^ zone and the PIZ. If we 
picture the NLR as an ensemble of clouds with very low covering factor symmetri- 
cally distributed around a nuclear ionizing source, the number of cloud seen from 
the back can be postulated to be equal statistically to that of clouds seen from 
the front. To a first order the resulting line ratios for such an "open" geometry is 
obtained by adding first with equal weight the line fluxed from the "f" and "b" 
perspectives and then taking the ratios of the summed fluxes. 

In § 2.3, we will first explore the effect of varying the internal dust-to-gas ratio 
fi on the Lya/H/3 and Ha/H/9 line ratios emerging from a symmetric distribution 
of clouds (i.e., the case represented by "f" + "b"). After comparing models (§2.4.1) 
with Seyfert 2 data of Kinnety et al. (1991) and after discussing the pure dust 
hypothesis (§2.4.2), we will look into the possibility of having more clouds seen 
from perspective "b" as well as the limiting case of a spherically closed geometry 
(§2.4.3). 

2.3 Results 

Our main purpose is to explore the effect of internal dust on the HI line ratios in 
the context of the standard photoionization model of the NLR. For this reason, 
we have limited the scope of the calculations to restricted plausible values of input 
parameters; we consider, for instance only solar abundances (but with progressive 
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depletion of gas phase metallicity Z gas with increasing dust content u) and canon- 
ical a = —1.4 ionizing power law [] (cf. Ferland & Osterbrock 1986; Kinney et al. 
1991). The choice and role of the gas density is first discussed (§2.3.1) as well as 
that of incomplete opacity to ionizing photons (§ 2.3.2) and the selection of the 
ionization parameter (§2.3.3). We then discuss the results of nebular models with 
internal dust §2.3.4. 

2.3.1 Selection of a characteristic cloud density 

The choice of a representative cloud density in the calculations is much more 
problematic than the selection of other parameters. There is evidence of a wide 
distribution of cloud densities in the NLR following the work of DeRobertis & 
Osterbrock (1984,1986) who have shown in many Seyferts the tendency of the 
forbidden line profile widths to correlate with the respective deexcitation critical 
density of the transition involved. This lends support to the picture of a wide 
gradient in cloud densities (10 3 < n e < 10 6 cm -3 ) as well as in velocity dispersion 
as a function of distance from the nucleus. In this picture, the bulk of a line 
luminosity would be contributed by clouds of density approaching the line's critical 
density. The few available direct measurement of the NLR electron densities have 
been mostly provided by the ratio of the [SII]AA6716,6731 density-sensitive doublet 
(Osterbrock 1989) which show densities which are relatively low (200-2000 cm -3 ). 
These low densities do not necessarily contradict a wide range of cloud densities 
within the NLR because the critical density of [SII] red lines is relatively low and 
moreover the S + zone coincides with the partially ionized zone of the emitting 
cloud where the electron density can be expected to be lower than in the fully 
ionized zone. This can be demonstrated with constant pressure photoionization 
calculations with densities covering; the range 80< n f H < 6000 cm which show 
that the electronic density inferred from the [SII] line ratio is a factor 4-5 smaller 
than the mean < n e > of the fully ionized zone. Finally, there also exists a 
significant fraction of Seyferts in which line widths correlate with line excitation 
(cf. Pelat, Alloin, & Fosbury 1981; Filippenko 1985; DeRobertis & Osterbrock 
1984,1986) which could indicate that the density range within these objects is 
narrower. 

Following the indications given by the [SII] densities and given that we will 
not integrate the line contribution of clouds over different densities, we opt for 
a single density of n^=5000 cm -3 which we expect to characterize the NLR gas 
that emits the bulk of the HI lines. Before discussing models with dust we first 
summarize the dependence of HI lines on density. 

The direct dependence of Ha and H/5 intensities on the gas density is very 
weak when the line is produced by recombination alone (Hummer & Storey 1987). 
However, in a power law photoionized plasma, both lines become indirectly coupled 

1 A much harder distribution (flatter power law) would give a rather poor fit to the low 
excitation lines assuming the photoionized clouds are radiation bounded. 
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to density when substantial collisional excitation takes place within the PIZ. The 
reason is that increasing the density reduces the effectiveness of cooling of many 
important forbidden lines as a result of collisional deexcitation, thereby increasing 
the equilibrium plasma temperature which in turn increases the rate of collisional 
excitation of HI in the PIZ. In the fully ionized zone, there is too little HI to 
excite and the Balmer lines are produced only by recombination independent of the 
plasma temperature. In the case of Lya, the dependence of the line intensity on the 
electron density is much more direct. When the density is increased much above 
10 3 cm -3 , there is an increasing conversion of the level 2s population otherwise 
responsible for the 2q continuum emission into Lya emission by excitation into the 
2p level (see Osterbrock 1989; Gaskell & Ferland 1984). In the low density regime 
case B in the absence of the Lya sinks (e.g., dust absorption), the fraction of 
energy emitted as Lya versus 2q emission from recombination alone is oc^J /ocb ~ 
0.67(T e /10 4 fT)" a054 and tends toward unity for densities above 10 5 cm" 3 . 

To illustrate quantitatively the effect of density, we have calculated a sequence 
of constant pressure models covering the range 50< n f H < 2xl0 6 cm . The results 
of integrated nebulae calculations as function of the mean slab electron density are 
presented in Figure 5. The case of a pure recombination-dominated HII region 
model {dotted line) is also shown for comparison with the power-law model {solid 
line). In Figure 5a, it is evident that both ionizing energy distributions result in 
a qualitatively similar functional dependence of Lya/H/3 on < n e >, except that 
Lya/H/3 is consistently higher for the power law due to collisional excitation within 
the PIZ. For the Balmer decrement (Fig. 56), at low densities (< n e 5 x 10 4 
cm -3 ), however there is a sharp increase in Ha/H/3 owing to the enhancement 
of collisional excitation that occurs when the gas temperature increases due to 
gradual suppression of the forbidden lines. 

Subsequent models with internal dust have all been calculated with the single 
representative density of njj= 5000 cm~ 3 . This reference density coupled with 
the adopted Uf introduced below provides us with a comparative zero point (i.e., 
dust- free case) in our Lya/H/3 versus Ha/H/3 diagrams (as represented by the filled 
star). If one is interested in estimating the line ratios different zero point in density 
space, Figure 5 can be used to determine very approximately the differential change 
in the ratios due to density alone and apply it to the models which contain dust. 

2.3.2 Asymmetric emissivity and internal dust absorption 

Because line scattering increases the path length of Lya photons, the efficiency 
of internal dust to absorb line photons is drastically higher for the optically thick 
resonant Lya line than for the optically thin Balmer lines. Hummer & Kunasz 
(1980) have derived scaling laws for the energy loss, the mean path length, and 
the mean number of scattering suffered by resonant line photons in the presence 
of continuous absorption (by dust and/or photoelectric absorption). (See Neufeld 
1990 for a discussion of these scaling laws using analytical methods). The fraction 
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Figure 5: 

Sequences of dust-free solar metallicity photoionization calculations with varying gas 
density n H . Results for two different energy distributions of the ionizing continuum are 
represented: the power law with a = —1.4 (solid line) adopted in this work and a black 
body of 40,000 K (dotted line) for which recombination is the only process generating 
HI lines. The ionization parameter Uj is constant along both sequences and takes on 
the values of 0.0015 and 0.010 for the power law and for the blackbody sequences, 
respectively. Panel (a) shows the Lya/H/3 ratio ("f" + "b") as a function of the mean 
electron density < n e > characterizing the isobaric photoexcited zone, while panel (b) 
shows the corresponding behavior of Ha/H/3. The star represents the standard model 
(f7 /= 0.0015, n f H = 5000 cm" 3 ) discussed in §2.3. 
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of energy from a resonant line generated within a midplane layer which escapes 
unabsorbed from a slab of total line opacity T depends essentially on the parameter 
P = k c /k L which is the ratio of continuous opacity to line opacity. The results of 
Hummer & Kunasz (1980) concerning the transfer as well as the interplay between 
escape probability and dust absorption have been incorporated in MAPPINGS as 
described in detail in Appendix of BWZM. 

One aspect of the transfer particularly relevant to thick nebulae is the asym- 
metry which characterized the Lja flux emerging from the front or from the back 
of the ionized slab. This is taken into account in the code through the definition of 
an emissivity whose asymmetry depends on the local two-sided escape probability. 
The asymmetry factor A^(0 < A? < 2) in our transfer solution is defined as follows 

A' = 2e(T,P)/[e(T,j3) + e(r-T,j¥)] [4] 

where r is line opacity of the layer under consideration, T(= / dr) is the 
line opacity of the whole slab, e(r, f3) is the escape probability from the back. 
(See Appendix A in BWMZ for further details on the escape asymmetry). The 
implication of the two-escape scheme is that the emissivity of the resonant line is 
asymmetric with the forward emissivity given by j h L = j'l(2 — A*'). The observer, 
in this simple two-stream transfer scheme, is placed either in front of the slab or 
behind it and the emergent line flux from either perspective "f" or "b" is obtained 
by integrating the appropriate emissivity (j[ or j b L ). 

An important consequence of using equation (4) in integrating the line emis- 
sivity is that for a slab illuminated on only one side by the ionizing radiation, the 
Lja intensity emerging from perspective "f" is much higher than that from "b" 
whenever the line is reasonably thick (r > 10), a condition almost always satisfied 
except for extremely low column density matter-bounded slabs. For a slab of total 
line opacity T, the depth Nh at which line photons have equal probability of escap- 
ing from the front ("f") as from the back ("b") occurs at the depth Nh{t = T/2), 
i.e., N H = N h {t = T/2). Since dr/dN H oc n H o/nn is a monotonically increasing 
function of Nh, the depth Nnij = T/2) is quite generally located much deeper 
than the depth that divides equal Lja production. As a result, Lja photons pref- 
erentially escape from the front ( "f " ) . This is true even for soft impinging ionizing 
continua as occurs, for instance, in HII regions. For a power-law ionizing contin- 
uum, to the extent that recombination still dominates Lja production, most of 
Lja is produced within the fully ionized zone, i.e., at Nh < 0.8N H+ which is much 
less than N H (r = T/2) ~ N H H* owing to the extended PIZ (cf. § 2.2.3). In the 
opposite case of the extremely low column density nebulae, the neutral fraction 
n H°/ n H ~ constant so that Nh{t = T/2) pa N^ ab /2, which coincides approxi- 
mately with the depth dividing equal Lja production. In this case, the Lja flux 
emerges symmetrically from the front ("f") and the back ("b"). 

To illustrate how one progressively goes from symmetrically emerging Lja 
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in the very thin slab case to the opposite asymmetric situation in the radiation- 
bounded slab case, we have computed a sequence of matter-bounded models in 
which N^ ab is gradually increased in lock steps. Figure 6 gives the fraction of 
Lya which escapes from perspective "b" as a function of the matter bounded 
slab's thickness (expressed in a different way in each panel). Note that internal 
dust has been included in these models, but it does not affect significantly the 
degree of asymmetry except when N^ ab approaches Njj*, the radiation-bounded 
case. Two of the sequences correspond to Uf— 0.015 with fx—1.0. The results 
indicate that the fraction of Lya which escapes from the back is close to half only 
when the matter-bounded slab is indeed very thin (i.e., for r^- (912A ) < 10 or 
N siab/ N s + < 0.5 -0.8). 

The effect of perspective need not be considered in previous work on the 
intrinsic NLR Lya/H/3 and Ha/H/3 line ratios since it was implicitly assumed that 
perspective effects averaged out due to the geometry of the ensemble of NLR clouds. 
This might be a valid assumption in the dust-free case where no Lya destruction 
takes place despite the high number of scatterings. When dust is present, however, 
the effect of perspective cannot be ignored since the Lya transfer with internal dust 
present is nonconservative with respect to perspective with most of the surviving 
Lya photons escaping from the front. 

Perspective can also affect significantly the integrated Balmer lines of a sys- 
tem of clouds with internal dust because the dust within the extensive PIZ signi- 
ficantly absorbs the Balmer lines which emerges from the back (producing at the 
same time the steeper Balmer decrement). This effect is particularly important in 
understanding the results presented in § 2.3.4. To illustrate the dominant effects of 
both perspective and dust, we picture a thin fully ionized region which generates 
the bulk of Lya, and, depending on how opaque the region is to dust [t v (H + )\, a 
substantial fraction of the resonant Lya escapes from the front while half of the 
Balmer lines' flux is required to transfer across the larger and dustier PIZ before 
escaping from the back (the other half emerges from the front after suffering little 
extinction). 

To illustrate the role of internal dust on the total emerging Lya and Balmer 
lines, we present in Figure 7 the line ratios resulting from adding together the fluxes 
from perspective "f " and "b" (before taking the ratio) as a function of the Lya limit 
opacity r^ (912A ) for the matter bounded sequences presented in Figure 6. The 
larger destruction rate of Lya for the high- ionization parameter {dotted line) model 
is apparent and is the result of a higher value of f3 throughout the nebula (cf. § 
2.3.3). We also find for all three sequences that for t^ (912A ) > 100, because dust 
in the PIZ starts absorbing much of the R/3 flux escaping from the back, the "f+b" 
Lya/H/3 ratio rebounds up (the Lya flux actually remains steady). The somewhat 
higher Balmer decrement in the matter-bounded regime in the higher U f model is 
generally caused by reddening [larger Tv{H + )]. When r^ (912A ) < 10, then in 
addition to reddening, for (3 sufficiently high, the normal conversion of the higher 
Lyman lines into Balmer or higher series is reduced because of dust absorption of 
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Figure 6: 

Matter-bounded calculations in which the ordinate represents the fraction of Lya flux 
which escapes from the back (perspective "b") as a function of different quantities which 
measures the ionized slab's thickness. The abscissae represent, panel by panel, the 
following: panel (a): the fraction of unabsorbed ionizing photons which leaks out from 
the back panel (6): the slab's column density to that of the Stromgren depth 

Ntf+ = clIf/aBj where as is evaluated at the slab mean electronic temperature; panel 
(c): the hydrogen ionized fraction within the last layer at the back; panel (d): the slab 
Lyman-limit continuum opacity. The solid and short dash curves were both calculated 
using the reference value of Uf = 0.0015 (and n^ H = 5000 cm -3 ; see §§ 2.3.2 and 2.3.3) 
but with an internal dust content of fi = 1.0 and 0.2 respectively. The dotted line 
corresponds to calculations with L r /=0.015 and /x=1.0 
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Figure 7: 

In panel (a) , Lya /H/3 as a function of the Lyman- limit opacity of the matter-bounded 
photoionized slabs. As in Fig. 6, the solid and short dashed curves were calculated with 
Uf=0.0015 and an internal dust content of /i=1.0 and 0.2 respectively. The dotted line 
corresponds to calculations with ?7/=0.0015 and fx = 1.0. Panel (b) plots the Balmer 
decrement as a function of the Lyman-limit opacity. 
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these Lyman line (this effect simulates something resembling case A). 

2.3.3 Selection of the ionization parameter Uf 

The most simple geometrical model of the narrow-line region consists of a symmet- 
ric distribution of ionized clouds surrounding the ionizing source. Any reasonable 
assumed electron densities leads us to infer a very small volume filling factor for 
the emitting gas. The ionization parameter which is inferred from the line ratios 
is surprisingly uniform within a given class of AGN. 

To determine Uf, we have simply required the radiation-bounded calculations 
to approximately fit the "mean" optical Seyfert 2 line spectrum derived by Ferland 
& Osterbrock (1986, hereafter F086) in which a ratio of [OIII]A5007/H/3 of ~ 10 
is observed. With the other parameters specified above, this leads us to adopt 
Uf =0.0015 for the reference ionization parameter. 

The choice of the ionization parameter affects strongly the role played by 
dust. For instance, internal dust can absorb a fraction of the ionizing radiation 
and therefore shrink the Strongrem depth due to competitive absorption with H°. 
In the case of our reference model with Uf =0.0015, only 3% of 4> H is lost to dust 
absorption when /i = 1. However, with an ionization parameter 10 times larger 
[recall that Ty(H + ) scales with N^+, we find that as much as 24 % of <p^ H is 
absorbed by dust when u,=l. This effect can be seen in the Uf =0.015 sequence 
shown in Figure 6b where the dotted line curve is shifted to the left because scaling 
is relative to the dust-free value of Nfj+ . 

Another aspect where Uf is important is the amount of destruction of resonant 
Lja by dust. To show this, let us define using equation (2) the pure absorption 
opacity by dust at depth half way within the fully ionized zone: 

i(1216A l -N s H+ ) - 2.2r-< = mtxUf [5] 

The justification for choosing this position is that most of the Lya is generated 
as well as destroyed within the fully ionized region. Assuming Uf = 0.0015 and 
pi = 1.0, we obtain T^ s st (121QA ) 0.1. As a crude estimate to illustrate the effect 
of dust on Lya, at this depth, the line (scattering) opacity is ~ 10 5 which allows 
us to estimate (3 ~ 10~ 6 and derive the fraction of Lya which escapes unabsorbed 
from this layer at ~ 50% by using the scaled results of Hummer & Kunasz (1980). 
On the other hand, for the higher value of Uf=0.015, the line opacity at a similar 
relative position corresponding to r^ s st (121QA ) ps 1 (/i=l) is higher (~3x 10 6 ), 
and we estimate that with (3 ~ 3 x 10~ 7 , only 2% of Lya escapes the layer. For 
calculations with the adopted ionization parameter Uf =0.0015, we see that Lya 
destruction due to dust should be moderate. 
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2.3.4 Results for a symmetric cloud distribution 

We approximate the typical NLR emission cloud as a slab of gas of constant internal 
pressure with an outer density of n H =5000 cm~ 3 (cf.§ 2.3.1). Using MAPPINGS, 
we compute the line spectrum from photoionization calculations using the ioniza- 
tion parameter L r /=0.0015 (cf.§ 2.3.3) and an ionizing power-law of index a = —1.4. 
In all the calculations, the gas phase abundances are modified consistently with 
the dust content /x in such a way that the total metallicity (gas+dust) remains 
always solar: Z = Z gas + Z dust =1 (cf.§ 2.2.2 and Appendix A). The assumed 
symmetric distribution of the clouds is taken into account in the calculated spec- 
trum by adding together spectral "f" and "b" of MAPPINGS calculations which 
is equivalent to supposing that there are as many clouds seen face-on as are seen 
from the back. 

Using this procedure, we have computed three sequences of models with the 
dust content \x covering the range < /i < 1.25. These three sequences are 
presented in Figure 8 (solid lines) and correspond to three different total slab 
column densities N^ ab . The intermediate sequence corresponds to the radiation- 
bounded case Nfjj ab = Nh* (see § 2.3.3 for the operative definition of the depth of 
the photoexcited region Nh*), a second sequence represents the matter-bounded 
case in which Nfj ab corresponds to the depth where hydrogen is 10% ionized (i.e., 
jysjab ^ 2.4 x 10 20 cm" 2 ), a third one is the case in which the NLR clouds are 
more massive than the radiation-bounded case and has Nfj ab = Nh* + N H0 = 5 x 
10 21 cm -2 , which is about 6 times the Nh* of the dust-free case 0. The zero-point 
dust-free reference model with L r /=0.0015 is represented in all figures by the filled 
star. We emphasize that our simple geometry assumes a very low volume filling 
factor of the emitting clouds so that as many "f" clouds as "b" clouds are seen by 
the observer without any intervening obscuration. 

The most important result of the calculations shown in Figure 8 is to prove 
that despite large amounts of internal dust, the Lja/H/3 calculated using canon- 
ical input parameters is not radically different from recombination case B, and, 
if anything, it can even be larger when /i — > 1. The expected range for Lya/H/5 
in the dust-free recombination case B is between 23 and 34 depending on the gas 
density. If, in addition, collisional excitation is present due to the hardness of the 
continuum or to low gas metallicity, both the Balmer and the Lyman decrement 
become larger (see Gaskell & Ferland 1984), a situation not unlike that shown in 
Figure 8 except that here it is the result of internal dust. 

The two main factors explaining the relatively high Lya/H/? are that, first, 
with Uf =0.0015, Lya destruction is modest in the clouds since the column density 

2 Note that iVy* progressively increases with fi, from 8 x 10 20 cm~ 2 when /z=0 to 2.6 x 10 21 
cm~ 2 when /j=1.25. The main reason has to do with the criterion requiring that hydrogen be 
only 1% ionized at the inner boundary of the photoexcited region. Because the recombination 
coefficient depends on T e , Nh* becomes larger when the temperature of the PIZ increases as a 
result of depletion of the gas phase with increasing \i. We have found that T e increases from 1200 
K in the dust-free case to 3300 K at /i=1.25 
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Figure 8: 

Lya/H/3 as a function of the Balmer decrement for the three sequences with Uj =0.0015 
in which /j, is monotonically increasing (0< fi < 1.25). Tick marks from left to right 
along the solid lines indicate values of fi of 0.1, 0.2, 0.4, 0.6, and 1.0, respectively. 
The three sequences (solid lines) differ only by the total column density of the slab 
N'lj ab . The lowermost sequence is matter-bounded using the criterion that the hydrogen 
ionized fraction is 10% at the back layer. The intermediate sequence corresponds to 
the radiation-bounded case without any additional neutral gas zone at the back. The 
uppermost sequence includes a neutral back zone such that the total column density of 
the slab is constant at 5 x 10 21 cm -2 . The neutral and photoexcited zones are assumed 
to have the same dust content [[i(H°) = /j,(H*)]. The filled star represents the reference 
dust-free case (/ti=0). Also shown in the figure is an ionization parameter sequence 
(dotted line) for radiation-bounded slabs (with N^ = 0) of constant dust content /i = 
1.0 and covering the range 10~ 4 < Uf < 10~ 2 . In all above calculations, the ionizing 
energy distribution is a power law with index a =-1.4, the gas density is n^ H = 5000 
cm' 3 at the face of the cloud, the metals are depleted from the gas phase in accordance 
with the dust content \i, and the line intensities are given by the sum of the "f" and "b" 
perspective. 
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of ionized gas and therefore Ty(H + ) is not large, and second, the portion of the 
Balmer lines' flux which escape from the back is subject to considerable absorption 
by the dust inside the PIZ or within the neutral core of the cloud. It is therefore not 
surprising that the models which depart most from case B Lya/H/9 are those of the 
matter-bounded sequence, comparatively little absorption occurs for the Balmer 
lines while destruction of Lya photons (which occurs mostly in the fully ionized 
zone) remains significant (cf.§ 2.3.2). To illustrate to what extent much higher 
values of Uf would be catastrophic to the escape of Lya photons, the dotted line 
in Figure 8 represents an ionization parameter sequence with models of constant 
dust content (/x=0.1). Of course, these models are of limited relevance since they 
result in a forbidden-line spectrum incompatible with observed optical line ratios. 

Since in the dust-free case we used solar metallicity for the gas, little collisional 
excitation takes place within the PIZ and the resultant Balmer decrement is not 
significantly above that of a dust content of fx—1, however, the Balmer decrement 
in the radiation-bounded case is shown in Figure 8 to approach 3.3. On the other 
hand, moderate amounts of internal dust /i ~ 0.2 — 0.4 would favor an intrinsic 
Balmer decrement in the range 3.0-3.1 assuming the clouds are radiation-bounded 
(the corresponding range the Lyman decrement is 32-37). This range (3.0-3.1) 
corresponds to that proposed by various authors (Gaskell & Ferland 1984; Halpern 
& Steiner 1983) for the intrinsic NLR Balmer decrement except that it is here more 
the result of the presence of internal dust rather than of a hard ionizing continuum 
or of deficiencies in the gas metallicity relative to solar. Binette et al. (1990) 
pointed out the apparent lack of Seyfert 2 with Hct/H/3 < 3.4 if one considers only 
objects with good S/N in the H7 line. If this high value of 3.4 was exclusively 
caused by internal dust present in all Seyfert 2's, then large values of pL (~1) 
are implied from Figure 8. Alternatively, a portion of the extinction may be the 
result of intervening gas containing dust, in which case the amount of internal 
dust implied is lower. There are two possibilities for how this intervening gas/dust 
may arise. In one case this dust is contained in a detached screen. Such models, 
however, may not be physically plausible as discussed in §4.2. Another possibility 
is that the intervening dust+gas arises from a large covering factor of the NLR 
clouds themselves. This case is presented in §2.4.3. 

Radiation-bounded calculations with our selected parameters are known to 
over predict by a factor ^2 the [OI]A6300/H/3 ratio. Using a moderately matter- 
bounded cloud such as the one shown in Figure 8, one can resolve this discrepancy 
by cutting the cloud at some depth to make the [OI] fit the observed value. As 
previously shown in Figure 7, the matter-bounded model /i—l results in a lower 
Lya/H/9 (than radiation-bounded calculations) somewhat below case B and in a 
moderately reddened Balmer decrement. The procedure of cutting the clouds to 
adjust [OI] to the observed value is, however, quite arbitrary. A more realistic 
approach is the dual gas components model of Viegas & Prieto (1992) which was 
proposed as a solution to explain the observed high HeII/H/3 ratio observed in the 
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spatially extended NLR (ENLR) (see also Morganti et al. 1991). This model con- 
sists of a matter-bounded (MB) and a radiation-bounded (RB) component. The 
MB component has a higher volume filling factor (to compensate for its much lower 
emission measure) and a higher ionization parameter (Uf >0.0015). Although we 
did not consider this level of complication in our modeling, we can still establish 
some of its general properties in relation to dust using the results of §§ 2.3.2 and 
2.3.3. Under most instances, the effects of dust on the HI lines cannot be larger in 
the case of the dual- component model than with our simpler single RB component 
because the MB component, despite its higher Uf, is nevertheless extremely thin 
in column density and therefore generally results in little Lya destruction []. Fur- 
thermore, the RB component should have U RB < 0.0015) (Viegas & Prieto suggest 
jjRB < 0.0006) and therefore much lower ry than in our single RB component. In 
conclusion, a working dual-component model would in all likelihood be character- 
ized by smaller line absorption effects due to dust at equal values of \x than in our 
simplified single component model. 

2.4 Discussion 

After comparing calculations with internal dust with recent Lya/H/? data on 
Seyfert 2's (§2.4.1), we discuss various interpretations of the low values observed. 
We look at the possibility of pure line scattering or pure dust absorption (§2.4.2) 
as well as that of mutual covering of Lya clouds (§2.4.3). On §2.4.4, we consider 
the possibility of dust absorption of the UV continuum for explaining the apparent 
deficit of ionizing photons. 

2.4.1 Observational data of KAW3 

The available data on the NLR Lyman decrement is still rather limited. The 
most reliable information on the NLR comes from Seyfert 2's since, as pointed 
out by F086, only for these objects can one deduce an optical to ultraviolet line 
spectrum which is totally free of contamination by BLR components. It is likely, 
however, that the physical information gathered from the NLR of Seyfert 2 is 
relevant to that of Seyfert 1 and quasars. One of the motivation of the current 
work was the recent publication by Kinney et al. (1991) of a Seyfert 2 data set 
for which the optical spectra were taken through an aperture which matched the 
large IUE satellite aperture. The matching of aperture size for the different lines 
is a crucial point in the study of the NLR Lyman decrement because the emitting 
line region is spatially extended and in the likely event that Lya, Ha, H/5 present 

3 The dual component model of Viegas & Prieto (1992) suggest a MB component with opacity 
t^ (912 A ) < 1. Adopting t^ (912 A ) =0.55 and a similar U MB = 0.01 to what they use, 
we obtain Nfj ab = 0.067 fi. Even with high dust content fi=l, our calculations indicate a Lya 
destruction of only 0.13 dex with a difference in flux between "f" and "b" perspectives of 15%. 
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a different spatial distribution in surface brightness, the interpretation of ratios of 
lines derived from different apertures becomes very complicated, if not impossible. 
By comparing for instance objects common to both the KAW3 and F086 data 
set, one finds many objects with significantly different line ratios between the two 
data sets presumably as a result of the much smaller aperture of the F086 optical 
observations. For this reason we have limited our analysis to the data set of KAW3 
which are presumed free of this problem. 

In Figure 9 (which is an expanded scale version of Fig. 8), we present the 
data of KAW3 (filled squares) with their estimated error bars. The earlier model 
sequences of Figure 8 now appear as broken dash lines in the upper left. What is 
immediately striking is that internal dust absorption for an isotropic distribution 
of emitting clouds does not account for the low Lya/H/5 ratio which is significantly 
below case B. 

In order to explain the overall low Lya/H/3 ratio in Seyfert 2's, either an 
intervening absorption/scattering screen is required or we must modify our simple 
NLR geometry. The first possibility, discussed in detail in §4.2, would require pure 
line scattering from the screen to explain the near-case B Ha/H/3 objects. The 
second possiblity is that the covering factor of the emission clouds is not small as 
was assumed so far. This idea is explored in §2.4.3 along with the possiblity that 
if the covering is extreme, a closed spherical geometry might be preferable to that 
of an isotropic distribution of slab-approximated clouds. 



2.4.2 The pure line scattering and the pure dust models 

There are broadly two types of Seyfert 2 objects in the KAW3 data: the objects 
with near-case B Hcc/H/3 and the objects with reddened B.a/B.(3. The first cat- 
egory might be explained by pure Lja scattering due to an intervening screen 
of H° which is supossed to be covering only the observer's view of the NLR but 
not the complete sky as seen from the NLR. Assume a slab geometry for the 
screen. Then for a pure scattering screen of thickness r H o scat , where T H a scat = 
5.9 x 10~ 14 Nffj een (T HI /10 4 K)~ 1 / 2 is the Lyct resonant scattering optical depth at 
line center, the fraction of photons impinging on the screen from the side facing 
the source which penetrates to the opposite side (the side facing the observer) is 
oc T H0 ' scat for T H o scat 3> 1 (Slater, Salpeter, & Wasserman 1982). Owing to the 
large magnitude of r H o scat , only modest values of Nffj een are required to efficiently 
suppress the Lja flux reaching the observer. Thus, for instance, for Nffj een ~ 10 15 
cm" 2 (10 17 cm -2 ), only 10% (1%) of the incident Lja photons reaches the ob- 
server; the rest are reflected back toward the source. At such column depths, the 
dust opacity within the screen is entirely negligible, and so the extinction vector + 
shown in Figure 9 is purely vertical. 

The lack of objects with near-case B Ha/H/5 in the F086 sample has led 
these authors to propose pure dust absorption (ignoring the destruction of Lja 
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Figure 9: 

The squares represent the homogeneous data set of Seyfert 2 observations of Lya/H/3 
vs. Ha/H/3 by Kinney et al. (1991) in which the aperture in the optical closely matches 
the large-aperture IUE Lya observations. Two reddening vectors corresponding to the 
ISM and the Orion extinction curve, respectively, are drawn to indicate the expected dis- 
placement due to foreground pure dust. The vectors' length correspond to an extinction 
r = 1 in the visible (5500 A ). The vertical solid line represents the effect of scattering 
of Lya by a foreground neutral hydrogen screen at the same redshift as the NLR (see 
§ 4.2). The dash-dotted lines represents the \x sequence with Uf =10~ 2 ' 8 and = 
but as seen separately) from perspective "f" and "b". The three dash lines correspond 
to the same ji sequences of Fig. 8 (i.e., with "f" + "b") 
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after multiple resonant scatters) since their objects all broadly lie along the ISM 
reddening vector (see Fig.6 F086). Within the sample of KAW3, many objects 
still lie in the general direction of the ISM reddening vector. The pure dust screen 
(without any line scattering; otherwise, the reddening vector would be steeper 
than shown above) remains a possible interpretation for the observed ratios, but 
only for the reddened Ha/H/3 objects. Several limitations exist, however with the 
picture of the pure dust screen which we now proceed to analyze. 

In the pure dust screen model, a layer of dust and gas is presumed to lie 
outside the NLR toward the line of sight to the observer. Photons that pass 
through this layer are assumed to suffer solely absorption by the dust, that is, 
interaction between the photons and the gas are ignored. While this is a good 
approximation for non resonant line photons such as Ha and H/3, it is a poor 
approximation for resonant line photons such as Lya. Owing to the large resonant 
scattering cross section of Lya photons, resonant scattering by these photons with 
neutral hydrogen atoms as they traverse the layer cannot be ignored in general. 
For example, to obtain a reddening of Ay = 1 require Nh = 2 x 10 21 cm~ 2 for 
/j, — 1. For a gas at T e =10 4 K, this corresponds roughly to Nhi ~ 10 20 cm~ 2 , 
which implies a line center resonant scattering optical depth fort Lya photons of 
reseat ~ 6 x 10 6 . 

One possible way to avoid resonant scattering so that only pure dust absorp- 
tion is operative is to have the layer move with sufficiently large bulk velocity 
so that the line is shifted out of the Doppler core. How large a column density 
in neutral hydrogen can be tolerated without significant resonant scattering then 
depends on how fast the layer is moving away from the line source. We plot in 
figure 7 this critical neutral hydrogen column density, N^} Uc , against the layer 
bulk velocity for gas at T e = 10 2 , 10 4 , and 10 6 K where N^f lc is defined as the 
Nhi (measured in the comoving frame) that gives T H o scat =l (see Appendix B, § 
1). Thus, resonant scattering may be ignored for layers with a neutral hydrogen 
column density lying below the curves. The "steps" in the curves correspond to 
where the line is shifted from the core into the wings. It is evident that N^f tc 
is very low < 10 14 cm -2 ) unless the bulk velocity is very high. In order to have 
Ay > 1, we require N H i > 10 20 cm~ 2 ) (assuming T e = 10 4 K and // w 1), which 
requires V > 500 km s _1 ) from figure 7. Aside from the issue of how to accelerate 
the layers to such high bulk velocites, there remains the problem that these clouds 
must be constantly replenished from some unknown source. This replenishing must 
be very efficient and continuous because the narrow-line gas in the dust screen is 
at T e =10 4 K) and has Nhi > 10 20 cm , the screen itself should have significant 
emissivity in lines such as Ha and H/3, unless the screen is placed sufficiently far 
outside the NLR so that geometric dilution reduces the impinging ionizing flux to 
insignificant levels. Otherwise, one should see highly blueshifted lines (from the 
screen) in addition to the usual lines from the NLR. This is not observed. 

Without invoking bulk motion, another way to make the pure dust screen 
viable is to make the gas in the dust screen very cold relative to the width of 
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Figure 10: 

The critical neutral hydrogen column density, Nfjf , of the dust screen vs. the screen's 
bulk velocity, V. The quantity Nfj'f is defined as the column density (measured in the 
comoving frame) giving Lya resonant scattering optical depth of unity (T H o scat = 1). 
The region below the curves have r H o scat < 1. The curves correspond to a dust screen 
gas temperature of 10 2 K (dashed), 10 4 K (solid), and 10 6 K (dot-dashed). The "step" in 
each of these curves correspond to where the line is shifted from the core into the wings. 
The impinging Lya line is assumed to be monochromatic at the line center (as measured 
in the line source frame). 
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Figure 11: 

The fraction of energy, /, in the Lya line impinging on the dust screen that suffers 
interactions with the screen's gas vs. the neutral hydrogen column density, Nhi, of 
the screen. The horizontal dotted line correpond to /=0.01. The curves correspond to 
k = SVs/SVa = 10 (solid), 10 2 (dashed), and 10 3 (dot-dashed), where 5Vs is the width 
(in velocity units) of the impinging Lya line and 5V a is the width of the Lya resonant 
scattering line profile for the gas in the dust screen. 

the line impinging on the screen. In this manner, only a vey narrow band about 
line center will suffer resonant scattering within the screen. In Figure 11 we plot 
the fraction of the impinging energy in the Lya line (assumed to have a gaussian 
profile) that suffers interactions with neutral gas atoms in the screen as a function 
of the neutral hydrogen column density of the screen (see Appendix B, § 2). The 
horizontal line shows where 1% of the incoming line energy suffers interactions. 
The curves are labeled by the width of the impinging line (expressed in velocity 
untis) relative to the width (in velocity units) across which the line will interact 
with the gas in the screen. Typical line widths for the NLR Lya emision line is ~ 
10 3 km s _1 . It is evident from Figure 11 that even for (screen) gas at T e =10 2 K so 
that k = 5v s /Sv a = 10 3 , the critical neutral hydrogen column density (to have 1% 
interaction or less) is only ~ 10 16 cm~ 2 ~ Nh, corresponding to Ay ~ 5 x 10 -6 //. 
Thus, a cold screen is unfeasible because the low hydrogen column densities it 
requires (to avoid resonant scattering) imply negligible absorption by dust internal 
to the screen. 

We conclude that either with bulk motion or with a very cold cloud, the pure 
dust screen model is not physically plausible. The only scenario we can conceive 
of where a pure dust screen might work is a "coronal" model wherein the dust is 
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embedded in a very highly ionized, and hence very hot (T e ~ 10 6 K), and tenuous 
(n e ~ 10 cm" 3 ) medium. Assume for now that dust can exist in significant quatities 
in this coronal gas. To have A v ~ 1 require N H ~ 10 21 cur 2 for /i=l, or N H ~ 
10 22 cm" 2 for fj,— 0.1. With n e ~ 10 cm -3 , this correpsonds to a coronal region 
with linear size L ~ 10 2 pc (/i=l) or L ~ 10 3 pc (fi = 0.1). In the former case, 
the size is comparable to the size of the NLR. While the existence of such a region 
external to the NLR cannot be ruled out, it is somewhat ad hoc. In the latter 
case, the large volume required for the coronal region almost surely implies that 
the corona would cover the whole NLR. In this case, one would expect the same 
dust in the coronal gas that affects the narrow lines to redden the starlight from 
the central galactic buldge. This reddening is generally not observed. 

The properties of the coronal gas is probably not all that different from those 
of the intercloud medium in the NLR, if we assume a two phase medium for the 
NLR with the two phases (cloud and interclouds) being in pressure equilibrium. If 
so, then postulating a separate corona outside the NLR seems superflouos; it would 
be more natural to assume dust to be present in the NLR intercloud medium itself. 
If, however, one is to consider dust in the intercloud medium, then one should also 
consider dust in the clouds themselves. The narrow lines produced within the NLR 
cloud would then be affected in general by dust in both the cloud and intercloud 
medium. With regard to dust survival, it is likely that both evaporation (owing 
to the high gas temperature) and sputtering will render the dust life-times to be 
shorter in the coronal-like intercloud medium than in the clouds. The main dust 
destruction process inside the clouds is sputtering which gives a dust lifetime of ~ 
10 5 yr (assuming n e ~ 10 3 cm" 3 and T e ~ 10 4 K). In the NLR model we adopt, we 
assume that dust is able to be replenished on such time scales inside the clouds to 
enable a significant dynamic steady state abundance of dust to be present. This 
replenishment is probably more difficult to achieve in the coronal-like intercloud 
medium and in our model, we assume that dust is not present in any significant 
quantities in the intercloud medium. 

2.4.3 The partially covering models 

One obvious advantage of a simple interpretation for the position of all the objects 
in Figure 9 rather than the two separate pure dust and pure line scattering inter- 
pretations is that it might be more informative of the NLR geometry. Although 
the Lya/H/5 observed is much below that predicted by the "f" + "b" models, some 
of the objects are not that far from the pure perspective "b" models (see dot- 
dash line in Fig. 9). This naively suggests a nonnegligible covering factor for the 
emission region in which individual clouds partly cover one another in the line of 
sight to the external observer. Although plane-parallel calculations do not allow 
us to replicate very accurately such a geometry, a first order estimate of the effects 
involved can nevertheless be arrived at. For instance, we can safely expect the 
Lya luminosity to be dominated by the few uncovered clouds which are directly 
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seen from the front (i.e., "f" clouds). The Balmer lines in our simplified scheme 
would essentially come from three different components: the few "f" clouds seen 
directly, the many "b" clouds seen from the back, and finally the "f" clouds seen 
through other clouds and which we postulate are seen through a depth statistically 
equivalent to a single cloud (we neglect multiple cloud covering). By assuming that 
the transmitted spectrum of each covered "f" cloud is approximately equivalent 
to that of a "b" cloud Q, we can readily estimate the main effects of a partially 
covered geometry by simply giving more weight to the "b" clouds, i.e., "f"+7"b". 

In Figure 12, we present the loci of "f" + "b" calculations in which the weight 7 
takes on the values of 1, 3, 10, 30, 100, and 300 using photoionization calculations 
in which the dust content u.(H*)=0.2. For a given value of 7, the solid line tracks 
result from varying the depth of neutral gas N^ (beyond the PIZ). Since the dust 
content within the neutral zone n{H°) (see eq. [3]) might be plausibly be different 
from /u,(H*), we specify N^o in Figure 12 in terms of the dust opacity in the 
visible, Ty(H°) [leaving (jl(H°) unspeciefied; see tick marks]. Considering the size 
of error bars (dotted lines) which characterizes the data, we find very encouraging 
that the tracks pretty much cover the whole region occupied by all the Seyfert 2 
observations. 

Even if the tracks are only crude estimates and are not unique in terms of the 
parameter n(H*), they succeed in illustrating that a very simple and self-consistent 
mixture of internal line scattering and dust absorption can explain the observed 
ratios without the need for ad hoc intervening screens. We infer from this exercise 
that the NLR HI line ratios points toward a significantly self-covered distribution 
of gas which is very patchy and where all the Lya originates from ionized gas in 
direct view of the observer. 

If this conclusion is valid, it is possible to study one limiting case of this picture 
by supposing the covering factor to be so high that spherical geometry becomes 
a preferable representation of the NLR (and which at least addresses the effects 
of multiple scattering inside the geometry, a possibility excluded when considering 
only independent slabs). 

Adopting spherical geometry, we will represent the NLR as a Stromgrem 
sphere uniformily surrounded by a layer of neutral gas of opacity Ty(if°) except 
for a a small area where the ionized sphere is not covered as depicted in Figure 13. 
As in the simple partial covering model dicussed above, Lya escapes exclusively 
from the region which is uncovered. The Balmer lines, on the other hand, emerge 
through both the uncoverd area and the neutral dust outer layer. 

To compute in spherical geometry and yet retain an adequate treatment of the 
transfer through dust, we integrate the emissivity of volume cells within either the 
near-side hemisphere "b" or the far-side hemisphere "f". For the line transfer we 

4 A zero-order aproximation based on the fact that the fully ionized zone is only a small 
fraction of the whole cloud and therefore any "b" spectrum is really seen through a screen of 
dust of depth comparable to that given by the whole slab. 
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Figure 12: 

Same observations as in Fig. 9 but compared with the partially covering model using 
n{H*) = 0.2 (see eq.[3]) and with the line fluxes obtained by giving a higher weight 7 to 
perspective "b" ("f" + 7"b"). For a given value of 7, the solid line illustrates the effect 
of increasing the neutral zone depth, N^ , characterized by a dust content n{H°) =1.0. 
The tick marks from left to right correspond to opacities of the A^ zone of ry = 0.05, 
0.19, 0.5, 1.2 and 1.9 [for example, with fi(H°) = 1.0 these correspond to N^ = 10 20 , 
4 x 10 20 , 2.5 x 10 21 and 4 x 10 21 cm -2 , respectively]. The three dash lines correspond 
to the same /x sequences drawn in Fig. 8 ( "f " + "b" ) 
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Figure 13: 

Schematic depiction of the spherical geometry assumed for the NLR. The photoexcited 
sphere of radius R* and dust content fJ>(H*) is divided into two hemispheres "f" and 
"b" . The integral of the emissivity of each hemisphere is integrated along 10 concentric 
cylinders whose axis defines the observer's perspective. The photoexcited sphere is sur- 
rounded by a neutral shell of gas of dust content ^(H°) and opacity t v (N^ ). A small 
area of the sphere is uncovered from which some Lya escape probability is determined 
radially starting from the inner gas radius i?o- 

use the dust opacity occurring between each cell and the surface of the sphere along 
the direction of the observer's line of sight. The ionization and thermal structure 
are calculated, as usual, radially from the inner gas radius Ro up to the limit of the 
photoexcited region R*. For the line opacity, the integral extends out to include 
the outer neutral layer. The line emissivity cells are denned by the intersection of 
the photoionized spherical layers (used in the calculation of the radial ionization 
structure) of either hemispnere : "f " or "b" with concentric cylinders oriented along 
the lines of sight to the observer. In this scheme, we are in effect simulating the use 
of 10 concentric ring apertures of an observer which is looking at the spherical NLR. 
The global model is the sum of 10 "f"+ "b" cylindrically integrated contributions. 

Supposing no external neutral gas layer, we find that for the same amount 
of dust f^(H*) and column density of ionized gas N^+, the spherical geometry 
results in a reduced Lya/H/5 compared to the slab "f" + "b" case, resembling in 
fact more the "b" perspective of the slab case. This is understandable given that 
the direction of the gradient ionization structure (n#o/n# increasing with radius) 
seen by an observer external to the H + sphere is similar to that of a slab looked up 
from perspective "b". We emphasize that the effect of dumpiness of the ionized 
gas (which we do not consider in the transfer) would be more consistent with the 
geometry of a large number of clouds superposed along the line of sight. This 
would allow a larger fraction of Lyct to escape the ionized region as argued by 
Neufeld (1991), although the implementation of such effects is beyond the scope 
of this work. 

To derive the emergent spectrum for the partially uncovered case, we simply 
combine the spectra of the fully covered sphere -excluding the contribution of a 
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certain number of aperture rings- with the correponding missing aperture rings of 
a spherical model calculated without any external gas layer covering the "b" hemi- 
sphere. Note that in the dust-free case, such procedure would underestimate the 
emerging Lja by a factor of 2. With the presence of internal dust, the treatment 
is only approximate and rests on the reasonable assumption that most of the Lya 
scattering is local to the point of emission an approximation which is increasingly 
valid as \i increases since a photon which would scatter too many times over a 
length scale of order will likely be destroyed. Furthermore, the resonant pho- 
ton does not know a priori where the hole in the outer layer is and will likely be 
absorbed in the much larger neutral outer zone which covers the ionized sphere. 

The results are presented in Figure 14 where the amount of uncovered area is 
expressed as a fraction of the total projected emitting area 7ri? 2 . The parameters of 
the models were derived as follows. As for the slab case, we adjusted the ionization 
parameter of the spherical radiation-bounded models until [OIII]A5007/H/3 ~ 10 
(by varying the source luminosity in the center and/or the volume filing factor of 
the infinitesimal gas filaments). The ionizing continuum is the same as before, the 
density at Rq is 5000 cm~ 3 and the gas pressure is maintained constant with radius. 
The resulting photoexcited sphere in our calculations has R*/Ro ~ 2.4 and N H * ~ 
6.3 x 10 20 cm~ 2 . The electron column density in N e = 2.1 x 10 20 cm -2 which is 
exactly the same as for the reference slab model U =0.0013. Because more Lya is 
destroyed in spherical geometry, we have assumed smaller values for (/i(H*), [See, 
for instance, the long dash line in Fig. 14 which corresponds to models with fi(H*) 
= 0.4, a clearly excesive value.] 

The solid lines in Figure 14 correspond to models with a small internal dust 
content (fi(H*) =0.1). The dust opacity of the external layer takes on the values of 
Tv(H°) = 0.5 (leftmost solid line) and Ty(if°) = 2.4 (rightmost solid line). From the 
position of the near-case B objects on the left, we infer that their NLR is probably 
less deeply embedded in dust [tv(H°) ~ 2-3]. The fraction of the uncovered area 
runs from 30% to as little as 3.0%-0.3%. Overall, the models bracket quite well 
the region occupied by the Seyfert 2 observations. 

2.4.4 Dust and the UV continuum 

It is not clear how representative of Seyfert 2 the sample of KAW3 is since as many 
as five objects of the 12 objects in Figure 9 show evidence of properties assigned to 
Seyfert 1. These include either a very weak BLR either seen directly as in Mrk 348 
and NGC4388, or seen in polarized light as in Mrk 348, Mrk 3, and Mrk463 (also 
NGC4388), or variability of the featureless continnum as in Mrk 477. Unification 
models of Seyferts have been proposed in which Seyfert 2's present a very thick 
torus which hides direct view of the BLR as a result of its inclination (see, for 
instance, Antonucci & Miller 1985; Krolik & Begelman 1988). Our conclusion 
from §2.4.3 is that the ionized gas regions responsible for the emission lines are 
at least partly covered by large amounts of gas and dust with Lya emerging from 
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Figure 14: 

Same observations as in Fig. 9 but compared with the spherically closed geometry cal- 
culations. The ionized sphere is surrounded by an outer shell of neutral gas N^ [with 
dust content n(H°); see eq.[3]) which we express in terms of dust opacity Ty(H°). The 
curves represent the effect of having a hole in this shell which directly exposes part of the 
Stromgren sphere and allow some Lya to escape out. Models with internal dust content 
H(H*)=0.1 and either rv(H°) = 0.5 or 2.5 are shown as solid lines. The long dash line 
correponds to a model with a higher internal dust content in the ionized region with 
n(H*)=0A. Tick marks along each curve represent the percentage of the uncovered area 
reltive to the sphere's projected emitting area (ttB%). 
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the least covered regions. According to this picture, we would expect the nuclear 
continuum to be reddened as well. We will present briefly some arguments which 
favor continuum reddening, a possibility earlier studied by Boisson & Durret (1986) 
(see also Carleton et al. 1987), although the idea of an extinction free continuum 
is more commonly accepted (e.g., F086, KAW3). 

KAW3 performed an important test on the photon budget of the nuclear 
engine by extrapolating the power laws which they had fitted to their IUE spectra 
and deriving the number of ionzing photons which originates from the nucleus 
(assuming an isotropic source). After comparing with the number of recombination 
photons deduced from the Wf3 luminosity (using the dereddened rlf3 flux and 
assuming a covering factor of unity), they found that most objects showed a clear 
deficit of ionizing photons {N^/N°^ > 1), which they concluded as favoring the 
"occultation/reflection picture" of Antonucci & Miller (1985). According to this 
picture, the ionizing photons stream out freely only along the axis of the torus axis 
and are sufncently numerous to account for the luminosity of the HI recombination 
lines (although the continuum observed in our direction [off-axis] is too weak to 
account for them). 

Although the "occultation/reflection picture" is supported by other indepen- 
dent arguments, there remains, however, an unresolved paradox (see Binette, Fos- 
bury & Parker 1993): if H[3 of the large NLR is so much reddened (and considering 
a closed geometry for the NLR), why is the nuclear continuum of comparatively 
pointlike size not reddened at all? Interestingly, using other IUE Seyfert observa- 
tions, Boisson & Durret (1986) as well as Carleton et al. (1987) concluded that 
there is no deficit of ionizing photons if dereddening of the nuclear continuum is 
taken into account. The imporant argument against continuum reddening given 
by KAW3 is that the 2175 A (10 1514 Hz) dust absorption feature was not present 
in their featureless continua. Although this seems to be clearly the case for objects 
like NGC 1068, we are not strongly convinced that this feature is totally absent 
from some of their spectra given the noise level present. In particular, objects like 
Mrk 3 and possibly NGC 4388 show plausible hints of such a feature. Furtermore, 
the abscence of the 2175 A feature is not so compelling if we consider the possiblity 
that the properties of the dust in AGNs might differ from that in the Galaxy. After 
all, we know that dust grain properties depend strongly upon the environment in 
which they form (Whittet 1992). 

One way of addressing the issue of extinction of the nuclear "nonstellar" con- 
tinuum is to look at possible hints of correlation of the inferred deficit of ionziing 
photons with the IUE spectral index as well as with the amount of reddening in- 
ferred from the Balmer decrement. We present in Figure 15 the results of such 
an attempt based upon the KAW3 data set. Note that we plotted in abscissa the 
quotient Nf^/N°^ which differ from KAW3 in that it is based on the observed 
lip not corrected for reddening (since we want to model all the observed quantities 
as a function of extinction). Although the correlations of Figures 15a and 156 are 
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Figure 15: 

In panel (a), the power-law index a vs. the ratio N°^/N°^. The spectral index a 
corresponds to the power-law fit by KAW3 of the observed continuum energy distribution 
between 2600 A and 1300 A . In panel (b), the observed Balmer decrement is plotted as 
a function of N°^/N°^. The ratio N*£/Nj*Z differs from KAW3 in that the observed 
value of the H/3 flux is used to derive the number of recombination photons instead of 
using the reddening-corrected H/? flux. The solid line is a model representing the effect of 
increasing dust extinction on an intrinsic power-law ionizing continuum of index a = —0.3 
assuming a covering factor of 15% (panel [a]) and an intrinsic Balmer decrement of 3.1 
(panel [6]). The apparent index a of the calculated transmited spectrum is derived by 
simply using the flux ratios at 2450 A and 1310 A . The star represents the locus of a 
model with a shell dust extinction of ry=1.3. 
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not statistiaclly significant, given the magnitude of the errors, the data is not in- 
consistent with a simple reddening model (solid line) described below and applied 
to an intrinsic power law of index a = —0.3 (assuming the gas covering factor and 
the continuum's shape to be intrinsically the same in all the objects). 

To obtain the theoretical results represented by the solid line in Figure 15, we 
illuminate isotropically a slab of dust with opacity Ty with an a=-0.3 power-law 
continuum and a one-parameter sequence of slab models in Ty. For eah model, we 
compute the emergent reddened a (defined by the flux ratios at 2450 A and 1310 A ) 
in the calculated transmitted dust extinction (t v on the observed start- 
ing with an unreddend value (N^ /N°^) =0.15, corresponding to a gas covering 
factor of 15% (varying the covering factor simply slides the curve horizontally). In 
figure 126 we present the intrinsic ratio of 3.1. The standard extinction curve has 
been assumed in our model. Details on the transfer solution for the continuum 
and for the lines are given in Appendix C. 

We have assumed the impinging flux to be isotropic in order to approximate 
a shell layer of dust illuminated by a pointlike central source; a more appropriate 
treatment which distinguishes the isotropic optical lines from the radio continuum 
flux will be the subject of subsequent work. Because both the continuum and 
Wf3 are reddened by the same amount of dust, the equivalent widht of Wf3 is 
independent of Ty in our slab treatment. Although the index a = —0.3 is rather 
flat, it is not at odds with the range -0.08 to +0.11 denoted by Shuder (1981) (for 
the case of a covering factor of 10%) in his study of the proportionality between 
continuum and Ha luminosity in various clases of AGNs (see also Binette, Fosbury 
& Parker 1993). 

Comparing Figures 15a and 156, the slope of the model in the case of the 
Balmer decrement panel seems too low. This could be a consequence of an in- 
homogeneous covering of the nuclear dust as proposed in §2.4.3 for the HI lines 
(Figs. 13 and 14). Since Njg/NgZ and a depend mostly on the UV part of the 
spectrum which is very sensitive to small amounts of dust, these quantities are ex- 
pected to be strongly weighted towards the lest obscured regions, while the optical 
Balmer lines may merge from more obscured regions. 

Despite the recognized validity of the "occultation/reflection picture" in ex- 
plaining cones of ionization and hidden polarized BLR in some Seyfert 2's, it 
remains important to test further whether the nuclear continuum is as reddening 
free as generally assumed and, therefore, whether part of the apparent deficit in 
the ionizing radiation does not simply result form dust extinction. 
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Chapter 3 



Calcium depletion and the 
presence of dust in large scale 
nebulosities in radio galaxies (I). 

Villar-Martm & Binette 1995, A&A, in press 

Abstract 

We show that the study of the Calcium depletion is a valid and highly sensitive 
method for investigating the chemical and physical history of the very extended 
ionized nebulae seen around radio galaxies (EELR), massive ellipticals and 'cooling 
flow' galaxies. By observing the near IR spectrum of nebular regions characterized 
by low excitation emission lines (LINER-like), we can use the intensity of the 
[CaII]AA7291, 7324A doublet -relative to other lines, like Ha- to infer the amount 
of Calcium depletion onto dust grains. The presence of dust in these objects - 
which does not necessarily result in a measurable level of extinction- would favour 
a 'galactic debris' rather than a 'cooling flow' origin for the emitting gas. Before 
applying such test to our data, we study four possible alternative mechanisms to 
dust depletion and which could have explained the absence of the [Call] lines: a) 
ionization of Ca + from its metastable level, b) thermal ionization of Ca + , c) a high 
ionization parameter and/or a harder ionizing continuum than usually assumed 
and d) matter bounded models associated to a hard ionizing continuum. We show 
that none of these alternative mechanisms explain the absence of the [Call] lines, 
except possibly for the highly ionized EELR where a high ionization parameter 
is required combined with a soft power law. We thus conclude that for the other 
low excitation emission regions (cooling flows, liners, low excitation EELR), the 
absence of the Call lines must be due to the depletion of Calcium onto dust grains. 
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3.1 Introduction 

The study of the interstellar medium (ISM) of external galaxies provides important 
information about the global kinematic (inflow, outflow) of such gas, its chemical 
composition and the implied star formation history, its mass distribution, etc. This 
gas forms a vital part of the record of the formation of the parent galaxy, and the 
evolutionary processes involved. 

How can this material be studied in details? One way is to have it illumi- 
nated or excited by a powerful AGN, giving rise to the phenomenon of extended 
emission line region (EELR). The drawback of course is that it only allows us to 
look at a restricted class of galaxy. The large scale EELR phenomenon is observed 
in a majority of the most powerful radiogalaxies with EELR extending to radial 
distances of up to 100 kpc from the nuclei (Tadhunter 1986; Baum et al. 1988), 
much larger radii than the stellar population distribution of the parent galaxy. 
The morphologies and kinematics of such regions cover the full range from regu- 
lar disc/ring systems to chaotic systems for which no pattern can be discerned. 
Their spectra show strong emission lines, covering a wide range in ionization. It 
is generally accepted that the EELR are ionized by some mechanism connected 
with the nuclear activity, but there is no full consensus on the excitation mech- 
anism. Some objects show evidences for an interaction between the radio jets, 
which transport energy to the outer radio lobes, and the gas in the outer region. 
Maybe this interaction is responsible of the excitation of the gas through some 
kind of shocks (Sutherland, Bicknell & Dopita 1993). For other objects, that do 
not show evident spatial coincidence between the radio structures and the EELR, 
the excitation might be due to direct photoionization by the nuclear ionizing radi- 
ation field (Robinson et al. 1987: hereafter RBFT87). A reduced scale version of 
the EELR is the one observed in many Seyfert galaxies (Haniff, Ward & Wilson 
1988) where the ionized gas may extend up to a few kpc although it is brightest 
within the central 100-300pc. The EELR morphology tends to be conical (Wilson 
& Tsvetanov 1994). Its detailed observation is complicated by the presence of the 
very luminous stellar background of the bulge. In normal bright ellipticals, quite 
weak extended nebulosities of low excitation are a common phenomenon (Buson 
et al. 1993, Goodfroij 1994). The very large scale gas around radio-galaxies, on 
which we focus here, presents the advantage that the lines are observed against 
the sky background rather than against the bright parent galaxy bulge. 

The origin of the gas making up these large scale extended nebulosities re- 
mains unknown. Furthermore, the distinction between these and the filamentary 
nebulae seen in some clusters and around some massive ellipticals and often iden- 
tified with cooling flows, is unclear. We do know, however, that this gas in every 
case is chemically enriched as compared to primordial gas. Emission line analyses 
(RBFT87) show common element abundances to be within a factor of a few (<) 
of Solar and also to be rather uniform over all the objects observed. 

The two most likely explanations for the origin of the material are: 
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a) debris from recent tidal interactions and mergers 

b) gas cooling from the hot (~ virial) phase which from X-ray observations 
(Forman, Jones & Tucker 1985) has been shown to exist around massive ellipticals 
and inside galaxy clusters 

The main arguments behind these explanations are the following: 

a) Heckman et al. (1986) showed that a large fraction of powerful radiogalax- 
ies have morphological features -shells, tails, loops, etc- similar to those produced 
in numerical simulations of galaxy interactions (e.g. Toomre and Toomre 1972, 
Quinn 1984). This could indicate that the activity has been triggered either be- 
cause fresh gas has been accreted from outside or because preexisting gas in the 
galaxy has been caused to collapse to the core as a result of the interaction. This 
interaction scenario is also supported by observations of a few nearby radio galaxies 
which, apart from morphological peculiarities, show large misalignments between 
the stellar and the gaseous rotation axes, indicative of an external origin for the 
gas. 

b) Hot X-ray corona (T > 10 4 K) are a common feature of bright early-type 
galaxies. Within some critical radius, radiative cooling becomes important, leading 
to the development of the so-called 'cooling flow' hypothesis (Nulsen, Stewart & 
Fabian 1984, Thomas et al. 1986) Eventually, condensations or filaments could be 
formed, dense and cool enough to radiate detectable optical emission lines. Most 
powerful radio galaxies are too distant for the characteristic X-ray emission to be 
currently detectable, but it is quite plausible that the EELR gas has condensed 
out of a surrounding cooling flow. 

Distinguishing between these alternative hypotheses has been attempted using 
gas kinematic measurements (Tadhunter, Fosbury & Quinn 1989) which show that 
the radio galaxy EELR generally have a high specific angular momentum which is 
difficult to reconcile with the cooling flow picture. An alternative approach is to 
look for the presence of dust associated with the EELR gas. If the gas has cooled 
directly from a hot phase, there will have been no opportunity for dust to form, 
according to the standard cooling flow theory. Any dust introduced from galaxies 
into the hot intracluster medium will be rapidly destroyed (Draine & Salpeter 
1979). We would need a mechanism to produce this dust once the gas has cooled 
down (Fabian, Johnstone & Daines 1994). If, on the other hand, the material has 
fallen in during a merger, the dust /gas ratio is expected to have a value appropriate 
to the gas chemical composition found in normal galaxies. 

Determining the presence or absence of dust is important because: 

-of its relevance in deriving chemical composition which takes into account the 
effects of depletion of metals unto dust grains 

-of the implication for the star formation history and the ISM evolution: when 
does dust form? (and so stars?) 

-of the effects of dust on the apparent morphology of continuum and line 
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features: pure absorption (reddening) and/or scattering (blueing/polarization) 

We discuss below how the absence of the forbidden [Call] lines can be used to 
infer the presence of dust mixed in with the emission gas. We assess in detail all the 
most plausible alternative explanations to that of internal dust for explaining the 
absence of [Call] lines. As no acceptable alternative solution is found, we conclude 
in favour of the validity of the method initially proposed by Ferland (1993). We 
however adapt and optimize the [Call] dust detection technique to the context of 
the EELR studies in which we are involved. The observational results and their 
interpretation will be presented in a subsequent paper. 

3.2 Outline of the method 

To investigate how the forbidden Calcium doublet of [Call] in the infrared is af- 
fected by the physical condition encountered in photoionized plasma, we outline 
first the computer code which we used and then proceed to illustrate how the [Call] 
lines might be used to infer dust in the ISM of galaxies and thus our interest in 
securing these conclusions by closing up the possibility of alternative interpreta- 
tions. 

To compute the emission lines used in our study, we have used the multipur- 
pose photoionization-shock code MAPPINGS described in Chapter 2. The Ca + 
ion is considered as a five level atom. Its structure is shown in Fig. 16. 

3.2.1 The input parameters 

We now proceed to define the parameters employed in the calculations. Most of 
these are derived from our observational knowledge of EELR although we recognize 
that it is incomplete and can be biased by particular diagnostic tools which are 
employed. 

Metallicity 

For definiteness, we adopt a set of solar abundances (Anders & Grevesse 1989) 
for the trace elements, in line with the finding of RBFT87 who indicated values 
for the EELR not radically different from solar. The solar abundance of (not 
depleted) atomic Calcium is 2.2 x 1CT 6 by number relative to hydrogen. Depletion 
into interstellar dust grains is known to reduce this abundance in the local ISM 
by a factor of ~ 5000 (Whittet 1992) so that if depletion was taking place in the 
partially ionized EELR plasma, even a small dust-to-gas ratio might be sufficient 
to eliminate any detectable trace of atomic Calcium. 
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Figure 16: 

Energy level diagram for Call showing the radiative transitions among the five levels 
considered here. Permitted transitions are represented with dotted lines and forbidden 
transitions with dashed lines. Fl and F2 are the forbidden doublet lines considered in 
this work. 



Geometry 

It's been accepted for a long time that the EELR is formed by individual clouds that 
are ionized by the central source. In most models presented here, each emitting 
cloud is considered (see section 2.2.1) like a radiation bounded slab (optically 
thick to the Lyman continuum) which comprises a fully ionized region, closer to 
the illuminated face and responsible for the high ionization lines and a partially 
ionized zone (PIZ), where the low ionization lines like Call are emitted. Matter 
bounded clouds will be considered in Sect. 3.3.4. As we want to concentrate on 
the depletion phenomena, we will simplify the calculations by not considering any 
neutral region (see Fig. 4, Chapter 2) beyond the PIZ which would contain dust, 
the effect of which would be to cause additional extinction for an hypothetical and 
unfavorably placed observer looking from the back of the slab. We have verified 
that the results reached here using line ratios of similar wavelengths are not altered 
by the presence or not of this neutral absorbing zone. 



Gas density 

We adopt a representative density for the EELR clouds of 300 cm -3 . Typical elec- 
tron densities values measured vary between a few tens (or less) to a few hundreds. 
The densities derived from forbidden line ratios might not apply to every subre- 
gion, but they are sufficiently low to consider the low density limit as a generally 
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valid and very good approximation to the physical conditions affecting the Call 
lines of large scale EELR. In the low density regime, the effects of density varia- 
tions on the emission line ratios which are considered in this paper are very small 
by comparison to the effects of other parameters like the ionization parameter or 
the hardness of the ionizing continuum. 

The calculations consider the gas pressure to be constant within the cloud 
(isobaric models) with the density behaviour modulated with depth into the cloud 
by the behaviour of the temperature and by the ionization fraction of the gas. 

The ionizing continuum and the ionization parameter U 

We implicitly assume that the dominant ionization mechanism of the ionized gas 
in radio galaxies is photoionization. It has been shown by RBFT87 that a hard 
continuum extending well down into the soft X-ray region, despite some discrep- 
ancies with the observed spectra, can be considered to reproduce generally well 
the measured line ratios. This continuum can take its source in the active nucleus 
or be locally generated (e.g. fast shocks: Binette, Dopita, Tuohy 1985). We have 
considered power law (PL) distributions of various values of a (F v oc v +a ) as well 
as hot blackbodies (BB) of temperature T^ ps 10 5 fT in order to study the effects 
that hardness has on our conclusions. 

The ionization parameter, a measure of the excitation level of the ionized gas, 
was defined in section 2.2.1. 



quotient between the density of impinging ionizing photons and the density 
of the gas cloud. 

We find that the parameters having the strongest effect on the line spectrum 
are the ionization parameter U and the mean ionizing photon energy {i.e., hardness 
of the continuum). 

3.2.2 How to detect dust 

How can we detect dust within the gas associated with EELR? We here summarize 
different techniques and compare their sensitivities. 

a) Reddenning: given the nature of the gas distribution and the fact that it is 
ionized externally (unlike HII regions which are internally excited), the extinction 
may not necessarily be sufficiently high to be easily detected using optical obser- 
vations. Indeed, the EELR ratios studied by RBFT87 show little or no reddening. 
Furthermore a small enhancement of the Balmer decrement over recombination 
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case B might be interpreted as resulting from collisional excitation rather than 
from reddening. 

b) Scattering: polarization measurements of high redshift radio galaxies (di 
Serego Alighieri et al. 1989, Januzi & Elston 1991, Tadhunter et al. 1992) show 
conclusive evidence for scattered nuclear light over large volumes and, although 
there are reasons to believe that the scattering medium is dust, it is difficult to rule 
out entirely Thomson scattering by electrons. Detailed studies of the low redshift 
galaxy PKS2152-69 do, however, show polarized continuum radiation from highly 
excited extranuclear gas cloud with an energy distribution which is so blue that it 
must arise from dust scattering (Fosbury et al. 1990). 

c) Infra-red thermal dust emission: the IRAS satellite has shown that many 
galaxies radiate significant fractions of their energy in the far infrared sprectral 
region. Significant masses of dust at temperatures of around 40K are responsible 
for this radiation at wavelengths of 60fim and beyond. In many cases the FIR 
spectral energy distributions is still rising at 100/im, out of the spectral range 
detectable by IRAS. The cool dust can only be detected at millimetre and submil- 
limetre wavelengths. A strong limitation of this technique is the very poor spatial 
resolution of the IRAS satellite. Groundbased studies of the far IR emission of 
galaxies in the sub-mm range also exist (e.g. Clements, Andreani & Chase 1993), 
but still with poor spatial resolution. Although it is possible to infer masses and 
temperatures of the warm and cool dust components (dependent on models), the 
spatial distribution of the dust is not known. 

d) Indirect effects on the line spectrum. There are several ways this can 
happen: effects of dust on the gas temperature (photoionization of dust grains may 
raise the temperature of the plasma), effects of dust on the ionization structure 
(dust grains selectively absorb ionizing photons of lower energies), and influence on 
apparent chemical composition via the depletion of refractory elements onto dust 
grains. The first effect does not provide a unique interpretation for the unusually 
high temperatures seen in some EELR (Tadhunter, Robinson & Morganti 1989) 
while the second effect cannot be discriminated against reliably since even dust- 
free photoionization models are still too uncertain to be used as absolute reference 
point. For these reasons, the last effect is the only clearly promising one and is 
looked into details below. 

Calcium is very sensitive to the presence of dust as it is always found to be de- 
pleted in the interstellar medium (e.g. Crinklaw, Federman & Joseph 1994). Pho- 
toionization calculations appropriate to LINERS -a hard ionizing spectrum with 
a relatively low ionization parameter- invariably predict the [CaII]AA7291, 7324A 
doublet (4s 2 S-3d 2 D) to be very strong (Ferland 1993). These two forbidden lines 
(hereafter Fl = A729lA F2 = A7324A ) have a high critical density ~ 10 6 cm" 3 . 
The latter line, 7324A is the weakest of the doublet and is furthermore blended 
with the [OII]A7325A multiplet. The other line, A729lA lies some 30A shortward 
of [Oil] and is therefore straightforward to isolate given reasonable spectral reso- 
lution. The fact that any of these doublet lines are generally not seen in LINERS 
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Figure 17: 

MAPPINGS photoionization models showing the strenght of [CaII]A7291 relative to Ha 
(left panel) and [ArIII]A7135 (right panel) as a function of the ionization parameter U. 
The models assume solar abundances but with different amounts of dust-to-gas ratios 
H expressed in solar neighborhood units. The gas phase abundances of Ca/H for the 
H = 0, 0.1 and 1.0 models are 2.2 x 10~ 6 (no depletion), 9.2 x 10" 8 and 4.4 x 10~ 10 
(Whittet 1992), respectively, after taking into account depletion. The hydrogen density 
is 300 cm" 3 and the ionizing power law index is a = —1.4 

but are so in some novae when their envelopes have reached the appropriate ion- 
ization level can be interpreted as evidence of Calcium depletion onto dust grains 
in the former objects. Since EELR which are the subject of our investigation are 
often seen in the ionization parameter region of the line ratio diagnostic diagrams 
occupied by LINERS (10~ 4 < U < 10~ 3 ), we can similarly use the [Call] doublet 
measurements to infer whether or not there is depletion taking place in EELRs 
and thus conclude whether dust is also mixed with the gas as is thought to be the 
case in LINERs. 

The above arguments have already been used for several objects with 'cooling 
flow' filaments by Donahue & Voit (1993) to infer the presence of dust mixed with 
the ionized gas. Ferland (1993) has shown the great sensitivity of this method 
to the presence of dust under NLR conditions. We now show it to be also the 
case under EELR conditions. We present two diagnostic diagrams in Fig. 17, the 
ratio [CaII]/Ha (7291/6563) and the ratio [CaII]/[ArIII] (7291/7135) as a function 
of the ionization parameter U. The variable parameter distinguishing the three 
different sequences in U is fi, the dust-to-gas ratio of the plasma expressed in units 
of the solar neighborhood dust-to-gas ratio. The dramatic difference in line ratios 
between the grain depleted Ca/H (/x > 0) and the undepleted case (/i = 0), shows 
the sensitivity of this method, particularly for low values of U. 
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3.3 Possible alternative explanations to deple- 
tion. 

Before we carry the conclusions of the current analysis to the interpretation of our 
observations (Chapter 4), we report first on our effort in investigating other possi- 
ble alternative mechanisms to dust depletion. If the warm Ca + region predicted by 
standard models does exist, then Calcium depletion becomes the only reasonable 
explanation for the non detection of the doublet lines. What we consider in this 
section is the possible NON EXISTENCE of the emitting [Call] region by inves- 
tigating different mechanisms which could eliminate it. During our investigation, 
we require however that successful models do not result in important discrepancies 
with other observed line ratios. The mechanisms we have considered to eliminate 
the [Call] region are 

- Ionization of Ca + * by Lja and soft continuum photons from the metastable 
level of Ca + 

- Thermal (collisional) ionization of Ca + 

- Photoionization with a much harder continuum or a much higher U than 
usually assumed 

3.3.1 Ionization by Lya and soft continuum photons. 

Wyse (1941) proposed that the ionization of Call from the metastable level by 
Lja photons, could explain the fact that the IR lines of Call at 8498, 8542 and 
8662A appear in emission near the maximum phase of Me variables, whereas the 
H and K lines only occur in absorption. Trapped Lja photons could also play a 
part in ionizing metastable Ca + * as suggested by Wallerstein et al. (1986). 

We investigate here if this process is important under the conditions found 
in EELR clouds. In order to do this, we add two terms to the ionization equi- 
librium equation of Call. One which considers photoionization of excited Ca + * by 
the impinging UV continuum. The other is photoionization of excited Ca + * by 
the nebular Lya photons. The statistical equilibrium equations give the relative 
population of the mestastable level, which turns out to be, under EELR condi- 
tions, ^ ~ 10~ 7 , being the density of Ca + ions in the ground level. With 
such a negligible population, the density of ionizing photons must be very high to 
increase the ionization rate to a non negligible level as compared to the ground 
state ionization rate. A simple estimate presented in Appendix D demonstrates 
this level to be out of reach. 

In summary, the very diluted radiation fields and the low densities appropriate 
to the EELR implies an extremely small population for the excited levels which 
prevents the ionization of Ca + * by Lja and soft continuum photons from being of 
any significance. 
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3.3.2 Thermal ionization of Ca + . 

We investigate here the possibility of collisional ionization by thermal electrons 
of Ca + to Ca ++ , a process which is important when the electronic temperature 
becomes higher than 20000K (Jordan 1969). In order to establish a comparison 
in U, we have considered two extreme cases in our calculations, logC/ = —4 & — 
2. To illustrate how a much harder continuum will result in much higher gas 
temperatures, we also use two different PL of index a = —1.4 and —0.4. Note that 
such a hard continuum as a = —0.4 is probably quite unrealistic. However, our 
intention here is simply to test whether very high temperatures can be achieved 
with photoionization models and specifically near the Ca + region. The results are 
shown in Fig. 18 as a function of depth in the photoionized slab. Of the eight 
plots, the four upper ones correspond to a = —1.4 while the four at the bottom to 
a = —0.4. The four plots on the right, have logU= —2, and the four on the left, 
logU= —4. Two plots therefore are shown for each pair of [U, a] values: one is the 
the temperature T4 in units of 10000K and the one immediately underneath is the 
intensity of Fl (erg.s~ 1 .cm~ 2 ), both as a function of the depth in units of 10 20 cm 
into the slab. These plots allow us to see the corresponding electron temperature 
to the position where the bulk of the [CalljFl emission takes place. 

For the traditional PL of index a = —1.4, the electron temperatures are not 
anywhere near high enough for the process of thermal ionization to be relevant. 
Harder continua and increasing U values do produce higher temperatures, but, 
unless U are unrealistically high (even logU= —2 is not enough), the gas is never 
hot enough. Thus, realistic photoionization models are not able to heat the gas 
sufficiently to thermally ionize Ca + . We might conjecture that there could be an 
additional heating source like shocks which could raise the temperature of the gas. 
This would be an interesting point for further investigation. 

3.3.3 Effects on Fl of varying U and the continuum hard- 
ness. 

We now investigate how the Fl line might become undetectable by simply vary- 
ing arbitrarily the ionization parameter U or the continuum hardness. Fig. 19 
shows six diagnostic diagrams with the abscissae always representing the line ratio 
[OI]A6300/[OIII]A5007. The [OI]/[OIII] ratio monotonically increases with de- 
creasing gas excitation (i.e., with decreasing U) and is therefore a good measure of 
the excitation level of the gas. Each diagram shows in the ordinate a different line 
ratio which can be related to a given gas property. [OI]A6300/Ha, for instance, 
might measure the hardness of the continuum. In the last diagram, the ordinate 
corresponds to the quotient F1/[OI]A6300. The three sequences of models shown 
in each diagram differ by the slope of the power law which takes on the values of 
a = —1.4, — 1 and —0.4. The values of Log U covered by each curve is in the range 
—4 to —1. Our aim is to look for models which can decrease the Fl intensity below 
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Figure 18: 

Distribution of temperature and integrated intensity of Fl (in ergs s _1 cm~ 2 ) line across 
an EELR cloud. Dotted lines (top diagrams) correspond to a power law of index —1.4; 
solid lines (bottom diagrams) to a power law of index —0.4. Each temperature diagram 
must be compared with the one just below, to see the temperature that correspond to 
each point in the nebula. 
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Figure 19: 

The curves correspond to ionization parameter U model sequences for different ionizing 

continua: power laws (small circles used as fiducial marks) of index a = —1.4 ( ), — 1 

(- - ), -0.4 ( ) and a black body continuum of T bb = 1.2 x 10 5 K (*-*-*). The 

region represented by vertical dotted lines in the last diagram shows the area where Fl 
would fall under the detection limit defined in Sect. 3.3.3. All the models within this 
region are successful in accounting for the absence of the [Call] lines without requiring 
depletion by dust. Line ratios of extended emission line regions observed by RBFT87 
(uncorrected for reddening) are shown with solid squares. 
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the detection limit. Let's look at how we might define a practical detection limit. 
The open squares in the diagrams of Fig. 19 represent line ratios measured by 
RBFT87 in several EELR. The faintest line they measure is typically HeIIA4686. 
The mean ratio of HeIIA4686/[OI]A6300 observed is 10" a4 = 0.4 for the large scale 
nebulosities. We establish our 'artificial' detection limit in the following way: since 
Hell is one of the weakest line successfully measured by RBFT87, we will assume 
that any line fainter than 0.4 below the [OI]A6300 flux is not detectable. In the 
last diagram, the region where Fl falls below this detection limit is shown by a 
dash line grid. Any model found in this area is deemed successful in explaining 
the non-detection of Fl without requiring depletion. 

We see that for the standard PL (a = —1.4), only models with high U (log 
U> —2) decrease log(Fl/[OI]) below -0.4. These models, as we can see in the 
diagrams, would therefore be valid only for the high excitation EELR, but not for 
LINERs, cooling flow filaments, or EELR of low and intermediate excitation. On 
the other hand, increasing the hardness of the continuum (flatter power laws), helps 
Fl to get fainter with respect to [OI]A6300, but the discrepancies with observed 
line ratios in other diagrams become totally unacceptable (see top two diagrams 
of Fig. 19). 

It is interesting to compare a BB sequence (1.2 x 10 5 K) with the canonical 
PL sequence a = —1.4. We see that both ionizing continua reproduce rather well 
the observed line ratios as was earlier shown by RBFT87. From these line ratios 
alone, there are no reasons to favour power laws over hot blackbodies. A similar 
conclusion was reached by Binette, Robinson and Courvoisier (1988) for the mean 
NLR spectrum of Seyferts. 

A BB produces a much stronger Fl compared to [OI]A6300 than any of the 
power laws considered. One reason for this is that the fraction O /O + in the PIZ 
is completely controlled by the charge exchange reactions of 0° and + with H + 
and H°, respectively, and not by direct photoionization of 0°. This is not so for 
Ca + /Ca° which is free to respond to the different amount of hard photons (the 
only one to make it to the PIZ) available in a PL or a BB. We conclude from the 
last diagnostic diagram that the BB models could not explain the absence of Fl 
from the observed spectra without invoking depletion. 

3.3.4 Effects on Fl of truncating clouds 

We showed in the previous section that a harder continuum can potentially bring 
Fl under the detection limit but results in important discrepancies with the low 
excitation lines. The reason is that the partially ionized zone (PIZ) where most 
of the low excitation lines are generated gets larger and larger with increasing 
hardness of the continuum. If the clouds were truncated, the smaller PIZ would 
generate weaker low ionization lines, thus improving the overall fit. 

We have investigated models with a = —0.4 which were truncated at a depth 
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Figure 20: 

Same diagrams as in Fig. 19 but for truncated gas clouds and only for hard power laws 
(a = —0.4). The open circles (connected by the dashed line) correspond to matter 
bounded clouds (5.2 x 10~ 5 < U < 5.2 x 10~ 4 ) truncated at a depth which produce 
[OIII]/H/3 = 10 (see text). The open triangles (connected by the solid line) are similar 
models truncated at a depth such that F1/[OI]<0.4 (3.7 x 1(T 3 < U < 2.7 x 10~ 2 ). 
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which satisfies a given criterion based on a specific line ratio. This has been done 
in two ways: 

1) The criterion in this case is to truncate the calculations when OIII/H/9 has 
reached a value of 10 which is the typical ratio for the high excitation EELR. 

The sequence of models shown in Fig. 20 are separated by a factor of 0.14dex in 
U. Altough there are still discrepancies, there is a notable improvement compared 
to the radiation bounded models of Fig. 19. The predicted line ratios are now 
located closer to the observed data (same scale as in Fig. 19). We conclude that 
models with a hard ionizing continuum must be matter bounded in order to fit 
acceptably most observed ratios. 

What happens now with Fl/[OI]? In the last diagram of Fig 20, we see that all 
these models produce Fl above the detection limits and cannot therefore explain 
the non detection of Calcium. 

2) The second method consists on imposing that the models produce Fl under 
the chosen detection limit (F1/[OI]=0.4) and check if such models agree with the 
position of the observed line ratios in the rest of the diagrams. 

The models in the sequence which satisfy this criterion are found in the range 
3.7 x 10~ 3 < U < 2.7 x 10~ 2 . They are represented in Fig. 20 as open triangles 
connected by a solid line. As we see in the top left diagram, the ratio [OIII]/H/9 is 
not any more defining a simple trend with excitation (which is represented by the 
ratio [OI]/[OIII]). Furthermore [01] /Ha remains a discrepant ratio as in Fig. 19. 
So although it is in principle possible to satisfy the F1/[OI]=0.4 criterion with a 
hard power law, the truncation must be done at a specific yet ad hoc depth and 
furthermore the previous trend of excitation with U has disappeared. 

Without rejecting the possibility of a more complicate mixture of matter and 
radiation bounded models, we believe that simply truncating clouds does not con- 
vincingly solve the problem of the weakness of the Call doublet and dust depletion 
remains the most likely interpretation. 

It is interesting to note that truncated clouds adjust better the HeII/H/5 ratio 
(bottom left diagram) as proposed before by Morganti et al. (1991) and Viegas 
and Prieto (1992). 

3.4 Conclusions 

This work is based on the method proposed by Ferland (1993) to investigate the 
presence of dust mixed with the gas of the Narrow Line Region of active galax- 
ies. Because photoionization models predict remarkably strong forbidden lines 
[CaII]AA7291,7324A assuming reasonable abundances of atomic Ca, the basic idea 
is to infer a systematic depletion of Calcium onto dust grains whenever the infrared 
[Call] lines are observed very weak or undetected. This test of the dust content 
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was applied to cooling flow filaments by Donahue & Voit (1993) who concluded on 
the presence of dust. 

We have shown here that this sensitive method is also applicable to the con- 
ditions found in the EELR of radiogalaxies. In order to make more secure any 
inference about the presence of dust based on [Call] lines, we have investigated 
alternative explanations for their absence: ionization of Ca + * to Ca ++ by Lja 
photons and soft continuum photons from the metastable level of Ca + , thermal 
ionization of Ca + , ionization of Ca + due to either a very high U value (ionization 
bounded case) or to a hard continuum (with truncated clouds). Except for the 
highly excited EELR which might not possess any Ca + region due to their high 
ionization level, the results are negative: none of the alternative mechanisms or 
models studied can explain the absence of the [Call] lines without dust depletion. 

Our conclusion is that the dust content test appears generally valid for the 
EELR of radiogalaxies (unless the excitation level of the gas is extremely high). 
This will allow us to make important conclusions about the origin of such gas, 
discriminating between galactic debris and the standard cooling flow theory. 

A BB ionizing continuum characterized by a temperature of 1.2 x 10 5 K can 
reproduce the observed line ratios of the EELR at least as well as a PL of index 
a = —1.4. On the other hand, the Fl/[OI] from a BB is higher than that of a PL 
so the case in favour of depletion is even stronger. 

In a follow up paper, we will present long slit spectra of EELR, cooling flow 
filaments and Seyfert 2 NLR, all taken in the region of the [Call] doublet. The 
goal will be to apply the test of the Calcium depletion described above in order 
to conclude whether or not the gas in these nebulosities is mixed with dust. This 
will be our starting point for deciphering the origin of the emitting gas. 
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Chapter 4 



Calcium depletion and the 
presence of dust in large scale 
nebulosities in radio galaxies (II). 

Villar-Martm & Binette 1995, A&A, submitted 

Abstract 

We investigate here the origin of the gas observed in extended emission line 
regions surrounding AGNs. We use the technique of calcium depletion as a test 
to prove or disprove the existence of dust in such a gas in order to discriminate 
between two main theories: (1) a cooling process from a hotter X-ray emitting 
phase surrounding the galaxy, (2) merging or tidal interaction between two or 
more components. We have obtained long slit spectroscopy of a sample of objects 
representative of different galaxy types although our main interest focus on radio 
galaxies. The spectral range always includes the [CaII]AA7291,7324 doublet. The 
faintness or absence of such lines is interpreted as due to the depletion of calcium 
onto the dust grains and, therefore, is a proof of the existence of dust mixed with 
the gas in the EELRs. 
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4.1 Introduction 

In Chapter 3, we have demonstrated that the test of proof of the existence of inter- 
nal dust based on the depletion of calcium is valid (and very sensitive) under the 
conditions found in EELR studied here. We studied in detail all the most plausible 
alternative mechanisms to that of internal dust for explaining the absence of [Call] 
lines. No acceptable alternative solution was found and we concluded in favour 
of the validity of the method. The observational results and their interpretation 
are presented in this chapter. We have applied this test to a sample of objects 
of different types: radiogalaxies with extended emission line regions, cooling flows 
(hereafter CF) and starburst galaxies. 

We describe the observations, data reduction and analysis of the spectra in 
§4.2. In §4.3 we present the results of the comparison between measurements and 
model predictions. The implications on the origin of the gas are discussed in §4.4. 
The conclusions compose section §4.5. 

4.2 Observations and data reduction 

The observations were carried out on the nights 21-23 August on 1993. All the 
spectra were obtained at La Silla Observatory, Chile, with the 3.6 m telescope, 
using the EFOSC 1 spectrograph with a CCD detector (TEK#26) of 512 x 512 
pixels 2 of 27 /im 2 . The slit width was 1.5". The grism, R150, with a dispersion of 
120 A /mm, a wavelength bin of 3.3 A /pixel and covering a spectral range of ~ 
6870-8560 A . The observing conditions where photometric. Table 1 gives a log of 
the observations. 



Table 1: Observing Log 



Name 


RA(1950) 


Dcc(1950) 


z 




Comment 


Exp time (s) 


PA 


NGC1052 


02 


38 


37.33 


-08 28 09 





.005 


Liner 


1800 


270° 


NGC6215 


16 


46 


47.0 


-58 54 30 





.005 


Liner 


2700 


240° 


NGC7552 


23 


13 


25.0 


-42 51 24 





.005 


Liner 


1200 


270° 


NGC7714 


23 


33 


40.59 


01 52 42 





.009 


Starburst 


2700 


270° 


PKS1404-267 (nuc) 


14 


04 


38 


-26 46 51 





.021 


RG 


2700 


270° 


PKS1404-267 (5" S) 
















2700 


270° 


PKS2014-55 


20 


14 


06 


-55 48 52 





.061 


RG 


3600 


190° 


PKS2152-69 (nuc) 


21 


52 


58 


-69 55 40 





.028 


RG 


2700 


270° 


PKS2152-69 (cloud) 
















3600 


290° 


PKS2158-380 


21 


58 


17 


-38 00 51 





.033 


RG 


2700 


270° 


PKS2356-61 


23 


56 


29 


-61 11 42 





.096 


RG 


3300 


285° 


PKS2300-18 


23 


00 


23 


-18 57 36 





.129 


RG 


2700 


240° 


2A 0335+096 


03 


35 


52 


09 48 10 





.035 


CF 


5400 


147° 


A2029 


15 


08 


30 


05 57 00 





.077 


CF 


2700 


270° 


A2597 


23 


22 


42 


-12 23 00 





.085 


CF+RG 


2700 


197° 
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a) Basic data reduction 

The reduction of the data was done using standard methods in IRAF. The spectra 
were bias subtracted and divided by a flat-field frame (dome flat-field). Illumina- 
tion corrections along the slit were found to be negligible. 

In general, we obtained three frames for each object and each slit position, 
allowing the direct removal of cosmic ray events. For a given object, all the frames 
corresponding to the same slit position were averaged together. In the cases where 
only two frames were available, the cosmic rays were removed visually, replac- 
ing the affected pixels by the mean of the surrounding region. The spectra were 
calibrated in wavelength using comparison spectra of an HeAr arc taken at the be- 
ginning and end of each night, and additionally before and after each object. The 
wavelength calibration was done very carefully in order to later subtract the sky as 
accurately as possible. The IRAF routine "background" was used to substract the 
contribution of the sky by interpolation of the background detected in windows 
close to and on both sides of the object. With this method the sky features were 
successfully substracted. 

b) Atmospheric extinction 

The spectra were corrected for atmospheric extinction with the aid of mean extinc- 
tion coefficients for La Silla. Molecular absorption bands of O2 (the A band at 7620 
A and the H 2 bands, bands in the region 7100-7450 A and 8100-8400 A ) were 
also evident in the spectra. They were removed separately from all of the observed 
spectra. These absorption bands are composed of many closely spaced absorption 
features, unresolved in our spectra. Many of the features which comprise the A 
band are optically thick and therefore nearly independent of zenith distance, but 
most of those in the H2O bands are not independent of time and zenith distance. 
In principle, the best method to remove these lines is to observe standard stars as 
close as possible in zenith distance and in time to that of the program object to 
obtain the best correction (c.f. Osterbrock, Shaw & Veilleux, 1990) . But as these 
authors pointed out the star closest to the objects's observing time and zenith 
distance is not always the one producing the best correcting spectrum. 

The procedure we used was the following: we observed three different stan- 
dard stars whose frames were reduced in the same way as described before. After 
modeling the atmospheric bands with them, we decided that the best results were 
obtained with Feige 110, a dwarf with no intrinsic absorption features in the spec- 
tral range observed. An 1-dimensional integrated spectrum of the standard star 
was obtained by adding all the light along the spatial direction. A smooth fit 
to the continuum was done and then the original spectrum was divided by this. 
The result was a normalized spectrum keeping the features due to the atmospheric 
absorption. Due to the dependence on zenith distance and time, we had to model 
the bands for each object, i.e. to construct an specific restoring spectrum for each 
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Figure 21: 

Spectra of observed Radiogalaxies, nuclear (first column) and extended emission 
line regions (second and third columns). The expected position of the Fl line is 
indicated with | . Flux is given in units of erg s~ x crrr 2 A~ x . 
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Figure 22: 

Spectra of observed Liners and Seyfert 2 galaxies. The liner spectra have been 
extracted from the whole spatial extension of the object. The expected position 
of the Fl line is indicated with | . Flux is given in units of erg s~ 1 cm~ 2 A~ 1 . 



object, multiplying the normalized spectrum by appropriate factors. 



c) Flux calibration 

The atmospherically corrected spectra of the standard stars were used to obtain 
the flux calibration. For each night we built a mean response curve from the two 
standard stars observed that night. 



d) Template galaxy subtraction. 

We checked that effects of stellar absorption features intrinsic to the galaxy could 
be neglected: the calcium doublet is on top of the raising part of a molecular band 
whose narrow line components are not resolved. Its only effect is to change the 
slope of the continuum under the doublet and we can easily correct for this effect, 
fitting a smooth continuum of certain slope in such region. Ha, a line that has 
been used as a reference to predict the flux in the Fl line, can be underestimated 
due to the underlying Ha absorption line. But in our objects the Ha emission 
is so strong that the difference is less than 5% of the emitted flux. This was 
computed by comparing the EW of the emission line in our objects and the EW 
of the absorption feature in a template elliptical galaxy. 
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Figure 23: 

Spectra ot objects with properties associated with cooling flows. The spectra of 
PKS 1404-267 and A2029 have been extracted from the whole spatial extension of 
the object. The expected position of the Fl line is indicated with | . Flux is given 
in units of erg s~ 1 cm~ 2 A~ 1 . 
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4.2.1 Data analysis 

a) Extraction of the spectra 

For those objects for which no extended emission lines were detected, we com- 
pressed the whole spatial extension to extract a 1-dimensional spectrum. For 
those objects showing extended emission lines, we separated it into two different 
components: the nucleus, and the extended emission line spectrum. In order to do 
this, we subtracted from the spatial profile in the brightest spectral line, the spatial 
profile in the nearby continuum. In this way, we are left with the spatial profile 
of the line emitting regions. In some cases blobs appear well detached from the 
nuclear emission. But sometimes the extended and nuclear emission were mixed. 
In order to avoid contamination of the EELRs by the nuclear spectrum, we fit a 
Gaussian (PSF) to the nuclear region. The residuals after subtracting this Gaus- 
sian were considered to be due to the EELR. The resulting spectra are shown in 
Fig. 21, 22 and 23. 

It is possible that we may be venturing close enough to the nucleus to sample 
higher densities, more typical of the classical narrow line region (with densities 
higher than 10 3 cm -3 ). To verify this, when it was possible, we used the [SII] 
doublet to estimate the density. The values were always below 100 cm -3 , which 
indicates that the spectra are dominated by low density gas. In any event, even 
if the NLR dominated the emission, as a result of its high critical density the 
line Fl (in the dustfree case) was calculated to remain stronger with respect to 
other emission lines (like Ha or [Arlll]) than in the EELR. Our conclusions about 
inferring dust or not are therefore not affected by high densities. (Higher density 
models would imply a stronger calcium depletion). 

b) Line measurements 

IRAF routines were used to measure the emission line fluxes. For the blends, 
decomposition procedures were used, fitting several Gaussians at the expected po- 
sitions of the components. The fluxes were the values estimated for the Gaussians 
that better fitted the given profile. 

The upper limits on [Call] emission were measured directly from the noise 
level in the data. 

4.3 Observed and predicted Fl fluxes. 
4.3.1 Prediction of Fl fluxes 

The second component of the [Call] doublet (F2) is the weakest and is blended 
with the [OII]AA7320,7330 multiplet. The first component is the strongest and lies 
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well aside, some 30 A shortward of [Oil]. It is straightforward to isolate it given a 
reasonable spectral resolution. Therefore, for simplicity, we have based our study 
on the measurement of this line. The observed values are compared with the flux 
that photoionization models predict. MAPPINGS (Binette et al. 1993a,b) is the 
multipurpose photoionization-shock code that we used for the modeling. It con- 
siders the effects of dust mixed with the ionized gas: extinction of the ionizing 
continuum and of the emission lines, scattering by the dust, heating by dust pho- 
toionization and depletion of heavy elements. The appropriate input parameters 
for the models were justified in Chapter 3. 

The predicted value of the Fl flux has been computed in the following way: we 
used as a reference line (ref line), the strongest and/or the simplest to measure, 
like Ha, [SII]AA6716,6731, the [OII]AA7320,7330 blend and/or [ArIII]A7136. The 
appropriate photoionization model allows the prediction of R = re J\ ine , from which 
we can easily deduct the Fl flux: 

Predicted Fl flux = R x Measured flux in the reference line. 
In table 2, we show the reference lines used for each object and for each spatial 
region, the measured flux, the predicted R= re J\ ine ratio, the measured Fl flux 
and its theoretical value. We have used different lines for different objects and 
different zones. For the nuclear regions, we always base our predictions on the 
blend [Oil] AA7320,7330. The density in the nuclear narrow line region spands a 
wide range up to 10 6 cm~ 3 . The critical density of the [Call] forbidden lines is 
very high, and the flux of the lines under the nuclear conditions are predicted to 
be far stronger than the [Oil] lines (Ferland 1993). [Oil] can be de-excited due 
to high density so that if Call remains much smaller than [Oil], it is actually a 
more stringent test since high densities would have actually helped to increase the 
Call/Oil ratio. To measure the flux in the [Oil] blend, we assume that there is a 
negligible contribution from the F2 component. This assumption is reasonable: it 
is fainter that Fl and this line is near or under the detection limit in all cases. We 
can therefore assume that the [Oil] multiplet is not contaminated by F2. 

4.3.2 Comparison with the observations 

The comparison between the measured and predicted Fl flux values shown in 
table 2 demonstrates that, whenever the calculations were possible (sometimes 
there was no a reference line available, however), it is always fainter (non detected 
in most cases) than expected. This result is common to all the regions considered 
here: nuclear, EELRs and cooling flow filaments. It demonstrates that calcium is 
depleted and, therefore, the gas is mixed with dust, both in the nuclear as in the 
extended emitting gas in radio galaxies and cooling flow filaments. 

Donahue & Voit applied this same test to emission line nebulae in cluster 
cooling flows (1993). They interpreted the observed lack of [Call] emission as 
depletion of calcium onto dust grains in the ionized filaments. Our work support 
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Table 2: Emission line fluxes for the sample of galaxies observed (in units of 
erg s~ 1 cm~ 2 ). The superscripts M and p indicate measured and predicted values 
respectively. The last two columns show the predicted and measured values for 
the Fl line. For those objects where the calculations were possible, the observed 
line is clearly fainter than photoionization models predict. 



Name 


Spatial reg. 


ref line 


Flux M 


Fl F 

ref line 


Flux{Fl) p 


Flux{Fl) M 


NGC1052 


Total 


None 








< 6.99e-16 


NGC6215 


Total 


None 








< 3.34c-16 


NGC7552 


Total 


None 








< 3.95e-16 


NGC7714 


Total 


[OII]A7325 


1.59e-14 


>1.00 


> 1.59c-14 


< 6.84e-17 


PKS1404-267 


Nucleus 
5" South 


None 
None 








< 1.79e-16 

< 5.23c-17 


PKS2014-55 


Nucleus 


[OII]A7325 


1.54e-15 


>1.00 


>1.54c-15 


< 3.34c-17 




EELR (E) 


Ha 


1.47e-15 


- 0.45 


6.62e-16 


< 2.60c-17 




EELR (W) 


Ha 


1.06e-15 


- 0.45 


4.77e-16 


< 3.18e-19 


PKS2152-69 


Nucleus 


[OII]A7325 


7.48e-15 


>1.00 


>7.48e-15 


1.92e-16 




EELR 


[SILU6731 


1.33e-15 


- 0.40 


5.32e-16 


< 3.03c-17 




Cloud 


[ArIII]A7136 


3.86e-16 


- 4.50 


1.74e-15 


< 2.24c-17 


PKS2158-380 


Nucleus 


[OII]A7325 


4.03e-15 


>1.00 


>4.03c-15 


3.15e-16 




EELR (E) 


[ArIII]A7136 


7.28e-16 


- 5.0 


3.64e-15 


< 2.55c-17 




EELR (W) 


[ArIII]A7136 


3.00e-16 


- 5.0 


1.5e-15 


< 5.00c-17 


PKS2356-61 


Nucleus 


[OII]A7325 


9.39e-16 


>1.00 


>9.39e-16 


2.78e-16 




EELR 


Ha 


5.36e-16 


- 0.18 


9.6e-16 


< 1.03c-17 


PKS2300-18 


Nucleus 


[OII]A7325 


1.80e-15 


>1.00 


>1.80c-15 


4.37e-16 




EELR 


Ha 


3.61e-16 


~ 0.40 


1.44e-16 


< 1.97c-17 


2A 0335+096 


cD galaxy 


[SILU6725 


7.46e-16 


- 0.40 


2.98e-16 


< 2.72c-17 




EELR (E) 


[SII]A6725 


2.24e-15 


- 0.40 


8.96e-16 


< 1.29e-17 




EELR (W) 


[SII]A6725 


4.22e-16 


- 0.40 


1.68e-16 


< 1.40e-17 




companion 


[SII]A6725 


7.89e-15 


- 0.40 


3.16e-15 


< 1.69c-17 


A2029 


Total 


None 








< 1.43c-16 


A2597 


Nucleus 


[OII]A7325 


1.79e-15 


>1.00 


>1.79c-15 


5.14e-16 




EELR 


Ha 


2.32e-15 


- 0.50 


1.16e-15 


< 5.86c-17 
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these results, in a different sample of objects: the NLR in AGNs and the filaments 
in objects classified as 'cooling flows' are well mixed with dust. 

The most important result is the confirmation that dust exists mixed with the 
ISM in low z radio galaxies. 

4.4 Discussion 

4.4.1 Implications on the origin of the gas 

As explained in Chapter 3, there are some pieces of evidence indicative of dust 
in the extended ionized gas in active galaxies, mainly derived from polarization 
measurements that show the existence of scattered nuclear light over large spatial 
scales, although it is difficult to discriminate between dust and electrons as the 
scattering agent. Our results confirm that dust exists mixed with the ISM in low 
redshift radio galaxies. 

As mentioned before, this produces discrepancies with the traditional cooling 
flow theory, which would now be required to explain the formation of dust in a 
shorter time than the cooling time! According to this theory (see Fabian 1994 for 
a review), the galaxies, groups and clusters of galaxies were formed out of gas that 
collapsed gravitationally. During this process, gravitational energy was released 
which heated the clouds and a hot X-ray atmosphere would remain. Inhomo- 
geneities in the gas would cause matter from the hot atmosphere to cool down and 
fall towards the center of the galaxy or cluster. The resulting filaments would even- 
tually be visible at optical wavelengths and would emit strong lines. For massive 
galaxies and cluster, the cooling process would have been slower than for normal 
galaxies and the hot atmosphere would still exist, with typical temperatures of 
several million K and emitting strongly in the X-ray band. 

If this is actually the origin of the gas in the EELRs (and cooling flow fila- 
ments), it should clearly be devoid of dust. Any dust introduced in the hot intra- 
cluster medium would be sputtered and rapidly destroyed (time scale of the order 
of 10 7 years) (Draine & Salpeter 1979), much before the filaments cooled down 
and became visible. Is there a way to introduce dust in the filaments during the 
cooling process? The common place where dust is formed is stellar atmospheres. 
Do stars exist in the accreted gas? Such gas must be deposited in some form and 
it has been generally thought that the more plausible fate for most of the gas is 
the formation of new stars (e.g., O'Connell & McNamara 1989; Fabian, Nulsen & 
Canizares 1991). Although there are some indications of star formation in cooling 
flow galaxies, the latter seems to be taking place in the inner parts (over the central 
few kpc) (Cardiel, Gorgas & Aragon Salamanca 1995, Fabian 1994). There is a 
lack of evidence of star formation (or any kind of stars) at large distances, out of 
the main body of the galaxies, where the EELRs filaments are still visible. 
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In high z RG there is an extended blue continuum aligned with the radio axis 
that has often been interpreted as due to young stars. However, as I mentioned 
already before, further results (e.g., Tadhunter, Fosbury & di Serego 1988; Cimatti 
et al. 1993) proved that this continuum is polarized and the contribution of a 
hidden quasar continuum scattered by dust and/or electrons in the ISM can well 
be an important component (in fact, this extended polarized UV continuum has 
been used as a proof of the existence of dust in the EELRs of high z radio galaxies). 
This effect has also been detected in some radio galaxies at low z (Cimatti & di 
Serego Alighieri S. 1995). Altough young stars may exist, we don't have a clear 
idea about their overall contribution. 

There is another mechanism which can form dust: if most of the cooled gas 
from a flow does not form stars with normal IMF, maybe it remains as cold clouds 
or as low-mass stars. There are evidences of X-ray absorption (White et al. 1991; 
Mushotzky 1992; Allen et al. 1993) in cooling flows in cluster of galaxies. The 
absorbing material could be in the form of cold gas embedded in the hot intra- 
cluster medium of the cooling flow (White et al. 1991; Daines, Fabian & Thomas 
1994). This cold gas, very slightly photoionized by X-rays from the surrounding 
hot corona, can become molecular (Ferland, Fabian & Johnstone, 1994). Fabian, 
Nulsen & Canizares (1994) propose that the conditions suitable for dust to form 
may occur in this cool gas, through the condensation of gaseous particles. The au- 
thors propose that the gas which is cooling towards the center could be a mixture 
of both the molecular gas and the very hot gas and, therefore, could contain the 
existing dust. However, even if such a scheme was possible, the dust grains will not 
remove the calcium which already existed in the hot gas; the temperature is too 
high for the condensation of calcium onto the dust grains. Therefore, we should 
observe the Call lines from the hot gas when it cools down, even if it contains 
dust. 

On the other hand, the interaction scenario (gas accreted from outside the 
galaxy, as a result of recent tidal interactions or mergers between two or more 
components) predicts the existence of dust mixed with the gas, the one already 
existing in the interacting objects. Some morphological evidences and theoretical 
ideas (see Chapter 3) support this scene. Heckman et al. (1986) showed that a 
large fraction of powerful radio galaxies have morphological features (shells, tails, 
loops, etc) similar to those produced in numerical simulations of galaxy interactions 
(e.g., Toomre and Toomre 1972, Quinn 1984). Kinematic measurements show that 
the radio galaxy EELR generally have a high specific angular momentum which is 
difficult to reconcile with the cooling flow picture (Tadhunter, Fosbury & Quinn 
1989). 

Therefore, if there is a common origin for the EELRs of all radio galaxies, our 
results suggest that it is mergers or tidal interactions. However, we don't exclude 
the possibility of an origin which is not universal, that is, which may differ from 
object to object. 



Internal dust in Large Scale Nebulosities in RGs (II) 



80 



4.5 Conclusions. 

We have confirmed that the gas in extended emission line regions in radio galaxies 
at low z is mixed with dust. 

Our results support the existence of dust mixed with the gas in the Narrow 
Line Region and in the cooling flow filaments. 

If there is an universal origin for the EELRs, the existence of internal dust 
favours mergers or tidal interactions as the most plausible scenario. 
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Chapter 5 



Effects of dust and resonant 
scattering on the UV spectrum of 
radio galaxies. 

Villar-Martfn, Binette & Fosbury 1995, A&A, submitted 

Abstract 

In the powerful, high redshift radio galaxies, it is believed that the dominant 
source of ionization for the gas is the hard radiation field associated with the ac- 
tive nucleus. The photon source is generally external to the clouds being ionized 
and so the geometrical perspective from which the gas is observed and the pres- 
ence and distribution of dust must be properly accounted for in the diagnostic 
process. In this paper, we examine the formation of three strong lines, CIVA1549, 
Lja and CIII]A1909 which are often observed in the nuclear and extended emis- 
sion from these sources. We find that the observed trends, in particular the high 
CIVA1549/Lya ratio, are often better explained by geometrical (viewing angle) 
effects than by the presence of large quantities of dust either within or outside 
the excited clouds. We show that condensations of neutral gas along the line- 
of-sight can increase the observed CIV/Lya ratio, by reflecting photons near the 
wavelength of Lya,. The existence of HI absorption clouds (i.e., mirrors) external 
to the emission region leads also to the presence of large, diffuse haloes of what 
appears to be pure, narrow Lya emission. 
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5.1 Introduction 

The strong, spatially extended, rest-frame ultraviolet emission lines observed in 
high redshift radio galaxies provide one of the principal diagnostics in establishing 
the state of the interstellar medium in galaxies at early epochs. The presence 
of a blue continuum and emission lines from regions aligned with the radio axis 
(McCarthy et al. 1987; Chambers, Miley & van Breugel 1987) warned us that much 
of the observed optical radiation might be associated with the nuclear activity and 
so may not be giving us a clear picture of the stellar processes which are of great 
interest in studies of galaxy formation and evolution. Subsequent work has shown, 
indeed, that much of the blue light is scattered, polarized nuclear radiation (e.g., 
Tadhunter et al. 1989; di Serego Alighieri et al. 1993; Cimatti et al. 1993) 
and that the emission lines have a high ionization state and cannot result from 
photoionization by normal stars (McCarthy 1993). It is clearly necessary, therefore, 
to reach a clear understanding of the physical processes involved in the formation 
of the various lines and continua to be able to disentangle the stellar and the 
AGN-related sources. 

For objects at highz, the UV rest-frame lines are shifted into the optical 
band and the spectrum is generally dominated by Lya, CIVA1549, HeIIA1640 and 
CIII]A1909. The strength of the high ionization lines suggests the presence of a 
hard photoionizing continuum which could originate at the AGN itself (Robinson 
et al. 1987) or be associated with fast shocks generated in extranuclear regions 
by the radio jets (Sutherland, Bicknell & Dopita 1993). The strong radio/optical 
asymmetries observed in these objects which exhibit the 'alignment effect' (Mc- 
Carthy, van Breugel & Kapahi 1991) may simply result from a one-sidedness in the 
distribution of material. It it is clear, however, that correlated line and continuum 
asymmetries could be produced by dust scattering and line fluorescence for sources 
where the radio axis falls significantly away from the plane of the sky. 

In this work, we concentrate on modeling the high excitation lines for which 
rather extreme ratios relative to Lya have recently been reported. The presence 
of dust has been universally invoked to explain the weakness of Lya which is 
a resonance line and therefore, due to multiple scattering, more susceptible to 
absorption. We explore the fact that any resonance line will be extremely sensitive 
to geometrical factors, an aspect of the problem which has so far been overlooked in 
modeling the UV lines. If in radio-galaxies the distant gas clouds are photoionized 
from the outside by partially collimated UV radiation emitted by the nucleus, the 
line formation process — particularly for the resonance lines — is very different 
from that of the internally ionized HII regions. The escape of resonance line 
photons is strongly influenced by the presence of voids between the line emitting 
clouds. 

We have collected from the literature the observed line ratios for a number 
of radio-galaxies at high z (z >1.5) in which no contribution from any nuclear 
BLR is apparent. We have built a diagnostic diagram consisting of the lines 
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CiVA1549/Lya vs. CIVA1549/CIII]A1909, in which we compare the position of 
the objects with photoionization models which not only consider the effects of 
internal dust but also those of the viewing perspective — the angle between the 
incoming ionizing radiation and the observer's line of sight. Our concentration on 
the particular class of radio galaxies is purely for pragmatic reasons. It is these 
objects, which we presume to harbour a powerful quasar which is hidden at opti- 
cal/ultraviolet wavelengths to our line of sight, which are most readily found and 
studied at high redshifts where we have access to the ultraviolet spectrum from 
groundbased observations. Our conclusions should be equally applicable to other 
classes of AGN. 

For some objects, Lja is observed to be fainter with respect to CIV than pre- 
dicted by dust-free photoionization models. The explanation previously proposed 
to explain the weakness of Lja with respect to Ha or H/3 has been dust destruction 
of resonant Lja photons. This is not borne out by our calculations in which we 
have used arbitrary amounts of dust and found that this cannot simultaneously 
weaken Lja while leaving the CIV/CIII] ratio relatively unchanged since resonant 
CIV suffers also from dust absorption. Alternatively, by varying the proportions 
of the illuminated and the shadowed cloud faces which contribute to the observed 
spectrum, we are better able to match the data. 

As we find that geometry alone (with or without internal dust) can in principle 
explain most of the specific line ratios observed (fainter Lja compared with either 
CIV or Hell), we also discuss the possibility of a patchy outer halo of neutral gas 
to account for the diffuse Lya seen in some cases to extend much beyond the CIV 
emitting region and even the outermost radio lobes. Reflection by cold gas of the 
brighter Lja emitting side of the ionized clouds would lead to a narrower profile for 
such a diffuse component. Another possibility is that part of the beamed nuclear 
continuum and BLR radiation might be reflected at the wavelength of Lya by thin 
matter-bounded photoionized gas at very large distances from the nucleus leading 
to a diffuse Lja component aligned with the radio axis. It appears to us that 
geometrical perspective effects are an essential component of the interpretation of 
the UV spectrum of radio-galaxies whether or not dust is present. Furthermore, a 
spectrum in which only Lja appears does not necessarily imply starburst activity, 
other lines must be observed before the existence of HII regions can be inferred. 
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5.2 Data sample and modeling procedure 

5.2.1 The data 

We have constructed a data sample containing galaxies at high z for which the 
CIV, Lya and, in most cases CIII], emission lines have been measured. Since 
very high densities such as those encountered in the BLR alter significantly the 
line formation and transfer processes, we have excluded those objects which show 
evidence of a BLR. In Table 3, we list the object names, the line ratios of interest 
to us here, the redshift and the reference to the observations. The larger fraction of 
the data are taken from the recent thesis by van Ojik (1995) which includes objects 
selected on the basis of a very steep radio spectrum. Probably by virtue of the radio 
selection, these sources populate the region of the line ratio diagram (Fig. 29) with 
lower CIV/CHI] ratios (lower ionization parameter) than the previously published 
objects. The line measurements refer to the integrated emission from the object 
collected with a long slit aligned with the major axis. 

5.2.2 The model and its parameters 

The data are compared to photoionization models computed using the code MAP- 
PINGS. The version described in Chapter 2 is particularly suited to this study since 
it considers both the effects of the observer's position with respect to the emitting 
slab and the ionizing source (see Fig. 4, Chapter 2), distinguishing between the 
spectrum seen from the back and from the front of the slab. Phenomena related 
to dust are also properly considered. The treatment of the escape of resonant CIV 
and Lja photons in a dusty medium (a fundamental proces that needs to be well 
understood) is described in Appendix B of Binette et al. (1993a) and is based on 
the results of Hummer & Kunasz (1980). The presence of dust implies depletion 
of refractory trace elements and this effect has also been considered, as described 
in Apendix A. The calculations consider the gas pressure to be constant (isobaric 
models) and so the density behaviour with depth in the cloud is determined by 
the behaviour of the temperature and the ionization fraction of the gas. 

5.2.3 Adopted physical conditions 

The detailed studies of the optical emission lines of low z radio galaxies provide our 
basic reference for the properties of the emitting gas and the ionizing continuum 
(e.g., Robinson et al. 1987, hereafter RBFT87). As a guide in choosing the input 
parameters, we have assumed that the excited gas of very high z radio galaxies 
(z > 2.5) has similar properties to that of the low redshift (z < 0.1) objects. As 
explained in Chapter 3, for the extended (EELR) as well as the (narrow) nuclear 
emission lines, photoionization by a hard continuum appears to best explain the 
various line ratio diagnostic diagrams: a power law of index a ~ — 1.4 (f v oc u +a ), 
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Table 3: Observed UV line ratios for several high z RG with not apparent broad 
component 



Name 


CIV/CIII] 


CIV/Lya 


z 


ref. 


Average RG 


2.054 


0.118 


0.1< z <3 


McCarthy (1993) 


MG1019+0535A 


2.12 


1.24 


2.76 


Dey et al. 1995 


F10214+4724 


3.68 


8.75 


2.29 


Elston et al. 1994 


TX0211-122 


3.33 


0.91 


2.34 


van Ojik et al. 1994 


3C294 


0.83 


0.10 


1.79 


McCarthy et al. 1990a 


0902+34 


— 


0.11 


3.40 


Lilly 1988 


3C256-3C239 


1.90 


0.14 


1.82 & 1.78 


Spinrad et al. 1985 


0200+015 


1.05 


0.241 


2.23 


van Ojik 1995 


0214+183 


1.67 


— 


2.13 




0355-037 


1.17 


0.24 


2.15 




0417-181 


— 


0.40 


2.73 




0448+091 


0.44 


0.10 


2.04 




0529-549 


0.22 


0.05 


2.58 




0748+134 


1.29 


0.29 


2.42 




0828+193 


0.95 


0.14 


2.57 




0857+036 




0.39 


2.81 




0943-242 


1.70 


0.19 


2.92 




1138-262 


0.62 


0.06 


2.16 




1357+007 




0.19 


2.67 




1410-001 


1.58 


0.16 


2.36 




1545-234 


0.61 


0.17 


2.76 




1558-003 


2.25 


0.18 


2.53 




2202+128 




0.25 


2.71 




2251-089 


2.20 




1.99 




4C23.56 


3.40 




2.48 




4C24.28 




0.23 


2.88 




4C26.38 


3.71 




2.61 




4C28.58 


0.17 




2.89 




4C40.36 


1.05 




2.27 




4C41.17 




0.05 


3.80 




4C48.48 


2.18 




2.34 




4C60.07 




0.27 


3.79 
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that has produced remarkably good agreement with the optical line ratios at low 
z, will be adopted as ionizing energy distribution. 

Concerning metallicity, RBFT87 indicated that the abundances cannot be 
much higher than solar values. They could however be lower by a factor of a few, 
a distinct possibility in the case of the extranuclear gas in very high z galaxies. 
Unless specified otherwise, we have adopted Z = 1 in our calculations and have 
verified that lower abundances do not in any way affect the conclusions reached 
here. 

The red [SII] doublet ratio in low z radio galaxies indicates densities for the 
extended emission line regions which are lower than 100 cm -3 . We therefore adopt 
the low density regime, specifically n H = 100cm" 3 , since the extranuclear gas 
generally dominates the line luminosities. 

5.3 Models and comparison with observed UV 
lines 

In this section we first extend the assumption of photoionization to modeling the 
line emission in the ultraviolet rest-frame of high z radio galaxies. We distinguish 
the effects of having a back and a front view of an externally photoionized slab 
(section 5.3.2). In 5.3.3, we introduce internal dust and discuss how, alone, it is 
insufficient to explain the low value of the Lya/CIV ratio in some objects. In 5.3.4 
we distinguish between the resonance line 'mirror' intrinsic to our photoionized 
slab and the possibility of having external cold gas — which is shadowed from the 
ionizing source — acting as a reflector of Lja photons. 

5.3.1 The photoionization assumption at low and high z 

At very high redshift, observational access to the normal optical plasma diagnostic 
lines is restricted and the available data set is small. Also, where measurements are 
available, few lines are measured in any given source. It is nevertheless interesting 
to compare the [OIII]A5007/H/3 value observed at very high z to that of low z radio 
galaxies as plotted in Figure 24. The curved line corresponds to the sequence of 
photoionization models of RBFT87 which reproduce reasonably well the optical 
(rest frame) line ratios of the extended emission line regions and their associated 
narrow nuclear line emission. The range in U is [10 -4 , 10 -1 ]. The three lower 
horizontal lines correspond to high z radio galaxies while the upper one corresponds 
to the "average" radiogalaxy spectrum as defined by McCarthy (1993). The high z 
objects are not substantially different from the low z sample although the apparent 
trend towards weaker [OIII]/H/3, if confirmed, might indicate substantially lower 
metallicities: low enough to overcome the higher ratios produced by the higher 
kinetic temperatures in moderately underabundant objects. Such an effect could 
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also be produced by insufficient spatial resolution at high z to separate the high 
from the low excitation regions. 
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Figure 24: Excitation diagnostic for the ionized gas in low z radio-galaxies 
([OIII]A5007/H/3 vs. [OI]A6300/[OIII]A5007). Open triangles correspond to extended 
emission line regions (EELR) while solid triangles correspond to nuclear regions. The 
three lower horizontal lines indicate the observed [OIII]A5007/H/3 ([OI] is not measured) 
ratio of very high z radiogalaxies (Eales & Rawlings 1993) while the long-dash line cor- 
responds to the "average" radiogalaxy spectrum as derived by McCarthy (1993). The 
dotted line correspond to the sequence of photoionization models of RBFT87. log U is 
in the range [-4,-1]. 



5.3.2 Effects of viewing direction on the UV lines 

To see if these models reproduce the observed CIV/Lya and CIV/CIII] ratios, we 
have presented the observed and predicted values in Fig. 25. The line ratios shown 
correspond to calculations in which the ionization parameter is varied, generating 
a sequence in U . The three sequences differ in density or metallicity and logC/ is 
in the range [-4,-1]. 

We find that somewhat higher values of U than are tipically used in the optical 
are required to reproduce the CIV/CIII] ratio observed. We have investigated 
possible effects due to variations in metallicity and density (see Fig. 25), but none 
of them is able to explain the high values of the line ratios, unles a higher ionization 
parameter than in the optical is considered. This is most probably a consequence 
of the different way in which the objects at high and low z are selected: the distant 
sources are all very powerful radio sources with luminous AGN. For a power law 
of index -1.4, the optimum value is U ~ 0.1. Hereafter, diagrams will only cover 
the range of [0.01,0.1] in U. 
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Figure 25: Observed and predicted UV line ratios. Filled triangles are the ratios of 
objects observed by several authors. Open circles correspond to data taken from van 
Ojik's thesis, selected on the basis of a very steep radio spectrum. The open diamond 
is the average radiogalaxy spectrum of McCarthy (1993). Dotted lines represent models 
in which the front and back spectra have been summed, log U is in the range [-4,-1], 
as in Fig. 24. It is clear that high values of the ionization parameter U are required to 
reproduce the high CIV/CIII] value observed. 

A notable observation is the distribution of observed points above the model 
loci in Fig. 25 even at the highest value of U. In the extreme objects at least — like 
F 10214+4724 and TX 0211-122 — this results from a weakening of Lya rather 
than from atypical values of CIV/CIII]. It is usually claimed that the destruction 
of Lya photons by resonance scattering in the presence of dust is the explanation 
for its faintness, but why does not the same process reduce CIV which is also a 
resonance line? Is there an alternative explanation for this selective dimming of 
Lya? 

To answer this question, we examine the geometrical aspects of the line forma- 
tion process. Each emitting cloud is approximated as a plane parallel slab which 
contains a fully ionized region and a partially ionized zone where low ionization 
species co-exist (e.g., O , S + , etc) with a mixture of H° and H + . In principle there 
can be an additional neutral zone which does not contribute to the emission line 
intensities (see Fig. 4, Chapter 2). 

In this section we consider the dust-free case. For most lines (like CIII], 
Hell, H/3, [OIII], etc) line opacity is negligible and the line is emitted isotropically 
with photons escaping freely in all directions. However, when the line opacity is 
important as it is for CIV and Lya, line scattering occurs which increases the path 
length. Another important effect of large optical depths is that the line photon 
will not escape isotropically. A resonance line photon following many scatterings 
must statistically escape in the direction of highest escape probability which can 
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Figure 26: Influence of viewing geometry on the UV line ratios. The dotted line corre- 
sponds to the line spectrum seen from the back of the slab (cf, Fig 4, Chapter 2) while the 
dashed line corresponds to the spectrum seen from the front (UV irradiated face). It is 
apparent that perspective plays a very important role on the CIV/Lya ratio. The solid 
line represents models obtained by summing the back and front spectra which would 
represent symmetric case where equal numbers of clouds are observed with shadowed 
and illuminated faces. The data are the same as in Fig. 25. 

be shown to be the front for the photoionized slab depicted in Fig. 4 (Chapter 
2). In the case of Lya, the reason is that — while Lya photons are generated 
more or less uniformly within the slab (except within the PIZ) by recombination 
- the neutral fraction and therefore the incremental line opacity oItl/oIx increases 
monotonically as a function of depth as discussed in more detail by Binette et 
al. 1993b (BWVM3). This means that for an open geometry like that shown in 
Fig. 4 (Chapter 2), the zone of equal escape probability of front vs back occurs far 
beyond the point where half the luminosity of Lyct is produced. For the collisionally 
excited CIV line, the tendency to escape from the front also exists although it is 
less pronounced. It arises mainly because the emissivity of CIV is larger towards 
the front due to the temperature gradient across the C +3 zone. Note that while 
CIV is emitted and scattered within the rather limited zone containing C +3 , Lja 
remains subject to scattering outside the region where it is produced. While most 
Lja emission is produced in the fully ionized zone, most of the line opacity occurs 
within the PIZ. The presence of a layer of neutral gas beyond the PIZ will increase 
the anisotropy of Lya escape. 

The effects described qualitatively above are shown in Fig. 26 using detailed 
photoionization calculations. We present the same sequence of dust-free models as 
in Fig. 25 but distinguish between the spectrum seen from the back — equivalent to 
observing the clouds through the PIZ — from that seen from the front — equivalent 



Resonance scattering and dust in HZRGs 



92 



to seeing the UV irradiated side. The differences are striking: both Lya and to a 
lesser extent CIV are fainter when seen from the back. The CIII] line is isotropic 
in the dust-free case. The fact that Lya is more affected by perspective is due 
to the significant amount of neutral hydrogen (i.e., large line opacity) in the PIZ 
which acts as a mirror. Although we have considered a very simplified geometry 
in our calculations, the method nevertheless treats properly the essential physical 
effects and indicates how important the viewing direction is in this open geometry. 

Fig. 26 suggests that perspective effects alone (without any dust) are sufficient 
to explain the weak Lya seen in some objects. Ionization bounded calculations with 
U — 0.1 imply total hydrogen column densities (H + region + PIZ) N H ~ 10 22 cm -2 
(of which about 60% is ionized). Adding a modest neutral zone beyond the PIZ of 
~ 4 10 21 cm" 2 would double the CIV/Lya ratio without affecting in any way the 
CIV/CIII] ratio. It seems, therefore, that a geometry where we see preferentially 
the ionized gas from the side of the PIZ gives us an explanation for the weakness 
of Lya. 

How would this apply to the EELR of powerful radio galaxies? In a very 
simplified scheme, we can imagine (see Fig. 27) that the clouds seen from the 
nearside cone are seen from a direction which we approximate as the back perspec- 
tive in our slab calculations while clouds on the far side would be seen from the 
front. The studies of McCarthy, van Breugel & Kapahi (1991) which emphasized 
the one-sideness of the line brightness distribution suggest that the observer with 
limited spatial resolution at very high z will be biased towards either a back or 
front dominated perspective depending on whether it is the near or the farside 
illuminated cone which is intrinsically brighter. Our proposed interpretation of 
the high CIV/Lya ratio of some objects in Fig. 26 is that they correspond to the 
case where the brightest clouds are seen from behind. Note that, because of the 
intensity weighting, there is rather little difference between the pure front perspec- 
tive and the sum of equal numbers of front- and back-clouds, a small effect which 
is not easily distinguished from that of slightly reducing the value of U . 

Our conclusion from this section is that perspective effects in an open - 
externally illuminated cloud — geometry go a long way towards explaining the 
behaviour of the CIV, Lya, CIII] line ratio diagram which implies that, when 
spatially resolved spectra become available, marked asymmetries in the CIV/Lya 
ratio could arise when the axis of the illuminated cones forms a large angle with 
the sky plane. 

5.3.3 The effects of internal dust 

There is evidence for the existence of dust in some high z radio galaxies. The 
detection of 4C41.17 (z = 3.8) and B2 0902+34 (z = 3.5) (Chini & Kriigel 1994; 
Dunlop et al. 1994) and 8C1435+635 (z = 4.26) (Ivison 1995) in the mm spectral 
range is attributed to warm dust. Also, the IRAS galaxy F10214+4724 (z = 
2.29) has been shown, from the far infrared flux, to contain ~ 1O 8 M of dust, 
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Figure 27: Front and back clouds within the ionized cones of a radio galaxy. 

although this figure may be reduced by the gravitational lensing amplification 
factor (e.g. Eisenhardt et al. 1995; Serjeant et al. 1995). In addition, the mid- to 
far-infrared measurements at intermediate z are explained as emission from warm 
dust (Heckman, Chambers & Postman 1992). We should emphasize, however, 
that we have no direct measurement of how this dust is spatially distributed. In 
the event that this infrared emission arises from the reprocessing of higher energy 
photons from the AGN by a dusty torus, it has no direct bearing on our modeling of 
ionized gas at tens of kpc. However, there is considerable evidence that aligned blue 
polarized continuum is the result of scattering of the anisotropic nuclear radiation 
field by dust (e.g., Cimatti et al. 1993). In this case it is very likely that at least 
some of the dust is internal to the extended line emission regions. In this section 
we consider the effects of including dust within the gas clouds which emit the UV 
lines. 

We illustrate the effects of dust mixed with the emitting gas in Fig. 28 where 
we plot sequences of models which correspond to the front perspective for three 
different dust-to-gas ratios: /i = 0, 0.3 and 1. The effects on the line ratios are 
evident. The resonance scattering suffered by CIV and Lya increases their path- 
lengths many times and, therefore, the probability of their being absorbed by dust 
grains is much higher than for CIII]. In the case of Lyct, however, the geometrical 
thickness of the H + region exceeds greatly that of the C +3 since we observe more 
than one stage of ionization of metals (e.g. C +2 ). Any reasonable parameters for 
the ionization structure of a photoionized slab with Z ~ 1 indicates that the opac- 
ity in Lya greatly exceeds that of CIV, which implies a larger pathlength increase 
for Lja than for CIV. This results in relatively more Lya absorption by dust. 
This effect explains how the ratio CIV/Lya increases somewhat with increasing fi. 
Dust absorption of resonant CIV on the other hand causes a comparable decrease 
in CIV/CIII]. What is important in these results is that, even with a concentra- 
tion of internal dust as high as fi = 1 (equivalent to that in solar neighbourhood 
cold clouds), it is not possible to reproduce the high ratio CIV/Lya > 1 which is 
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Figure 28: Effects of varying the amount of internal dust as seen from the front per- 
spective. As in the previous figures, the various lines correspond to sequences in U. 
Short-dashed line corresponds to models with jj, = 1.0, the long-dashed to models with 
H = 0.3 and solid line to dust free models. 
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Figure 29: Effects of perspective combined with small quantities of internal dust. The 
dust-free U sequences of Fig. 26 are repeated in this figure. The dotted line connected 
by open circles corresponds to the last model with U = 0.1 (but with metallicity Z = 0.3 
in which the dust content is increased in proportion up to \i = 0.017. Larger quantities 
of dust do not improve these results because dust severely extinguishes both the CIII] 
and CIV lines which are produced in the front layers. 
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observed in some objects and has been attributed to dust. 

When the clouds are viewed from behind, even with \x = 0.3 they are suffi- 
ciently opaque that all the UV lines are severely absorbed. With Ty = 5 f 0~ 22 [iNh 
~ 5 and \x = 1, all of the high excitation UV lines become absorbed within the PIZ. 
To produce an acceptable spectrum, without reddened CIV/CHI] and CIV/Hell 
ratios, as seen from the back of an ionization bounded slab with U = 0.1, we 
have to use much smaller amounts of dust like /i ~ 0.017 (2% of local ISM) in 
models. In Fig. 29, we plot the back and front sequences for such models. These 
can reproduce the weak Lya objects although this result is obtained only for the 
back spectrum, emphasizing that perspective is the dominant factor. It should be 
noted that the amount of extinction within the ionization bounded slab implied 
by /j, — 0.017 (i.e., Ay ~ 0.1-0.2) is consistent with the quantity of small dust 
grains needed to explain the extremely blue continuum of the "detached" ionized 
cloud in PKS 2152-69 (di Serego Alighieri et al. 1988; Magris & Binette private 
communication), supposing that the continuum energy distribution is the result of 
dust scattering of the nuclear radiation. 

Note that because the scattering/ absorbing dust is locally kinematically linked 
to the emission line gas, we require much smaller column densities of HI than the 
absorbing screen proposed by van Ojik et al. (1994) (~ 10 23 cm~ 2 ) to explain 
CIV/Lya ~ 1. The line widths and centroids of the PIZ and of the fully ionized 
gas are expected to be quite similar within each of the emitting clouds, the ensemble 
of which could have a greater 'turbulent' velocity dispersion. 

To reproduce the extreme case of the IRAS galaxy F10214+4724 with CIV/Lya 
~ 10, we need additional neutral gas beyond the PIZ. For instance, it requires only 
a column density of N H o ~ 1.5 10 21 cm -2 assuming \x = 0.017. Alternatively we 
might consider that u, within the PIZ increases with depth as the degree of ion- 
ization decreases. In this case, no additional HI gas would be required. If this 
IRAS galaxy is indeed an extreme Seyfert 2 as argued by Elston et al. (1994), 
then a closed, dust enshrouded geometry as proposed by BWVM3 to explain the 
Lyman and Balmer decrements in Seyfert 2 would be more appropriate than the 
open geometry adopted here which may be applicable only to the truly extended 
large scale gas. It is likely that objects with unusually strong NVA1240 emission 
(such as in F10214+4724 and TX021 1-122) are cases where the NV originates pre- 
dominantly from the inner NLR. A high NV/CIV ratio indicates very enriched gas 
which is not unexpected within the inner parts of an AGN (Hamman & Ferland 
1993) and is consistent with the lack of convincing evidence that NV is spatially 
resolved in any of these objects. 

5.3.4 Neutral gas mirrors 

So far we have considered the effects of scattering by gas which forms part of 
the line emitting clouds: we refer to this as an 'intrinsic' process. Similar effects 
could be produced by a large-scale distribution of predominantly neutral material 
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surrounding the emitting regions: we will refer to this as the 'extrinsic' case. Any 
extrinsic neutral gas component with a non-negligible covering factor could affect 
the observed spectra in a way which mimics that of the back perspective described 
earlier. Let us suppose that this outer material is broken up into cold gas clumps 
which are randomly distributed. In such a case, some Lja photons which leave the 
ionized cones will escape the region through the holes between the external clumps 
while others will strike the neutral clumps and be immediately scattered away to 
escape eventually through another hole in a different direction (see Fig. 30). If 
observed with sufficient spatial resolution, such a geometry would result in holes 
in the Lja brightness distribution due to reflection by intervening clumps as well 
as diffuse sources corresponding to reflection from clumps on the far side of the 
source. The bulk of the Lja luminosity would be preserved but redistributed 
on an apparently larger scale than the true line emitting clouds. Only a more 
closed geometry would result in a significant destruction of Lja, assuming that 
the interclump space does not contain pure dust segregated from the gas phase. 

The reflection efficiency of clouds will, of course, depend strongly on the rela- 
tive velocity fields of the emitting and the cold regions. If the extended gas has a 
large scale ordered motion but small 'microturbulence', the reflection effects would 
be localized. Broad Lja emission and continuum from the AGN could, however, 
be scattered by any of the extranuclear clouds (see below, 3C294). In general, we 
would expect the diffuse, scattered Lja to show a narrower line profile than the 
integrated emission profile. 

We will review the observational evidence for the existence of such large scale 
mirrors. 

The radio galaxy PKS2104-242 

PKS2104-242 (McCarthy et al. 1990b) shows very extended emission lines of Lja, 
CIV and Hell associated with continuum knots which lie between the radio lobes 
and are aligned with the radio axis. The Lja image shows emission resolved into 
three distinct clumps, two of them corresponding roughly with the two continuum 
knots. In addition, there appears to be a low surface brightness halo of emission 
surrounding the entire object. Is this halo a consequence of reflection by neutral 
material of either Lja emission or of the intense nuclear continuum, or are these 
photons emitted locally by H + recombination? A way of discriminating between 
these two possibilities is to look for the detection of any other emission line in 
the halo. If Lja is the result of reflection by cold HI, there will be no other lines 
(except possibly resonant Mgll). If it is instead produced by recombination, other 
lines should be detected in the rest-frame optical band like Ha, [OIII]A5007, etc. 
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Figure 30: The extrinsic case. Neutral material external to the ionized regions can 
strongly influence the appearance of the object in Lya. In this figure, the neutral clumps 
(black clouds) cover a significant fraction of the ionized regions (grey+black clouds). 
The solid arrows represent Lya photons emitted directly in the ionized regions. Dashed 
arrows represent Lya photons that after striking neutral clumps are then reflected in 
other directions. The reflection by the neutral H atoms is so effective that the dust has 
little chance to interact with the photons before they find a hole to escape. 
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The radiogalaxy 3C 294 

From images of this object, McCarthy et al. (1990a) report the existence of Lya 
emission extending over a region covering 170kpc. The Lya emission is elongated 
and well aligned with the inner radio-source. An intriguing aspect of the observa- 
tions is the one-sidedness of the CIV emission (obtained with an long slit aligned 
with the axis of elongation). The decreasing linewidth of Lya on the side where 
CIV is absent (South) suggests to us the possibility that Lya in the south corre- 
sponds to scattered continuum and BLR Lya photons by HI clumps lying along 
this direction provided the ionizing radiation has already been reprocessed (filtered 
out) within the nuclear regions into NLR or BLR emission. 

Failing that (scattering by HI clumps), an alternative possibility resides in 
very thin matter-bounded photoionized sheets of gas which can be very efficient at 
reflecting the intense (beamed) nuclear continuum as well as BLR Lya photons. 
Due to its large scattering cross section, very small column depths of H° (a trace 
specie within the ionized phase) are sufficient to scatter effectively the impinging 
flux within the core of the Lya thermal profile. If we suppose the existence within 
the cone of ionization of the radio-galaxy of a population of clouds of similar 
physical conditions to that thought to apply to Lya forest clouds, the energy 
reflected due to resonant scattering by HI typically exceeds that generated within 
the clouds by reprocessing of the ionizing radiation (i.e. by recombination). In 
effect, adopting N H o = 10 13,8 cm~ 2 as a typical column density of such cloud, we 
derive an equivalent width in absorption (i.e. the scattered intensity) of 

EWf y s a « 0.15^111(4.2 10- u N H o) A = 0.15A 

which corresponds to the saturated part of the Lya curve of growth. Adopting 
for definiteness U = 0.1 and a = —1.4, we obtain that the fraction of reprocessed 
ionizing photons, F MB , within a very thin photoionized sheetQ is given by 

F MB w 1.8 10- 18 A^o = 1.1 10~ 4 
The equivalent width of Lya in emission as a result of recombination is simply 

EW[™ a « 390 F MB C f A = 0.04 (7/ A 

where Cf is the covering factor. If we integrate the absorbed intensity over the 
same area as that given by Cf, the ratio of scattered versus emitted Lya is ~ 4 
(the ratio between the above two equivalent widths). Allowing for the fact that 
the cloud will be exposed not only to the nuclear continuum but to the BLR Lya 
which typically peaks at twice the continuum level, the ratio can be expected to 
increase up to 8 times the Lya emitted by recombination. To account for the 

1 With these parameters, the total column density is given by Nh = N H + w 2.3 10 4 N H o. 
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observed FWHM < 1000 kms 1 in 3C294 (South), one most rely on a turbulent 
velocity field for the the Lya clouds. 

We ought to consider seriously the possibility that Lya on the southern knot 
of 3C294 corresponds to large scale scattered light by either HI or even ionized 
gas unless other bona fide emission lines were detected. So far, no lines other than 
Lya are observed and this is consistent with our suggestion. The only certain way 
of discriminating between reflection and in situ emission would be the detection at 
this radius of any non- resonance line in either the optical ([Oil], [OIII], H/5 etc.) or 
the UV (CII], CIII]). To conclude, we believe that HII regions or starbursts should 
not be considered the default emission mechanism in cases where only Lya is 
present, especially along the axis where the AGN supposedly collimates its intense 
nuclear continuum + Lya(BLR) light. 

Lya absorption in 0943-242 

The radio galaxy 0943-242 shows a Lya profile with several absorption features 
which cover the whole spatial extent of the Lya emission (Rottgering et al. 1995). 
The strongest of the absorbers is at least as extended spatially as the emitting 
region (~ 13kpc) and its HI column density is 10 19 cm -2 . The absorption line is 
blue-shifted 250km s _1 with respect to the emission peak. The 'screen' is kinemat- 
ically distinct from the emitting gas and therefore clearly external to the emitting 
clouds). This is a clear demonstration that HI screens — or mirrors depending 
on the perspective — exist on galaxy scales at these early epochs. From Fig. 5 
of Rottgering et al., we estimate that this cloud would reflect ~ 30% of the Lya. 
Thicker clouds could exist around other objects but would be hard to detect when 
only the faint wings at the extremities of the emission profile were transmitted. 
In the case of 0943-242, even if the screen contained dust, its effects would be 
negligible. The reason, as stated before, is that Lya photons (seen by the cloud) 
are incident from the exterior and will be very effectively scattered away before 
being absorbed by the dust. Even the photons getting through (i.e., without 
any scattering) to us because they are sufficiently far in the wings of the profile 
would not see much dust in this particular cloud. With such low column densities, 
r v ~ 5 10~ 2 ViVH = 7.5 10~ 5 for /i = 0.015 and t v ~ 0.005 for fi = 1. Much thicker 
clouds may result in some non-negligible extinction but again this does not arise 
because of the multiple scattering of Lya. 

5.4 Conclusions 

If the large-scale extended emission line regions (EELR) in radio galaxies are pho- 
toionized predominantly by the collimated UV radiation emitted by a hidden AGN, 
the line emission and transfer processes are characterized by what we call an 'open' 



Resonance scattering and dust in HZRGs 



100 



geometry with externally illuminated clouds. This produces an emission line spec- 
trum which is distinctly different from an internally ionized HII region, especially 
in the ultraviolet. 

Our principal conclusion is that internal dust, even in proportions as large as 
that found in the solar neighbourhood, does not satisfactorily explain the large 
CIV /Lja ratios commonly seen in groundbased observations of high redshift radio 
galaxies. We have demonstrated that geometrical effects can, in most cases without 
invoking any dust, explain the observed trends in the data. While internal dust 
can quite effectively kill the Lja luminosity, it is essential to examine also its 
effects on the other UV lines, especially CIV which is also a resonant transition. 
Although some internal dust within the ionized gas is allowed by our modeling 
and is required in a few cases, we find that the effects of geometrical perspective 
dominate the observed decrease of Lja without affecting excessively the other lines. 
This implies that the faintness of Lja cannot necessarily be used as an argument 
to demonstrate the existence of dust in these objects. 

The influence of these geometrical factors on the line spectrum gives us, po- 
tentially, a natural explanation for some of the emission line asymmetries seen in 
the 'alignment effect' radio galaxies. 

We have also investigated the effects that material external to the ionized 
cones will have on the appearance of the objects in the light of Lja. We suggest 
that the existence of diffuse Lja halos observed in some high z radio galaxies 
could be due to the reflection of Lya photons by neutral clumps lying outside 
the ionization cones. In addition to reflecting Lja, such neutral clouds could be 
responsible for the spatially extended absorption seen within the Lja emission 
profile. Pure Lja emission clouds need not, therefore, necessarily be identified as 
star-forming regions. 
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Chapter 6 
Discussion 



The projects developed in this thesis, although selfcontained, provide results rel- 
evant, not only to the individual class of objects or gaseous regions considered, 
but also to the nature of active galaxies in general. In this concluding chapter, I 
want to summarize what we have learned about the family of active galaxies and 
connect our results with the current ideas on this field. 

Active galaxies have been subdivided into several subgroups defined accord- 
ing to the strength and breadth of their emission lines and their luminosity over 
the entire electromagnetic spectrum. Sometimes, different physical conditions and 
emission mechanisms are responsible for the differences between subgroups. In 
some cases, however, objects pertaining to distinct groups can be intrinsically 
similar, the differences being due to pure geometrical effects. This statement sum- 
marizes the basic ideas of the so called "unified scheme models" (e.g., Lawrence & 
Elvis 1982, Lawrence 1987, Barthel 1989), mentioned several times in this work. 
Such model propose that different classes of AGNs are the same kind of object 
but seen from different directions. The directionality is taken to be due to an 
anisotropic distribution of absorbing or scattering matter and a simple geometry 
which would produce this situation involves an absorbing torus around the nucleus 
(the proposed "molecular torus"). In such scenario, the continuum produced in 
the central, hidden nucleus is emitted anisotropically into two wide, oppositely 
directed ionized cones, as a result of the collimation by the torus. 

Where does the unified scheme enter our discussion? It is represented in the 
kind of geometry assumed in all our calculations, in which a hard ionizing contin- 
uum illuminates the gas externally: a different situation from the one existing in 
HII regions, or planetary nebulae, where the continuum source (the stars) is em- 
bedded within the gaseous region. In AGNs, numerous pieces of evidence indicate 
that the main ionization mechanism of the EELRs is photoionization by a hard 
continuum, which is produced in the central engine (shocks could also play a role, 
see section 6.2). Such situation is better described by an open geometry. 

The resulting models are able to explain the observed data which suggests 
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that the scene is basically correct. However, the relatively high number of input 
parameters necessary to construct the models allows an apparent high degree of 
freedom to fit the observational results. One might argue that, using another 
geometry, the same strategy would have driven us to the same results. The question 
is, how unique is the scheme that can reproduce the observed data? 

Different observational pieces of evidence restrict the n-dimensional space of 
valid input parameters, to a much smaller space than that of all the models able 
to reproduce the observed emission line ratios. In this thesis, each chapter where 
theoretical models are presented, justifies in detail the assumed input parameters, 
based on our own or other authors' knowledge. We select the input parameters in 
accordance with the nature of the objects and regions under study. 

This chapter is divided into several sections, each one referred to a specific 
aspect of the nature of active galaxies: dust, ionization mechanisms and geometry. 
Every section is completed with a brief overview of the current state of the field, 
as well as a comparison between other authors' results and ours. 
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6.1 Dust in the ionized regions. 

Dust has been the main object of investigation in this thesis. We have provided 
a detailed analysis of the interaction between the dust grains and the radiation 
impinging on the gaseous component. We have also studied its effects on the 
line emission from the ionized gas like changes on the ionization structure of the 
gas, electronic temperature, metallicity, etc. Taking these effects into account is 
necessary in order to interpret more correctly what is observed. 

This understanding is useful, for instance, to distinguish when dust must be 
considered and when its effects can be neglected, allowing a simplification in the 
interpretation. Robinson et al. (1987, hereafter RBFT87) showed that in several 
diagnostic diagrams involving optical lines, the observed position of the EELR of 
radio galaxies at low z were well represented by a sequence in which one simply 
varies the ionization parameter, U (a measurement of the excitation level of the 
emitting gas, see Chapter 2). Using the observed values of the line ratios, they 
inferred that the EELR spectra were not reddened and presented models which 
did not include dust effects. 

I used the same set of data employed by these authors in order to investigate 
how dust mixed with the ionized gas would influence the predicted line ratios and, 
therefore, the positions within the diagnostic diagrams (Villar-Martm, 1993). The 
models showed that internal dust does not affect severely the optical line ratios 
(unless species susceptible of depletion are considered!). The influence is negligible 
for low U values and increases with this parameter. For highly excited EELR, like 
a few of the objects in RBFT87's sample, we should be more careful. At such 
high U, the effects of dust are stronger (see Chapter 2) and it might be necessary 
to consider them when interpreting their observed ratios. However, for typical 
EELRs at low redshift, where the ionization parameter is in the range [-4,-3] in 
the logarithm, dust effects on the optical line ratios are not relevant. The absence 
of reddening effects does not imply that dust is not present, instead is probably 
indicative of an open geometry for the gaseous ionized regions, where the clouds 
are externally illuminated. In a closed geometry, dust effects are stronger and 
the high reddening of some NLR line ratios of Seyfert 2 suggests a much higher 
covering factor than previously assumed (see Chapter 2). 

The situation described in the previous paragraph has advantages and disad- 
vantages. On the one hand, if dust effects are negligible, we can simplify enor- 
mously our modeling when interpreting the optical line ratios or the diagnostic 
diagrams that have been traditionally used to learn about the nature of the ioniz- 
ing source and the physical conditions of the ionized gas. We are therefore confident 
that our diagnostics are not wrong because of dust effects which we did not take 
into account. On the other hand, these same optical line ratios, easy to measure, 
do not tell us about the properties of the dust in the EELR, not even whether 
dust exists or not. As explained in Chapter 3, we need to develop other strategies 
and tools which are highly sensitive to dust, like observations of emission lines 
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emitted by species susceptible to depletion, as shown in detail in Chapters 3 and 4. 
In these chapters we have studied (theoretically + observationally) the emission 
of the [CaII]AA7291,7324 lines, which are appropriate because 1) calcium is an 
element very sensitive to depletion, and 2) as our models demonstrate, these lines 
should be strong and easy to detect under the physical conditions of EELR whose 
the gas does not contain dust. Their weakness demonstrates that dust does exist 
in the EELR of low z RGs. This favours the idea that the origin of the EELR 
is interaction or merging of two or more galaxies rich in gas, rather than cooling 
from a hotter corona, a scenario which is difficult to reconcile with the existence 
of dust mixed with the emitting gas (see section 4.4.1). 

6.1.1 Dust at high z. 

The subsequent question is: can we extrapolate our conclusion to high z EELRs? 
Can we affirm that they also contain dust? The answer is NO. The techniques 
applied in this thesis, do not demonstrate the existence of internal dust mixed 
with the ionized gas of high z radio galaxies. However, as explained in Chapter 
5 and below, there is strong evidence that dust exists in objects at redshifts as 
high as 2 and even more. Our work does not contradict this, but demonstrates 
that the line diagnosis based on UV (or optical) line ratios does not help us to 
infer the presence of dust because other factors, like geometry, dominate what 
we see. We have studied the UV emission line spectrum of a sample of high z 
radio galaxies analyzing the position in the diagnostic diagram CIVA1549/Lya vs. 
CIVA1550/CIII] A1909. Both line ratios are very sensitive to geometry and dust due 
to the resonance character of the CIVA1549 and Lya lines. We have demonstrated 
that dust, even in proportions as large as the one found in our ISM, can not 
explain the observed data, unless the influence of geometry id considered. In fact, 
our study shows that geometrical effects can in most cases, without invoking any 
dust, explain the observed data. Our models allow some internal dust but it must 
be in very little amounts (~ 2% of the dust/gas ratio in our ISM) in order to 
not contradict other line ratios, like CIV/CIII, which is not reddened. Effects of 
geometrical perspective, and not dust, dominate the observed line ratios. 

This can have important implications. As explained in Chapter 5, the weak- 
ness of Lya has been used as an argument to demonstrate the existence of dust at 
high z but we have shown that there is an alternative explanation based on purely 
geometrical effects. The resonance character of the Lya line makes neutral hydro- 
gen a very efficient mirror for the Lya photons which, depending on the geometry, 
prevent photons from reaching the observer. In this case the line becomes much 
fainter than expected while no photons have actually been destroyed but rather 
deflected from our line of sight. Conversely this same effect can be responsible for 
the apparent strength of Lya emitted from regions associated with neutral hydro- 
gen and which reflects Lya photons back towards the observer. This could explain 
the huge Lya halos observed around some objects, where the line originates from 
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zones which seem to be located out of the ionizing cones (like PKS0943-242, see 
section 5.3). This reflection effect is very efficient, whether dust is present or not, 
which means that Lja can even be stronger than predicted in the absence of dust, 
if the geometry is favourable. 

We have shown that the Lja line (its intensity or ratio with other emission 
lines) is not a reliable indicator of dust, but, on the contrary, a good tool to 
learn about the gas distribution (ionized as well as neutral) and, thus, about the 
geometry of the object. I will discuss further this point in section §6.4. 

The existence of dust at high z is indicated by other observational pieces of 
evidence like the extended polarized blue continuum mentioned several times in 
this work. Also, as explained in Chapter 4, there are several objects at high z 
detected in the millimeter spectral range, like 4C41.17 (z—3.8) and B2 0902+34 
(z—3.5) (Chini & Kriigel 1994; Dunlop et al. 1994). Their fluxes imply the presence 
of dust masses of the order of 10 8 M Q . Also high z quasars have been detected in 
the millimeter region by Andreani et al. (1993), McMahon et al. (1993) and Isaak 
et al. (1994). 

The paragraphs above refer to objects where the UV line emission is dominated 
by the EELR. There can also be objects where the unresolved nuclear emission 
dominates over the extended emission (EELR). In such cases, as demonstrated in 
Chapter 2, a rather closed geometry would be more appropriate. A good example 
is the IRAS source F10214+4724 (z=2.3) (Downes et al. 1992), a galaxy where the 
emission lines seem to reveal a very strong reddening. IR measurements indicate 
the existence of large quantities of dust, due in part to the amplification effects 
of a gravitational lens (e.g. Eisenhardt et al. 1995, Serjeant et al. 1995) which 
amplifies the IR luminosity. We know that dust is present, also from polarization 
measurements (Lawrence et al. 1993; Goodrich et al. 1996). But is its emission 
line spectrum as reddened as some authors claimed (e.g. Elston et al. 1994)? Our 
results indicate rather little reddening (or a very patchy distribution) while Elston 
et al. (1994) reported strong line reddening by measuring Ha/H/3 >20, which 
would imply Ay >5.4. However, according to these same authors, the CIV/CIII 
ratio is quite high (~ 3.6) (~3.4 according to Goodrich et al. 1996). Also, the 
HeII/CIII~1.7 and CIV/HeII~2.0 ratios (Goodrich et al. 1996) are similar to 
the values observed for other high z objects (see Chapter 5 and Villar-Martm, 
Binette & Fosbury 1995). Why is Hct/H/3 so reddened and not the UV line ratios?. 
Goodrich et al. (1996) for instance detect broad components of several UV lines 
in the polarized flux. For CIII, however, a broad component can be seen in direct 
upolarized light. Maybe Ha has similarly a broad component in direct unpolarized 
light while the fainter H/3 has none. This would result in an apparent higher Balmer 
decrement. By themselves, the UV line ratios do not indicate high quantities of 
dust and as explained in Chapter 5, if this object is a Seyfert 2 (Elston et al. 1994), 
a more closed geometry might be more appropriate for describing the emitting gas. 
And, as shown in Chapter 2, only small amounts of dust are required to strongly 
weaken the Lya emission. 
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In the near future, I want to combine what we have learned from our pho- 
toionization models, with the knowledge that the powerful polarization techniques 
have already provided. Much of what is known nowadays about the nature of the 
alignment effect observed at z >0.7 (see Chapter 4, Introduction) is the result 
of polarization studies, which distinguish between the scattered and the "in situ" 
continuum emission and tell us about the nature and distribution of the scattering 
material. In this way, for instance, di Serego Alighieri, Cimatti & Fosbury (1994) 
were able to exclude the possibility of scattering in two radio galaxies as due to 
hot electrons in a gaseous halo. 

Comparisons of our models with the observations do not allow a precise es- 
timation of the amounts of dust existing in the EELRs, both at high and low z. 
Studying the depletion of calcium, we have demonstrated that dust exists, but we 
can not say how much. Calcium is so sensitive to depletion that very tiny amounts 
of dust will deplet it completely: a dust /gas ratio (/i) =0.1 decreases the abun- 
dance of calcium in the gaseous phase by a factor of ~ 24 and if /i=l, by a factor 
of 5000! (see chapter 3). On the other hand, although polarization measurements 
set strong constraints to the nature and distribution of scattering material, it is 
still difficult to distinguish between cool electrons and dust as the origin of the 
scattering (Dey, Spinrad & Dickinson, 1995 and di Serego Alighieri et al. 1996). It 
could be interesting to investigate if the low dust/gas ratios which we deduce, are 
able to explain the level of polarization and the polarized flux observed in some 
objects. 

6.2 The ionization mechanism. 

There is strong evidence which support photoionization by the hard continuum 
emitted by the central AGN as the main ionization mechanism of the EELRs 
(and NLRs) of most radio galaxies and Seyferts, at least at low redshift. At 
higher z the uncertainties due to the lower spatial resolution and sensitivity of 
our instrumentation makes such study much more difficult. Moreover, at high 
z, most objects have been observed in the UV rest frame, a spectral range that 
is not well known, and still not often used for diagnosis. The optical line ratios 
that have provided so much information about the nature of low z objects, are 
redshifted into the IR where the observations are still difficult. The existing IR 
instrumentation (MIRACLE MPE/ROE camera, UKIRT CGS4 spectrometer, the 
3D Near Infrared Imaging Array Spectrometer (MPE), etc) has already allowed 
the study of the rest-frame optical wavelengths of a few objects at high z (>1) 
(e.g. Eales & Rawlings 1993, McCarthy, Elston & Eisenhardt 1992), but much 
more effort needs to be dedicated. 

Our models support that photoionization by the hard continuum produced in 
the central engine is the main ionization mechanism both for the NLR and the 
EELR line ratios of radio galaxies and Seyferts. However, there is an important 
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limitation of the optical line ratios as a diagnostic: with photoionization models 
it is possible to diagnose the mean ionizing photon energy (RBFT87), but not to 
distinguish between hard power laws and hot black bodies. Ionization by normal 
stars is ruled out, but a cluster of superhot stars, so called 'Warmers' could also 
explain the optical line ratios. 

At this point, I want to center this discussion around shocks, which have not 
been discussed yet in this thesis and deserve special attention. 

6.2.1 Shocks in active galaxies 

In many active galaxies radio jets have been detected. Very likely the energy 
which triggers them is originated in the central AGN. The jet contains relativistic 
particles accelerated to near the speed of light which propagate through the gas 
of the host galaxy. In its way out the jet will interact with the external gas and 
the shocks produced during such an interaction will disturb the kinematics and 
morphology of that gas. The hot gas behind a fast shock emits cooling radiation 
and this could also ionize the gas. However, as explained before, we expect that the 
central AGN plays also an important role in the ionization of the gas. But, which is 
the dominant process? Do shocks only alter the morphology and density of the gas 
or do they also modify its ionization state? Jets in AGNs are still surrounded by 
many mysteries and a clear understanding of the jet/cloud interaction phenomena 
is needed. 

At least for most low z objects, AGN photoionization has had better agree- 
ment with the observations. Beautiful evidence supporting this point came from 
the ionization cones imaged in several nearby Seyfert galaxies, like NGC5252 (Tad- 
hunter & Tsvetanov 1989; Acosta-Pulido et al. 1996; see also the clearly defined 
ionization cones in the WFPC images of NGC5728, Fig. 31). Wilson, Ward & Han- 
iff (1988) suggested a possible broad trend for the radio axis to be aligned with the 
axis of the extended emission line regions in Seyfert galaxies on a similar spatial 
scale. Wilson and Tsvetanov (1994) confirmed this tendency, pointing towards a 
common collimation mechanism of the radio plasma and the ionizing radiation. 

In many objects, spectroscopic evidence favours AGN against shock photoion- 
ization. For instance, the already mentioned trend defined by the EELRs line ratios 
in the diagnostic diagrams of RBFT87, which is well explained if they are ionized 
by the central AGN allowing for a different ionization state (i.e. U) in each object. 
It is difficult to understand this behaviour with the excitation parameter in the 
frame of the shock theory. 

On the other hand, there have been also problems with photoionization models 
that shocks were able to explain, at least partially. Three were the main discrep- 
ancies: a) too weak high excitation lines, b) too low electronic temperatures, and 
c) too small range in HeII/H/5. A recent work by Binette, Wilson & Storchi- 
Bergmann (1995) demonstrates that a mixture of matter- and radiation-bounded 
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Figure 31: 

An HST WFPC (1992) image of the core of the barred spiral Seyfert galaxy NGC5728. 
The bi-conical structure of the ionized gas indicates that ionizing continuum emitted in 
the central AGN is beamed in two opposite cones (Wilson et al. 1993). 
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clouds (proposed earlier by Viegas and Prieto, 1992; although not successfully) 
can solve these discrepancies, without the need of additional heating sources, like 
shocks. Moreover, from comparing U at different positions within the EELR in 
many RG (e.g., the radio galaxy 3C227, Prieto et al. 1993), the observations indi- 
cate that, in general, the ionization (U) does not change dramatically with radius, 
as we would expect if photoionization is due to the central AGN. This is a difficult 
question with no answer yet within the frame of pure AGN photoionization. [] 

There are some objects, where the presence of shocks is clear: a clear example 
is the intermediate redshift (z=0.310) radiogalaxy PKS 2250-41 which shows an 
emission line arc circumscribing the western radio lobe (Tadhunter et al. 1994), 
which strongly indicates the influence of the radio jet on the distribution and, pos- 
sibly also on the ionization of the warm gas. Also, the highly disturbed kinematics 
and the low ionization state near the radio structures is indicative of the presence 
of shocks. However, the source of ionizing energy is not clear: are the shocks dom- 
inating over the AGN continuum and ionizing the gas or are they just compressing 
the clouds in the arc, which are then ionized by the AGN? 

Although at low redshifts, the line emission spectra of most active galaxies are 
well explained by photoionization by the central AGN it is important to understand 
the jet/cloud interaction physics to disentangle on a case by case basis which effects 
are dominant. 

In our models we have also assumed photoionization by a hard continuum to 
explain the UV line ratios of high z radio galaxies. The trend in U shown by the 
data in the diagnostic diagram CIV/Lya vs. CIV/CIII, is indicative of this kind of 
excitation. But the possibility that shocks could also play an important role in the 
ionization of the gas can not be excluded. At high redshifts, radio galaxies show a 
strong alignment between the optical and the radio structures. One interpretation, 
well supported by polarization studies, has already been mentioned here: scatter- 
ing of continuum emitted by the hidden AGN. Another interpretation is that the 
propagation of the radio jet through the ISM of the host galaxy compresses the 
gas and triggers intense star formation (e.g. Rees 1987, de Young 1989, Begelman 
& Cioffe 1989). In fact the optical structures revealed by ground based and HST 
images (Longair, Best & Rottgering 1995; see Fig. 33) point towards shocked gas, 
with much narrower spatial distribution than the expected diffuse photoionized 
cone. High z radio galaxies have been selected on the basis of their powerful radio 
emission and line spectra. At such high z we are biased towards very powerful ob- 
jects, where the jet-cloud interaction might be very strong. Moreover, the shocked 
material, more concentrated, emits stronger lines than the surrounding gas. This 
could mean that due to the lack of sensitivity of our instrumentation, we are only 
detecting the material that is interacting with the radio plasma. 

In order to reach a clear understanding of the jet/cloud interaction physics, 

1 If shocks are strongly influencing the spatial distribution of the excited gas, even if they 
do not ionize it, the variation of U may be different than that expected from the pure AGN 
presence, as often assumed. 
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Figure 32: 

WFPC2 HST images of gaseous jets from three newly formed stars. 
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it is necessary to use low and intermediate z objects as laboratories, or even more 
local events, like the jets ejected from stars. HST WFPC2 images of jets originated 
by newly forming stars (e.g. Stapelfeldt et al. 1995) show the complicated inter- 
actions that take place when the ejected gas collides with the interstellar medium 
(see Fig. 32). The general picture which could describe this scenario is a new star 
surrounded by a dusty disk, which feeds material onto the star. When the star 
is hot enough, it will blow away much of the disk. The high speed material will 
collide with the slower gas in the interstellar medium and a bow shaped shock 
wave is created. In this way, a cavity is formed surrounding the star. If the shock 
wave encounters an obstacle (denser isolated clouds), a bow shaped shock wave is 
formed around the obstacle. Such a process will not only influence the ISM sur- 
rounding the star, but also the future of the star itself. The jet can, for instance, 
modify the mass accretion rate. 

The geometry in this scene is quite similar to that which describes the AGN 
phenomena, where the central source (presumably a black hole) is surrounded by 
an accretion disk that feeds it. What we learn from stellar jets can be very useful 
in understanding the role jets play in active galaxies. 

6.3 Geometry 

This thesis provides a tool for the understanding of the geometrical effects which 
influence the emission line spectrum radiated by an "open" geometry appropriate 
for the AGN phenomenon. This study demonstrates how important it is to take 
perspective effects into account; the narrow line emission depends on the orienta- 
tion of the ionized cones with respect to the observer. Such geometrical effects, 
studied in detail by us for the first time, have been used to explain some observa- 
tional effects, like the UV spectrum of high z radio galaxies and the discrepancies 
between the predicted and the observed HI line ratios emitted by the Narrow Line 
Region of Seyfert 2 galaxies (and, therefore, indirectly of Seyfert Is according to 
the unification schemes). As explained in section 6.1, we have studied how neutral 
material external to the ionized cones can influence the "appearance" of an object 
in the Lya light. 

The excited gas that we see does not necessarily map the distribution of ma- 
terial in the galaxy, but rather, the distribution of excited gas. Outside the radi- 
ation cones, there must be important reservoirs of neutral material. The work on 
NGC5252 by Prieto & Freudling (1993) shows that the neutral hydrogen surrounds 
the nucleus but is located outside the radiation cones of ionized gas. Also, in many 
high z radiogalaxies, the Lya profile shows absorption features due to absorption 
of Lya photons by intervening hydrogen gravitationally bound to the radio galaxy 
and spatially extended over a region of ~50 kpc in some cases (van Ojik 1995). 
The HI absorption systems might be due to neutral gas located outside the ionized 
cones. It is therefore very important to distinguish between excitation and matter 
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distribution. 

As explained in Chapter 5, the rest-frame UV spectral range constitutes a 
powerful tool for exploring the 3- dimensional structure of active galaxies in general: 
the existence of resonance lines in this spectral range which are so sensitive to 
geometry has made this study possible. In the "open geometry" that describes the 
EELRs, the two opposite cones, with apex in the nucleus which are ionized by the 
central AGN will most often be forming a non negligible angle with the plane of 
the sky in such a way that one of the cones will be closer to the observer. It is not 
a trivial task to distinguish which is closer and which is farther from the observer: 
only when there is some kind of observable asymmetry is this distinction possible. 
For instance, in high luminosity FRII sources, one-sided jets are detected, while 
two-sided jets are detected only in low luminosity (FR-I type) radio galaxies. This 
jet asymmetry can be explained by an illusion caused by a doppler-boosting effect 
due to the bulk relativistic motion of the jet. This effect makes the nearest side 
(to the observer) brighter than the other one. Moreover, in powerful radio galaxies 
with one-sided jets, the depolarization is systematically stronger on the counter-jet 
side (Laing 1988). Garrington & Conway (1991) interpreted this effect as produced 
by an external halo of ionized gas (identified by the authors with the extended 
X-ray emission often observed in clusters). The depolarization asymmetry is a 
geometrical effect, due to the different column densities of hot gas that the radio 
emission from both sides has to cross, the counter-jet side being seen through a 
longer path length of gas. 

When these asymmetries are present, we can readily distinguish which cone 
is closer to the observer and which one is further. However, it is not often easy to 
find these. For radio galaxies with no visible jets (or which are symmetric) or with 
the radio emission from both lobes too faint for the polarization to be measured, 
the distinction must be based on another strategy. 

We have proposed a new method to solve this problem using the strong emis- 
sion lines that are expected in the UV rest frame domain. Our models (Chapter 5) 
demonstrate that the line ratios involving resonance lines (like Lya and CIVA1550) 
should be different on both sides of the nucleus, with the differences increasing with 
the inclination angle with respect to the plane of the sky. This suggested us a test 
to distinguish between the cone closest to the observer and the furthest one: the 
cone which emits the spectrum resembling our "back" clouds will be the nearest 
one, while the furthest one should emit line ratios more similar to our "front" 
models. 

In this way, we expect to be able to describe the 3-dimensional structure of the 
objects. Furthermore, this will also help to understand some of the anisotropies 
that are observed in radio galaxies (some of them, in fact, expected from the 
unified schemes). For instance, the 3C256 radio galaxy (Dey, Spinrad & Dickinson 
1995) has different levels of polarization in different regions. The optical images 
obtained with the Keck 10m telescope, show two clumps of emission in the south 
and a fainter more diffuse component extending to the north west. Both regions 
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Figure 33: 

WFPC2 HST images of three high z radiogalaxies 3C265 (z =0.81), 3C324 (z =1.21), 
3C368 (z =1.13) (left to right). VLA data are superimposed on the images of 3C324 and 
3C368. The alignment between the radio axis and the UV rest-frame structures is clear 
in these two cases. The radio axis is marked by lines in the case of 3C265. The rest-frame 
UV continuum of this object has an elongated morphology misaligned by about 35° from 
the radio axis. 
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have rather different polarization levels. As the authors suggest, maybe there is a 
different contribution of diluting radiation. However, we have to be careful with 
this interpretation. It could also be a geometrical effect. The cross sections for 
forward and backward scattering of dust grains are different. This means that the 
dust scattering produced in both cones (when they are not in the plane of the 
sky) is asymmetric, and, therefore, the relative polarization measured will also be 
different. I plan to develop this study in the near future. 

McCarthy, van Breugel & Kapahi (1991) showed that in essentially all ra- 
dio galaxies the extended optical line emission is brightests at the cone closest to 
the nucleus. The most plausible explanation that the authors propose is inhomo- 
geneities in densities and distribution of dense clouds, which means that one of the 
cones is intrinsically brighter than the other. Assuming this to be the situation, 
in some cases, the brightest cone must be nearer the observer and in other cases, 
it will be further away. Therefore, (independently of which side is brighter) some 
objects can be expected to be dominated by a "back" type spectrum and others by 
a "front" type spectrum. We find that this effect could explain the high dispersion 
observed in the Lya/CIV ratio in very powerful radio galaxies at very high z, while 
CIV/CIII remains nearly constant from one object to another. 

As we see, there is evidence to support our view. A very strong test will be 
provided by spatially resolved UV spectroscopy of active galaxies with extended 
emission on both sides of the nucleus (objects devoid of any broad components 
will simplify enormously the study). The data that we have used for our modeling 
were extracted from the literature and correspond to high redshift radio galaxies 
where the spatial resolution is generally insufficient to separate clearly the extended 
from the nuclear emission. A project for the near future is to obtain long slit 
spectroscopy (UV rest frame) of high z radio galaxies with extended Lja already 
reported in the literature and where other strong UV lines (like CIV, III and Hell) 
have also been detected. The aim is to analyze the spatial distribution of the 
different lines and how the line ratios vary on both sides of the nucleus. So far, the 
spectra have generally been presented in the literature with no spatial information, 
compressing the objects in the spatial direction. Moreover, there has not been an 
attempt to obtain narrow band image in other UV lines but Lya. Quotient maps, 
like Ly«/CIV will describe the geometry of the ionized gas and will reveal regions 
of pure Lya emission probably indicative of the existence of neutral hydrogen. 

But, of course, the best tool for this research, available nowadays is the HST, 
that, with the Faint Object Spectrograph, makes possible UV spectroscopy at 
high spatial resolution for low z objects, where the distinction between the nuclear 
region and the extended ionized gas is much easier. This is one of the main goals 
of our HST project (section 6.6). 
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6.4 Our HST project 

I have mentioned the HST project in which I am involved and which got observing 
time for the next observing cycle. I will describe it in more detail in this section. 

We plan to obtain spatially resolved UV spectroscopy with the Faint Object 
Spectrograph (FOS) covering the observed frame Lya-CIII]A1909, of four low z 
radio galaxies with clear extended emission line regions already observed from the 
ground in the optical spectral range. With the data to be obtained, we will apply 
our newly developed ultraviolet line diagnostic and should be able to describe the 
3-dimensional geometry of the objects as well as study the dust effects on the mea- 
sured lines. We will be able to distinguish between the material on the near and 
far side of the AGN and compare our results with other observed asymmetries (see 
section 6.4). The main goal is, therefore, to elucidate the nature of the geomet- 
rical effects which influence the line spectrum radiated by the "open geometry" 
(externally ionized clouds) that describe the EELRs. 

In order to understand the nature of these objects, a clear separation be- 
tween the different components contributing to the spectral energy distribution is 
needed: AGN + stars + scattered continuum + additional sources (like nebular 
continuum). There are many phenomena that are a direct, or indirect consequence 
of the AGN activity, like the extended scattered continuum or the line emission 
(both if the ionization is directly by the AGN or indirectly by shocks induced by 
the radio jet crossing the ISM). The study of these phenomena will provide in- 
formation, not only on the AGN, but also on the stellar populations. Once we 
have a clear understanding of how the AGN influences its environment, it will be 
easier to separate the different components. What we learn will be very useful to 
interpret phenomena observed at very high z, where powerful radio galaxies are 
good laboratories to study the formation and early evolution of galaxies. So far, it 
has been very difficult to establish the state of evolution due to the difficulties in 
interpreting the observations with a non clear understanding of the AGN influence. 
If we knew the age of the oldest stars we could constrain the formation redshift 
and test even the current cosmological models. 

When we were choosing our targets we searched for objects with known 
anisotropics and for which we could distinguish which cone is on the near or far 
side. This will provide a strong constraint to test the validity of the diagnostic 
method that we have proposed. A practical limitation appeared when we realized 
that not many objects exist in which the aforementioned distinction between cones 
has been possible. We looked for targets with a clear one-sided jet, indicative of 
the cone nearer to the observer. However, objects containing jets are rather better 
candidates of classes strongly influenced by shocks produced by the interaction 
of the said jet with the cooler gas of the galaxy, while our models are more ap- 
propriate and intended for objects where photoionization by the central AGN is 
(more clearly) the main ionization mechanism of the gas. Therefore, we decided 
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to prepare a more complete proposal which would allow us the study of both phe- 
nomena: the influence of geometry and the jet/cloud interaction. We chose one 
target (PKS2 152-69) which presents clear evidence of interaction between the ra- 
dio and optical plasmas. Moreover, we selected one source (PKS0349-27) which 
show EELR on both sides of the nucleus and with no evidences of shocks. A 
new UV line diagnosis to distinguish between shock and AGN photoionization will 
be developed. With our newly developed diagnostic method, we expect to infer 
important clues on the geometry of the ionized gas. 

References 

Acosta-Pulido J. A., Vila-Vilaro B., Perez-Fournon L, Wilson A.S., Tsvetanov Z.I., 
Astrophys. J., 1996, in press 

Andreani P., La Franca F., Cristiani S., 1993, Mon. Not. Roy Astr. Soc, , 261, 
L35 

Barthel P.D., 1989, Astrophys. J., 336, 606 

Begelman M.C., Cioffe D., 1989, Astrophys. J., 345, L21 

Binette L., Wilson S.W., Storchi-Bergmann T., A&A, in press 

Cimatti A., di Serego Alighieri S., Fosbury R.A.E., Salvati M., Taylor D., 1993, 
Mon. Not. Roy Astr. Soc, 264, 421 

Cimatti A., di Serego Alighieri S., Field G.B., Fosbury R.A.E., 1994, Astrophys. 
J., 422, 562 

Chambers K.C., McCarthy P.J, 1990, Astrophys. J., 354 L9 
Chini R., Kriigel E., 1994, A&A, 288, L33 

Dey A., Spinrad H., Dickinson M., 1995, Astrophys. J., 440, 515 
de Young D.S., 1989, Astrophys. J., 342, L59 

di Serego Alighieri S., Cimatti A., Fosbury R.A.E., 1994, Astrophys. J., 431, 123 

di Serego Alighieri S., Cimatti A., Fosbury R.A.E., Perez-Fournon, 1996, Mon. 
Not. Roy Astr. Soc, submitted 

D'Odorico S., Cristiani S., Fontana A., Giallongo E., 1995, ESO PR 11/95; 15 
Sept. 

Downes D., Radford S.J.E., Greve A., Thum C, Solomon P.M., Wink J.E., 1992, 
Astrophys. J., , 398, L25 

Dunlop J.S., Hughes D.H., Rawlings S., Eales S.A., Ward M.J., 1994, Nature, 370, 
347 

Eales S.A., Rawlings S., 1993, Astrophys. J., 411, 67 

Eisenhardt P., Soifer B.T., Armus L., Hogg D., Neugebauer G., Werner M., 1995, 



Discussion. 



119 



BAAS, 186, 5202 

Elston R., McCarthy P.J., Eisenhardt P., Dickinson M., Spinrad H., Januzzi B.T., 
Maloney P., 1994, Astron. J. 107, 910 

Garrington S.T., Conway R.G., 1991, Mon. Not. Roy Astr. Soc, 250, 644 

Goodrich R.W., Miller J.S., Martel A., Cohen M.H., Tran H.D., Ogle P.M., Ver- 
meulen R.C., 1996, Astrophys. J., 456 L9 

Isaak K.G., McMahon R.G., Hills R.E., Withington S., 1994, Mon. Not. Roy Astr. 
Soc, , 269, L28 

Laing R.A., 1988, Nature, 331, 149 

Lawrence A., Elvis M., 1982, Astrophys. J., 256, 410 

Lawrence A., 1987, Pub. Astr. Soc. Pac, 99, 309 

Lawrence A. et al, 1993, Mon. Not. Roy Astr. Soc, 260, 28 

Lilly S.J., 1988, Astrophys. J., 333, 161 

Lilly S.J., 1989, Astrophys. J., 340, 77 

Longair M.S., Best P.N., Rottgering H.J.A., 1995, Mon. Not. Roy Astr. Soc, 275 
L47 

McCarthy P.J., van Breugel W.J.M., Kapahi V.K., 1991, Astrophys. J., 371, 478 
McCarthy P.J., Elston R., Eisenhardt P., 1992, Astrophys. J., 387, L29 
McCarthy P.J., 1993, ARA&A 31, 639 

Macchetto F.D., Lipari S., Giavalisco M., Turnshek D.A., Sparks W.B., 1993, 
Astrophys. J., , 404, 511 

McMahon R.G., Omont A., Bergeron J., Kreysa E., Haslam C.G.T., 1993, Mon. 
Not. Roy Astr. Soc, 267, L9 

Prieto A., Walsh J., Fosbruy R.A.E., di Serego Alighieri S., 1993, Mon. Not. Roy 
Astr. Soc, 263, lOp 

Prieto A., Freudling W., 1993, Astrophys. J., 416 668 
Rees M.J., 1987, Mon. Not. Roy Astr. Soc, 239, lp 

Robinson A., Binette L., Fosbury R.A.E., Tadhunter C.N., 1987, Mon. Not. Roy 
Astr. Soc, 227, 97 (RBFT87) 

Serjeant S., Lacy M., Rawlings L.J., King L.J., Clements D.L., 1995, Mon. Not. 
Roy Astr. Soc, 276, L31 

Stapelfeldt K.R. et al, 1995, BAAS, 185, 1802 

Tadhunter C, Tsvetanov Z., 1989, Nature, 341, 422 

Tadhunter C, Shaw M., Clark N., Morganti R., 1994, A&A, 288, L21 

van Ojik R., 1995, Ph.D. thesis, Rijksuniversiteit te Leiden 



Discussion. 



120 



Viegas S.M., Prieto A., 1992, Mon. Not. Roy Astr. Soc, 258, 483 

Villar-Martm, 1993, Proyecto de tesis: "Influencia del polvo en modelos de fotoion- 
izacion aplicados a niicleos activos de Galaxias" . Ito. de Astrofisica de Canarias y 
Universidad de La Laguna. 

Villar-Martm M., Binette L. & Fosbury R.A.E, 1995, in Proc. for the Workshop 
" Cold Gas at high z" , Hoogeveen, August 1995 (in press) 

Wilson A.S., Ward M.J., Haniff C.A., 1988, Astrophys. J., 334, 121 

Wilson A.S., Braatz J.A., Heckman H.M., Krolik J.H., Miley G.K., 1993, Astro- 
phys. J., , 419, 61 

Wilson A.S., Tsvetanov Z., 1994, Astron. J. , 107, 1227 



Chapter 7 
Conclusions 



We have considered for the first time the combined effects of geometrical per- 
spective and dust on the narrow emission line spectra of active galaxies. This is 
necessary to make a more correct interpretation of the observations. The kind 
of geometry determines strong differences on the emission line ratios when dust is 
present, due to the much more destructive effects of dust when we approach a close 
(spheric) geometry. Moreover, when we study the line emission from the EELRs or 
objects dominated by such extended emission, the influence of perspective must be 
considered if resonant lines are involved, both if dust is or is not present. In gen- 
eral the ionized cones have an inclination angle with respect the plane of the sky, 
which results on different viewing perspectives for both cones. This implies strong 
asymmetries on the UV line ratios of regions at both sides of the nucleus, due to 
the resonant character of Lya and CIV and therefore, their strong dependence on 
perspective. 

We have demonstrated that the NLR HI line ratios of Seyfert 2 galaxies points 
towards a significantly self-covered distribution of gas mixed with dust which is 
very patchy. The combination of line scattering and dust absorption effects can 
explain the observed ratios. Therefore, for those objects where the line emission 
is dominated by the unresolved Narrow Line Region, a more close geometry than 
peviously assumed is appropriate. On the other hand, when the extended emission 
is dominant (EELR), an open geometry is more realistic. If such close geometry 
is valid for the NLR, we expect the central continuum source to be also reddened. 
We claim that the apparent deficit of ionizing photons seen in many Seyfert 2's 
could in part be due to dust extinction of the ionizing continuum, and not only 
to the presence of a circumnuclear torus that collimates the ionization radiation 
along its axis, as proposed by other authors. 

Our studies on the UV emission line spectra of high z radio galaxies demon- 
strate that the UV line ratios are strongly determined by geometrical perspective, 
rather than dust effects. Internal dust, even in proportions as large as that found 
in the solar neighbourhood, can not explain the UV line ratios measured for high 
z radio galaxies, unless perspective effects are taken into account. 
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We claim that Lya is not a good indicator of dust but rather a good tool 
to understand the geometry of the emitting gas. Moreover, the influence of geo- 
metrical factors could give a natural explanation for some of the emission line 
asymmetries seen in the 'alignment effect' radio galaxies. The reflection of Lja 
photons by neutral material can determine the morphology of the object in the 
light of Lja. We suggest that the diffuse Lja halos observed in some high z radio 
galaxies could be due to the reflection of Lja photons by clumps lying outside the 
ionization cones. 

We have demonstrated that the gas in Extended Emission Line Regions in low 
z radio galaxies is mixed with dust. The faintness of the [CaII]AA7291,7324 lines 
indicates that calcium is depleted onto dust grains. We have used this result to 
discriminate among the current theories about the origin of the extended ionized 
gas. If there is an universal origin, the existence of internal dust favours debris 
from mergers or tidal interactions as the most plausible scenario. 
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Apendix A 

Binette, Wang, Villar- Martm, Martin & Magris 1993, ApJ, 414, 535 



The term depletion refers to the underabundance of gas phase elements with 
respect to the solar standard abundances as a result of their assumed presence 
in dust. References to the studies of depletion in the local ISM can be found 
in Whittet (1992). We describe here a theoretical algorithm which attempts to 
deplete in a meaningful way the trace element abundances when computing line 
diagnostic for a plasma of arbitrary gas metallicity Z gas and dust content /i. Both 
quantities Z gas and fi (oc Zdust] see eq. [A5]) are known to depart significantly from 
solar neighborhood values in low-mass extragalactic systems or as a function of 
galactic radius in the Galaxy. For the purpose of modeling plasma conditions in 
a nonsolar environment, we defined an algorithm which scales in a self-consistent 
manner the depletion of the gas phase in proportion to the quotient [i/Z, where Z 
is the total metallicity of the plasma [Z = Z dust + Z gas ) expressed relative to local 
ISM values: 

ry _ Y,m X {{n X }dust + {nx}gas)pl asma r ^ 1 1 

and ji is the dust content of the plasma expressed relative to that of the 
standard ISM dust-to-gas mass ratio. The reference solar abundances set which 
we adopt is taken from the compilation of Anders & Greveese (1989) and is listed in 
Table 4 (see entry {nx}/{%} ). We define D(X), the depletion index of element 

X, as 

D(X) = £ ] = log{^} IS M - W{— }o [A2] 
H uhL riff 

where {nxL/nnL}isM represents the ratio of the column densities of the 
gaseous element X (measured along various lines of sight of path length L[< 6 x 10 20 
cm]) over the corresponding measured column densities of H. In the local ISM, the 
number density of atoms % of element X which are found either in gas form 
(or depleted into dust grains) is given relative to n# by /n^} 10°^ (or 
{nx/nf/}©(l - 10 D ( X )). The adopted set of local ISM depletion indices D(X) is 
listed in Table Al and was derived from Table 2.2 of Whittet (1992) except for 
carbon for which the value -0.5 was adopted as suggested by some authors. This 
value is also consistent with our dust model which includes graphite grains. 

For the local ISM, for which ji = 1 (and Z = Zdust + Z gas =1), the gas phase 
abundances are derived directly from the depletion indices of Table Al. As we are 
often interested in modeling line emission plasma with Z ^ 1, we could assume 
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Table 4: Solar Abundances, ISM Depletion Indices and Correlation Coefficient 
Ratios. 



Element 


{n x /n H }o 


D(X)ism 


{r n /r N y x tiM 


H 


1.000 






He 


0.1 


0.0 




C 


3.6 x 10~ 4 


-0.5 


0.88 


N 


1.1 x 10~ 4 


-0.1 


0.46 





8.5 x 10" 4 


-0.2 


0.42 


Ne 


1.2 x 10~ 4 


0.0 




Mg 


3.8 x 10~ 5 


-0.7 


1.11 


Si 


3.5 x 10~ 5 


-1.6 


1.19 


S 


1.9 x 10- 5 


-0.2 


0.94 


Ar 


3.6 x 10~ 6 


0.0 




Ca 


2.2 x 10~ 6 


-3.7 


1.42 



that the dust content fj, scales linearly with Z. That this is approximately true is 
supported by the results of Issa, MacLaren, & Wolfendale (1990) and Sauvage & 
Vigroux (1991) who compared the metallicity and extinction in the Galaxy with 
those of the Magellanic Clouds or of nearby spirals. In our calculations, if we 
adjusted the dust content such that always fi/Z = 1, assuming that the depletion 
indices of Table Al remain approximately valid, the procedure for depletion is 
straightforward and the gas phase abundance by number of element X is simply 
given by 

{n x /n H } gas = Z{n x /n H } Q 10 D{x) , [A3] 

This is equivalent to saying that when /i/Z =1 both the dust grains as well 
as the gas phase are made up of the same relative proportion of trace elements. 

For the more general case where \i ^ Z, we have no observational basis as to 
how the depletion indices might vary. If we use the local ISM indices as a guess 
of the depleted abundances, and then simply renormalize the whole abundance set 
by a constant factor until the depleted mass match the value implied by \x (see. 
eq. [A5] below), we obtain, in the case fi <^ Z, rather implausible gas abundances 
for some elements. To alleviate this problem, we have chosen instead to vary the 
depletion indices in some proportion to the ratio \ijZ before linearly renormalizing 
the total mass of the depleted metals. An indication of the relative variation of 
the depletion indices D(X) with the ratio \xjZ can be intuitively inferred from the 
measured quatities \r n jr x Yx M which represents the ratio for element X of the 
correlation coefficient of D(X) with density, r n , over the correlation coefficient of 
D(X) with column density, r N . The merits of the measured {r n /r N y^ M ratios 
is that they sample regions of somewhat different dust content /i. The algorithm 
we devised for each case fi ^ Z is the following: we first define modified depletion 
indices D(X, fi/Z) 1 which are dependent on ji/Z as follows 
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D(X,p/Z)' = D(X)(^y {r " /rN}x , [AA] 

By trial and error, we established that 7 = 0.25 gave a reasonable behaviour of 
the depleted elements even when p/Z — > 0. The values we adopted for {r n /r x Yx M 
(see Table 4) were derived from Figure 2.9 of Whittet (1992). Linear renormaliza- 
tion of the abundances implied by the modified indices D(X, p/Z') is required in 
order that the mass of the elements depleted onto dust grain remains consistent 
with the value of p. This is done in an iterative fashion by determining the value 
of 5 (as 1) until the following dust mass equation is satisfied 

Z dust E m x {n x } Q = 0.0089/i; " 1 " " 1 



J2mx{n x } 



ZV{^} [1 - K min(l,510 D ^ z '^], [A5] 
^ m H n H 

where the constant 0.0089 represents the dust-to-gas mass ratio relative to 
H (i.e., pdust/ Ph) of any of the extinction curves computed by P.G. Martin (this 
number is consistent with the ISM value if our dust composition within the grain 
model is valid). The constant k = 1.123 is a small correcting factor to allow for 
the fact that the summation is limited to the elements of Table 4 (considered by 
MAPPINGS) which does not include all the depleted elements (e.g., Fe). Once 5 
is determined, the gas phase abundance of any element X is 



{nxMZ = Z Ux/n H } e mm(l, 510 D ^ Z ">), [A6] 

The value of k is such that 5 becomes unity if p/Z = 1 in which case equation 
(A5) reduces to equation (A3). With the above depletion algorithm, as p/Z — > 0, 
the uncertainty in the absolute abundances of the highly depleted elements (e.g., 
Ca) is probably quite high but since no observational depletion data for p/Z <^ 1 is 
available, this situation cannot be helped. The main purpose of our algorithm was 
that it gave a cummulative depletion of the gas phase elements always consistent 
with the dus content p. Furtermore, as p/Z — > it also has the property of giving 
for all the elements cocerned much more plausible values than any alternative linear 
scheme which we could devise. The need for a depletion algorithm even if imperfect 
is called for in any photoionized gas model since it is quite conceibable that the 
dust content is lower than in the ISM as a result of destruction of dust grains by 
sputtering (although the timescales are relatively long compared to the lifetime of 
an HII region; see Osterbrock 1989) or as a result of segregative acceleration away 
from the cloud of much of the smaller grains' population due to radiation pressure. 
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1. CALCULATIONS OF N%f° 

We assume a slab geometry for the screen. For an isothermal screen with HI 
column density N^j een , the optical depth for Lyct resonant scattering is 

r(x) = N s H r en MT)<t>(x), [Bl] 

where x = {uj—ujq)uj D) <t (T)4>(x) is the thermally averaged resonant scattering 
cross section, 

a (T) = 5.9 x 10- 14 (T/10 4 iir)- 1 /2 cm" 2 , [B2] 
a f+ QO e~ y2 dy 



7T J-oo a 2 + [x — v) 



-oo a 2 + (x — y) 2 

a = T/2cu D = 4.7 x 10" 4 (T/10 4 K)" 1/2 is the damping constant with T" 1 the 
radiative lifetime of the excited state, and lud = ujq{2kbT /mHC 2 )^ 1 ^ 2 is the thermal 
Doppler width with ujq/2ti) = 2.47 x 10 15 Hz the line center (Lya frequency 
(Bonilha et al. 1979). For photons of frequency Ui and direction cosine rij injected 
into a screen moving at velocity v, we evaluate the optical depth (a Lorentz scalar) 
in the cloud (comoving) frame to obtain 

r(uj doud ) =N£ em a (T)<f>(u cloud ,T), [B4] 

where io doud = 7^(1 — v.n) and 7 = (1 — v 2 )" 1 / 2 . For line center injection, 
Ui = u , and the column density at which r{uj c i oud =1, N^}*, is given by 

Nfff = . . } = Nfffiy). [B5] 

a (T)(/)(uj doud ,T) 

2. CALCULATIONS OF / 

Let N(x) be the flux of Lya photons impinging on a screen which we model as 
a plane parallel slab. The width of the impinging line feature is 5V S (in km s _1 so 
that x = (uj — uq)/(uj5V s ). Let the width of the Lya resonant scattering line profile 
in the screen be given by SV a (in km s -1 ). The probability that a line photon will 
suffer an interaction when traversing the screen is 1 — e~ T ^ x '\ where r(x') is given 
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by equation (Bl), and x' = (ou — uj )/ (uj 5V a ) = kx with k = V s /5V a . The fraction 
of Lya photons that suffer interactions when traversing the screen is then 

f - ■l'-^ dxN ^ 1 ~ e T(X) ] _ f( N screen \ [n.l 

/_ J^dxN(x) - f ^ m ' K) - [m 

Taking N(x) oc e~ x2 gives 

/ = n^ 2 f +Q ° dxe' x2 [1 - e~ T{x \ [Bl] 

Equation (B6) also gives the fraction of the impinging Lya energy flux that 
interacts with the HI gas in the screen since the impinging energy flux of Lya is 
approximately huj N(x). 
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The transfer of the UV continuum across an arbitrary thick slab of dust which 
is discussed in 52.4.4 was carried out using the numerical solution described in 
detail by Magris (1985; see also Bruzual et al. 1988). Let us consider an arbitrary 
frequency with its corresponding absorption (a a b s ) and scattering (cr sca ) dust cross 
sections. We are interested in solving the following transfer equation 

p*l!^ = I(T,ri-S 8 (T,ri, [CI] 

where dr = (cr sca + a abs )dx is the differential optical depth of the medium. 
This equation which takes into account scattering as well as absorption by homo- 
geneously distributed dust particles is appropriate to the plane parallel geometry 
case. The scattering source function due to the dust grains is given by 



S s = ^j I{t^)-S s {t,h), [C2] 



where a = a sca /{a sca + a a bs) is the albedo and <E>(#) = (1 — g 2 )/(l + g 2 — 
2g cos9) 3 / 2 is the scattering phase function. The phase function is characterized 
by the asymmetry parameter g =< cosO > with 9 the scattering angle and g the 
Henyey-Greenstein asymmetry parameter characterizing our dust grain models and 
tabulated by P.G.Martin as a function of frequency along with a abs and a sca (see 
grain model in Martin & Rouleau 1991). The transfer equation was solved by using 
the finite difference solution method of Milkey, Shine, & Mihalas (1975) with the 
following boundary conditions: 

I(t,h) = for r = 0, fj, < 

[C3] 

J(t, fi) = I c for r = T exU fi > 0, 



where T ext is the optical depth of the slab and I c is the intensity of the incident 
continuum radiation which illuminates isotropically one side of the dust slab. To 
obtain the numerical solution, we use a logarithmically uniform mesh that is sym- 
metric about slab points. In order to test the accuracy of our numerical solution, 
we have compared our results with the analytical solution of Roberge (1983). In 
particular, we could reproduce very well the behavior of the mean intensity versus 
optical depth inside the slab when illuminated isotropically on both surfaces. 
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By comparing the estimated photoionization rates from the excited level 3d 
of Ca + *, n 3 d, to that from the ground state 4s of Ca + , n 4s , we show that pho- 
toionization of excited Ca + * is a negligible process. 

A fundamental parameter which determines the importance of the ionization 
rate from the metastable level is its relative population with respect to the ground 
state: being n 4s the density of Ca + ions in the ground level and n 3d the density 
of Ca + ions in the metastable level, 3d. To evaluate this fraction, we have solved 
analytically the statistical equilibrium equations. To simplify the calculations we 
have considered a three level atom, reducing the two 4p sublevels (see Fig. 16, 
Chapter 3) to a single level and the same for the 3p sublevels. We have taken 
into account all the processes (collisional and radiative) which can populate or 
depopulate each of the levels. 

The resolution of the system of three equations gives us the ratios ■^ L -, with 
L=3d,4s,4p, being n Ca + the total density of Ca + ions. The density and temperature 
we considered were 300cm -3 and 10000K, respectively. From this we deduced the 
relative population with respect the ground state. The results turn out to be: 

^ ~ io- 7 

n As 

and 

n 4s 

The negligible population of the upper 4p level prevents any contribution by 
cascade to the population of the 3d level, therefore collisional excitation is the only 
important mechanism populating the metastable 3d level. This is consistent with 
the fact that we do not observe the triplet of Ca+ (8498,8542,8662) (4p to 4s) in 
any EELR although it is observed in the broad line region of AGN where densities 
are higher by many order of magnitudes. 

a) Ionization of Ca + * by soft continuum photons. 

The soft UV counterpart of the ionizing continuum provides a source of ion- 
ization for both Ca+ (IP: 11.9eV) and Ca+* (IP: 10.2eV). To estimate the pho- 
toionization rates, we will make the following approximations: 

1) At the fairly large depth in the cloud where the specie Ca + becomes abundant, 
we only need to consider photons with energies < 13.6eV, the ionization potential 
(IP) of H°, because photons just above this energy have already been absorbed 
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and also because of the rapid decrease of the photoionization cross section with 
increasing energies. 

2) The ionizing continuum which reaches the PIZ is considered to be the soft UV 
counterpart of an ionizing PL of index —1.4 (but unattenuated since the opacity 
due to dust or trace elements below 13.6eV is relatively small). The continuum 
impinging the cloud is described by F v = F s v~ 1A . In number of photons this is 
F u /hu = F s /hv~ 2A . The constants F s /h will cancel out when taking the ratio 

rwn 4s . 

The atomic data was taken from Osterbrock (1987) for H° and from Shine & 
Linsky (1975) for Ca + and Ca + *. Tables 1, 2, 3 where we define a' v = a v * 10 18 
and v' = z//10 16 show the relevant atomic data. The threshold ionizing frequency 
of H° is labeled u . 

We estimate the quotient n 3rf /n 4s as follows: 

„ / u- 2A a„(Ca+*)dv 

^3d _ nzdyzd 



H4s n4s J v~ 2A a u (Ca + )dv 



"4s 



If we define a' u = a u * 10 18 and v' = z//10 16 , we have 



/ v'- 2A a' u {Ca + *)dv' 

II4C U4 S v 

J u'- 2A a' v {Ca+)dp' 

The values of a' u (Ca + *) at v'^v' u and a' v {Ca + ) at u' and u' As have been obtained 
by interpolation from (see Tables 2 and 3). Using the rule of the rectangle to 
approximate a given integral: 

XI 

J f(x)dx ~ (xi - x )f(x ) + (x 2 - xi)f(xi) + ... 
xo 

)f(Xn-l) 



[D2] reduces to 

U 3d n 3d 10.06 



n 4s n 4s 0.156' 



[£>3] 



If we take into account that ^ ~ 10 7 we obtain 

714s 
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Table 5: H° photoionization cross section 



v' 




v'- 2A a' u (H°) 


0.3 


6.4 


115.1 


0.4 


2.5 


22.54 


0.6 


1.15 


3.919 


0.8 


0.6 


1.025 


1 


0.3 


0.3 


1.2 


0.2 


0.1291 


1.4 


0.15 


0.06689 


1.6 


0.1 


0.03237 


1.8 


0.07 


0.01708 


2 


0.05 


0.009473 


2.2 


0.04 


0.006029 


2.4 


0.03 


0.00367 


2.6 


0.02 


0.002019 


2.8 


0.015 


0.001267 


3 


0.015 


0.001074 


3.2 


0.01 


0.0006133 


3.4 


0.01 


0.0005302 


3.6 


0.01 


0.0004622 




n 4s 


6 x 10~ 6 



which means that the ionization of Ca + into Ca ++ originates overwhelmingly from 
the ground level 4s. 

b) Ionization of Ca + * by Lya photons. 

Lya photons have energies slightly higher than the treshold energy of level 
3d of Ca + * and constitute undoubtedly an important source of ionization for this 
level. To compute the Lya emissivity, we assume the low density regime whereby 
— | of recombinations of H + lead to the emission of a Lya photon. Since Lya is a 
resonance line of large line scattering opacity, we consider that the density of Lya 
photons within the Ca + region result from the production of Lya either locally 
or from deeper regions. The justification for this is that resonance line photons 
generated from layers nearer the slab'surface would be reflected outward as a result 
of the increasing fraction of H° with depth (see Binette et al 1993b). Therefore the 
Lya flux potentially available to ionize level 3d is a small fraction rj of the total 
Lya flux emitted by the cloud. This fraction is of order 0.2 corresponding to the 
fraction of ionizing photons of H° not yet absorbed at the typical depth where the 
Ca + specie is abundant. 

We first compute the n^P/IL^ ratio taking into account that resonance scat- 
tering will increase the density of locally emitted line photons by a factor £ ~ 10 7 , 
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Table 6: Ca + 4s photoionization cross section 



\{A) 


v' a' v (Ca+) v'- 2A a' u (Ca + ) 


1044 


0.2873 (z4) 


0.2036 


4.063 


1000 


0.3 


0.2097 


3.772 


950 


0.3158 


0.2145 


3.412 




0.3288 (i/ ) 


0.2157 


3.114 


900 


0.3333 


0.2170 


3.031 


850 


0.3529 


0.2172 


2.644 


800 


0.375 


0.2149 


2.262 


750 


0.4 


0.2103 


1.896 


700 


0.4286 


0.2033 


1.554 


650 


0.4615 


0.1942 


1.242 


600 


0.5 


0.1830 


0.966 


550 


0.5455 


0.1700 


0.7282 


500 


0.6 


0.1554 


0.5295 


450 


0.6667 


0.1394 


0.369 


400 


0.75 


0.1225 


0.2443 


350 


0.8571 


0.1049 


0.1518 


300 


1 


0.0870 


0.08697 


250 


1.2 


0.0691 


0.0440 


200 


1.5 


0.0519 


0.0199 


150 


2 


0.0352 


0.00668 


100 


3 


0.0199 


0.001417 


50 


6 


0.0198 


0.0002685 
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Table 7: Ca + * 3d photoionization cross section 



A(A) 


v' 


a! v {Ca + *) v'- 2A a' u (Ca+*) 


1218 


0.2462 ( 


3dJ 


6.148 


177.70 


1200 


0.25 




6.086 


169.50 


1150 


0.2609 




5.907 


148.60 


1100 


0.2727 




5.716 


129.20 


1050 


0.2857 




5.511 


111.4 


1044 


0.2873 ( 


4s J 


5.403 


107.80 


1000 


0.3 




5.295 


95.22 


950 


0.3158 




5.066 


80.56 




0.3288 ( 




4 946 5 


71 39 


900 


3333 

V / • '.J 'J 'J 'J 




4.827 


67.41 


850 


0.3529 




4.576 


55.72 


800 


0.375 




4.315 


45.43 


750 


0.4 




4.044 


36.46 


700 


0.4286 




3.762 


28.74 


650 


0.4615 




3.47 


22.19 


600 


0.5 




3.168 


16.72 


550 


0.5455 




2.856 


12.23 


500 


0.6 




2.534 


8.634 


450 


0.6667 




2.202 


5.827 


400 


0.75 




1.862 


3.714 


350 


0.8571 




1.517 


2.196 


300 


1 




1.173 


1.173 


250 


1.2 




0.8424 


0.5439 


200 


1.5 




0.5401 


0.2041 


150 


2 




0.2865 


0.05428 


100 


3 




0.1051 


0.007524 


50 


6 




0.01471 


0.0001995 
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the mean number of scatterings before escape. Such a high number of line scatter- 
ing, however, characterizes only the locally produced Lya photons. This accumu- 
lation (or slowing down) effect which we want to estimate is only effective within 
a zone of optical depth of order a few. Let's take r scat ~ 10 = <j scat n(H ) 5X 
where 5X is the geometrical depth and n(H°) the local density of H°. Adopting 
a sca t ~ 6 x 1CT 14 cm 2 (cf Appendix B of Binette et al 1993b), the column den- 
sity of recombining H + we should consider in generating Lya is within a thickness 
N(H+) ~ N(H°) = n(H°) 5X = 10/6 x 10~ 14 = 1.7 x 10 14 cm" 2 . 

Taking these considerations into account, the problem reduces to estimating 
the quotient 



ivt aLya (Ca+*)N(H+) J v'-*X(H°W 



u 3d n 3d 



His fife u 'o 

J v'-^a' v {Ca+)dv' 

= —\^a Lya {Ca + *)N{H + )^ [DA] 

where aL ya (C& + *) is the photoionization cross section of Ca + * at the Lya energy 
(~6x 10 _18 cm 2 ). Taking again n 3 rf/n 4s ~ 10 -7 , we obtain 

^ 1.36 x 10" 4 -, [D5] 
Using again the rule of the rectangle we obtain 

^~1.5x 10- 2 . 

The effect of slowing down of resonance line photons is interesting but appears 
insufficient as it involves too small a fraction of the Lja photons generated within 
the nebula. Let's now estimate the importance of all the Lya photons generated 
within the deeper zones which, after having scattered far enough in frequency to 
escape, must still cross the Ca + zone (without appreciable scattering in that zone). 
The effect on the photoionization of level 3d is given in this case by 



Il4 s U4 S v 

J v>- 2A a' v {Ca+)dv> 



, [D6] 



n 4s 3 I 2 

(Note that a' is used here for Ca + *). Using the rule of the rectangle we obtain 
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M - 2.43 x 10" 6 



n 



r-NJ 



■Is 



In summary, Lya emission and/or trapping are insufficient to ionize Ca + . As 
in the case of soft continuum photons, this is basically the result of the extremely 
small population characterizing the excited level 3d. 
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